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The  National  Bureau  of  Standards'  was  established  by  an  act  of  Congress  on  March  3,  1901. 
The  Bureau's  overall  goal  is  to  strengthen  and  advance  the  Nation's  science  and  technology 
and  facilitate  their  effective  application  for  public  benefit.  To  this  end,  the  Bureau  conducts 
research  and  provides:  (1)  a  basis  for  the  Nation's  physical  measurement  system,  (2)  scientific 
and  technological  services  for  industry  and  government,  (3)  a  technical  basis  for  equity  in 
trade,  and  (4)  technical  services  to  promote  public  safety.  The  Bureau's  technical  work  is  per- 
formed by  the  National  Measurement  Laboratory,  the  National  Engineering  Laboratory,  and 
the  Institute  for  Computer  Sciences  and  Technology. 

THE  NATIONAL  MEASUREMENT  LABORATORY  provides  the  national  system  of 
physical  and  chemical  and  materials  measurement;  coordinates  the  system  with  measurement 
systems  of  other  nations  and  furnishes  essential  services  leading  to  accurate  and  uniform 
physical  and  chemical  measurement  throughout  the  Nation's  scientific  community,  industry, 
and  commerce;  conducts  materials  research  leading  to  improved  methods  of  measurement, 
standards,  and  data  on  the  properties  of  materials  needed  by  industry,  commerce,  educational 
institutions,  and  Government;  provides  advisory  and  research  services  to  other  Government 
agencies;  develops,  produces,  and  distributes  Standard  Reference  Materials;  and  provides 
calibration  services.  The  Laboratory  consists  of  the  following  centers: 

Absolute  Physical  Quantities2  —  Radiation  Research  —  Chemical  Physics  — 
Analytical  Chemistry  —  Materials  Science 

THE  NATIONAL  ENGINEERING  LABORATORY  provides  technology  and  technical  ser- 
vices to  the  public  and  private  sectors  to  address  national  needs  and  to  solve  national 
problems;  conducts  research  in  engineering  and  applied  science  in  support  of  these  efforts; 
builds  and  maintains  competence  in  the  necessary  disciplines  required  to  carry  out  this 
research  and  technical  service;  develops  engineering  data  and  measurement  capabilities; 
provides  engineering  measurement  traceability  services;  develops  test  methods  and  proposes 
engineering  standards  and  code  changes;  develops  and  proposes  new  engineering  practices; 
and  develops  and  improves  mechanisms  to  transfer  results  of  its  research  to  the  ultimate  user. 
The  Laboratory  consists  of  the  following  centers: 

Applied  Mathematics  —  Electronics  and  Electrical  Engineering2  —  Manufacturing 
Engineering  —  Building  Technology  —  Fire  Research  —  Chemical  Engineering2 

THE  INSTITUTE  FOR  COMPUTER  SCIENCES  AND  TECHNOLOGY  conducts 
research  and  provides  scientific  and  technical  services  to  aid  Federal  agencies  in  the  selection, 
acquisition,  application,  and  use  of  computer  technology  to  improve  effectiveness  and 
economy  in  Government  operations  in  accordance  with  Public  Law  89-306  (40  U.S.C.  759), 
relevant  Executive  Orders,  and  other  directives;  carries  out  this  mission  by  managing  the 
Federal  Information  Processing  Standards  Program,  developing  Federal  ADP  standards 
guidelines,  and  managing  Federal  participation  in  ADP  voluntary  standardization  activities; 
provides  scientific  and  technological  advisory  services  and  assistance  to  Federal  agencies;  and 
provides  the  technical  foundation  for  computer-related  policies  of  the  Federal  Government. 
The  Institute  consists  of  the  following  centers: 

Programming  Science  and  Technology  —  Computer  Systems  Engineering. 

'Headquarters  and  Laboratories  at  Gaithersburg,  MD,  unless  otherwise  noted; 

mailing  address  Washington,  DC    20234. 

'Some  divisions  within  the  center  are  located  at  Boulder,  CO    80303. 
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PREFACE 

Preparation  of  the  NBS  Self-Study  Manual  on  Optical  Radiation  Measurements  is  an  important 
part  of  our  effort  to  meet  the  needs  of  the  National  Measurement  System  with  respect  to  the 
measurement  of  optical  electromagnetic  radiation.   Significant  needs  exist  for  measurements 
with  uncertainties  of  about  one  percent,  but  the  uncertainties  actually  achieved  are  often   5 
to   10  percent,  or  even  more.   These  relatively  poor  accuracies,  as  compared  to  those  in  many 
other  types  of  physical  measurements,  result  to  a  great  extent  from  the  multidimensional 
character  of  optical  radiation;  radiant  power  is  distributed  and  may  vary  with  position, 
direction,  wavelength,  polarization,  and  time.   Measurement  results  are  also  affected  by 
various  instrumental  and  environmental  parameters.   In  addition,  many  of  those  who  make  meas- 
urements of  optical  radiation  have  little  or  no  training  or  experience  in  this  field  and  are 
limited  in  the  amount  of  time  they  can  devote  to  acquiring  the  needed  information,  under- 
standing, and  proficiency  with  measurement  techniques.   Moreover,  there  are  few  schools  that 
offer  courses  in  radiometry,  and  there  are  almost  no  adequate  texts  or  references  dealing  with 
this  entire  subject. 

The  idea  of  producing  a  Self-Study  Manual  at  NBS  to  try  to  fill  some  of  this  void  was 
developed  by  one  of  us  (HJK)  in  the  latter  part  of  1973.   Detailed  planning  got  under  way  in 
the  summer  of  1974  when  a  full-time  Editor  (FEN)  was  appointed.   The  first  nine  chapters  were 
published  in  a  series  of  four  NBS  Technical  Notes  (910-  series)  up  to  June  1979.   However, 
subsequent  progress  was  slowed  by  administrative  uncertainties  and  problems  that  affected  the 
availability  of  authors,  causing  several  changes  in  scheduling,  and  that  finally  ended  all 
direct  NBS  funding  for  the  Project  on  September  30,  1979.   Welcome  transfers  of  funds  from  a 
number  of  other  federal  agencies  were  obtained  to  carry  the  Project  to  September  30,  1980, 
with  the  hope  that  it  could  then  be  restored  to  the  NBS  budget.   When  that  failed  to  materi- 
alize, outside  funding  was  again  obtained  to  September  30,  1981  (through  FY  81). 

Meanwhile  the  Program  Director  (HJK)  retired  at  the  end  of  May  1981,  and  the  Editor  (FEN) 
on  September  29,  1981.   However,  remaining  FY  81  funds  were  used  for  a  grant  to  Catholic 
University  of  America  (CUA) ,  with  whom  NBS  negotiated  a  cooperative  agreement  under  which  CUA 
hired  the  full-time  Editor  (FEN)  who  continues,  as  before,  at  NBS  (same  address  and  telephone), 
working  closely  with  the  new  Project  Head  (JBS)  there.   CUA  is  also  engaging  the  services  of 
three  other  former  NBS  employees:   Dr.  Kostkowski  as  part-time  consultant  throughout  the 
Project  and  as  author  or  co-author  of  two  more  chapters,  Mr.  Richmond  as  author  or  co-author 
of  two  more  chapters,  and  Dr.  Venable  as  author  of  one  more  chapter,  as  outlined  in  more 
detail  below. 

This  Technical  Note  910-6,  with  Chapter  10  of  Part  I--Concepts,  is  the  first  departure 
from  classical  (geometrical-  or  ray-optics)  radiometry  in  this  Manual.   From  the  beginning, 
due  primarily  to  the  rapid  growth  in  the  use  of  lasers,  there  has  been  a  clear  and  increas- 
ingly urgent  need  for  a  treatment  of  the  radiometry  of  coherent  radiation,  where  geometrical 
(ray)  optics  is  often  an  inadequate  approximation.   However,  when  Dr.  Shumaker  started,  almost 
three  years  ago,  to  look  into  the  possibility  of  writing  a  chapter  on  coherence  in  radiometry, 
it  was  still  not  clear  whether  the  basic  theory  on  which  this  chapter  could  be  based  had  been 
developed  sufficiently  to  make  it  feasible.   Fortunately,  the  recent  rapid  progress  in  this 
area  has  taken  place  and  we  are  delighted  to  be  able,  at  last,  to  respond  to  a  number  of 
repeated  requests  for  the  material  in  this  chapter. 

Our  aim  for  the  entire  Manual  has  been  to  provide  a  comprehensive  tutorial  treatment  that 
is  complete  enough  for  self  instruction.   That  is  what  is  meant  in  the  title  by  "self-study"; 
the  Manual  does  not  contain  explicitly  programmed  learning  steps  as  that  term  often  connotes. 
In  addition,  through  detailed,  yet  concise,  chapter  summaries,  the  Manual  is  designed  to  serve 
also  as  a  convenient  reference  source.   Those  already  familiar  with  a  topic  should  turn 
immediately  to  the  summary  at  the  end  of  the  appropriate  chapter.   They  can  determine  from 
that  summary  what,  if  any,  of  the  body  of  the  chapter  they  want  to  read  for  more  details. 

The  material  in  the  Manual  is  presented  at  the  level  of  a  college  graduate  in  science  or 
engineering,  but  even  for  those  with  facility  in  college  mathematics  and  a  first  course  in 
physics,  it's  not  at  all  easy  reading  in  spite  of  our  best  efforts  at  clarity  and  simplicity. 
This  is  an  unavoidable  result  of  the  primary  aim  ("to  make  one-percent  measurements  common- 
place") coupled  with  the  fact  that  it  must  serve  the  needs  of  so  many  different  fields, 
including  astronomy  and  astrophysics,  mechanical  heat-transfer  engineering,  illumination 
engineering,  photometry,  meteorology,  photo-biology  and  photo-chemistry,  optical"  pyrometry, 
remote  sensing,  military  infrared  applications,  etc.   Chapter  10  is  probably  one  of  the  most 
difficult,  since  it  is  the  first  attempt  by  anyone,  so  far  as  we  know,  to  present  the  radio- 
metric implications  of  coherence  as  they  affect  those  making  such  measurements.   We've  done 
our  best  to  make  it  as  simple  and  readable  as  possible,  but  it  just  isn't  a  simple  subject. 
Most  readers  will  probably  find  that  they  can't  follow  it  all  on  first  reading,  even  taking  it 
slowly  and  carefully.   But  don't  be  discouraged.   Sleep  on  it  and  come  back  to  it,  more  than 
once  if  necessary.   We  wish  we  could  suggest  a  number  of  other  sources  to  consult,  a  technique 
that  often  helps  in  such  situations.   However,  the  whole  viewpoint  and  approach  needed  here  is 
so  different  from  that  of  most  of  the  other  publications  on  coherence  (that  aren't  primarily 
concerned  with  making  radiometric  measurements)  that  it  may  be  difficult  to  reconcile  them 
with  the  material  here.   Furthermore,  enough  care  and  thought  have  gone  into  the  preparation 


of  this  chapter  to  convince  us  that  it  is  sound  and  that,  if  you'll  persevere  and  give  enough 
time  for  each  new  idea  to  become  familiar  to  you,  you'll  find  that  it's  an  adequate  and  useful 

introduction. 

Apparently  it  is  very  difficult  for  those  who  have  not  been  directly  involved  to  realize 
the  full  implications  of  the  situation  just  discussed.   Each  of  us  tends  to  think  of  radiom- 
etry  and  radiometric  measurements  in  terms  of  our  own  immediate  experiences  and  requirements. 
What  each  of  us  would  like  to  have  is  a  set  of  simple,  carefully  designed  procedures  for 
making  our  own  particular  measurements,  with  appropriate  cautions  concerning  likely  sources  of 
error.   However,  the  next  reader  wants  the  same  thing,  but  for  entirely   different   measurements. 
The  desired  radiometric  quantities  to  be  measured  are  different,  the  instrumentation  is  dif- 
ferent, the  ambient  conditions  are  different  —  the  possible  ways  in  which  significant  differ- 
ences may  exist,  in  terms  of  the  radiation  parameters  (position,  direction,  spectrum,  time  or 
frequency  of  modulation  or  fluctuation,  and  polarization)  as  well  as  instrumental  and  environ- 
mental parameters,  are  so  numerous  that  any  attempt  to  cover  them  all  with  a  "cookbook"  treat- 
ment of  specific  measurement  procedures  would  be  impossibly  unwieldy  and  could  never  be  com- 
pleted within  any  conceivable  budget  limitations  short  of  Utopia.   The  only  way  in  which  we 
can  hope  to  effectively  assist  every    reader  who  needs  to  make  one-percent  measurements  is  to 
provide  you,  the  reader,  with  material  which,  with    sufficient    effort    on    your   part,    will  help 
you  to  develop  sufficient  understanding  and  grasp  of  basic  principles  to  solve  your  own  par- 
ticular measurement  problems.   That's  why  we  have  concentrated  on  the  basic  material  of  Part 
I--Concepts  and  why  it  is  not  easy  reading. 

As  an  exception,  we  are  currently  publishing,  in  NBS  TN  910-5,  just  one  chapter  of  Part 
III--Applications  (see  list  of  chapters,  past,  present  and  future,  below).   It  gives  a  detailed 
account  of  some  very  difficult  high-accuracy  field  measurements  of  ultraviolet  solar  terres- 
trial spectral  irradiance,  with  frequent  references  to  applicable  portions  of  the  earlier 
chapters  of  Part  I--Concepts  to  illustrate  how  they  can  be  used  to  achieve  improved  measure- 
ment results.   It  also  illustrates  clearly  the  complexities  of  such  high-accuracy  measurements 
that  make  it  impossible  for  anyone  other  than  you  yourself  to  determine  the  details  of  how 
best  to  make  your  particular  measurements  with  your  particular  instruments  and  constraints  and 
for  your  particular  objectives. 

The  present  chapter  on  coherence  is,  like  Chapter  6  on  polarization,  an  attempt  to  present 
the  subject  as  it  pertains  to  radiometric  measurements.   There  are  many  interesting  and  highly 
important  features  of  coherence  phenomena,  as  well  as  of  polarization  phenomena.   However,  we 
are  not  concerned  with  them  here,  as  such,  but  only  to  the  extent  that  they  affect  the  measure- 
ment of  the  flow  of  power  or  flux  in  a  beam  of  electromagnetic  optical  radiation. 

The  basic  approach  and  focal  point  of  the  treatment  in  this  Manual  is  the  measurement 
equation,  first  introduced  in  detail  but  limited  to  the  radiation  parameters  of  position, 
direction,  and  spectrum,  in  Chapter  5.   We  believe  that  every  measurement  problem  should  be 
addressed  with  such  an  equation,  relating  the  quantity  desired  to  the  data  obtained,  through  a 
detailed  characterization  of  the  instruments  used  and  the  radiation  field  observed,  in  terms 
of  all  of  the  relevant  parameters.   These  parameters  always  include  the  radiation  parameters 
(listed  above) ,  as  well  as  environmental  and  instrumental  parameters  peculiar  to  each  measure- 
ment configuration.   The  objective  of  the  Manual  is  to  develop  the  basic  concepts  required  so 
that  the  reader  will  be  able  to  use  this  measurement-equation  approach.   It  is  our  belief  that 
this  is  the  only  way  that  uncertainties  in  the  measurement  of  optical  radiation  can  generally 
be  limited  to  one,  or  at  most  a  few,  percent. 

The  original,  overall  plan  for  the  Manual  organized  it  into  three  Parts:   Part  I--Concepts, 
Part  II--Instrumentation,  and  Part  III--Applications .   That  was  our  rather  ambitious  plan  when 
we  started  out;  limitations  of  support  and  available  resources,  particularly  available  authors, 
have  determined  how  much,  or  how  little,  we  could  accomplish.   So  far,  as  indicated  above,  we 
have  concentrated  on  Part  I--Concepts.   The  Part  I  chapters  already  published  are: 

1.  Introduction,  by  F.  E.  Nicodemus ,  H.  J.  Kostkowski,  and  A.  T.  Hattenburg 

2.  Distribution  of  Optical  Radiation  with  Respect  to  Position  and 

Direction--Radiance,  by  F.  E.  Nicodemus  and  H.  J.  Kostkowski 

3.  Spectral  Distribution  of  Optical  Radiation,  by  F.  E.  Nicodemus 

and  H.  J.  Kostkowski 

4.  More  on  the  Distribution  of  Optical  Radiation  with  Respect  to 

Position  and  Direction,  by  Fred  E.  Nicodemus 

5.  An  Introduction  to  the  Measurement  Equation,  by  Henry  J.  Kostkowski 

and  Fred  E.  Nicodemus 

6.  Distribution  of  Optical  Radiation  with  Respect  to  Polarization, 

by  John  B.  Shumaker 


7.  The  Relative  Spectral  Responsivity  and  Slit-Scattering  Function 

of  a  Spectroradiometer ,  by  Henry  J.  Kostkowski 

8.  Deconvolution,  by  John  B.  Shumaker 

9.  Physically  Defining  Measurement-Beam  Geometry  by  Using  Opaque 

Barriers,  by  Fred  E.  Nicodemus. 

Chapters  1,  2,  and  3  were  in  TN  910-1;  4  and  5  in  TN  910-2;  6  in  TN  910-3;  and  7,  8,  and  9  in 
TN  910-4.   Now  in  publication,  as  mentioned  above,  is  TN  910-5  with  Part  Ill—Applications : 

Chapter  1.   Measurement  of  Solar  Terrestrial  Spectral  Irradiance  in  the 

Ozone  Cut-Off  Region,  by  Henry  J.  Kostkowski,  Robert  D.  Saunders, 
John  F.  Ward,  Charles  H.  Popenoe ,  and  A.E.S.  Green. 

And  in  the  present  TN  910-6,  we  return  to  Part  I — Concepts: 

10.  Introduction  to  Coherence  in  Radiometry,  by  John  B.  Shumaker. 

Contingent  on  the  availability  of  outside  funding  support,  we  hope  that,  during  the  next 
two  to  three  years,  we  can  continue  the  current  arrangements  with  Catholic  University  of 
America  as  described  above  and  complete  the  following  additional  chapters  of  Part  I — Concepts: 

Distribution  of  Optical  Radiation  with  Respect  to  Time,  by  Fred  E.  Nicodemus 

Linearity  Considerations  and  Calibrations,  by  John  B.  Shumaker 

Spectrophotometry,  by  William  H.  Venable,  Jr. 

Blackbody  Radiation  and  Temperature  Scales,  by  Joseph  C.  Richmond  and/or 
Henry  J.  Kostkowski 

Physical  Photometry,  by  A.  T.  Hattenburg  and/or  Fred  E.  Nicodemus 

Thermal  Radiation  Properties  of  Matter,  by  Joseph  C.  Richmond 

Spectroradiometry  of  Spectral  Lines,  by  Henry  J.  Kostkowski. 

As  stated  above,  the  measurement-equation  approach  is  central  to  the  entire  Manual  and 
the  material  presented  in  each  chapter  needs  to  be  related  to  the  approach  introduced  in 
Chapter  5.   Depending  on  how  well  this  is  accomplished  in  each  chapter,  there  may  be  need  at 
tne  2nd  to  have  a  final  summary  chapter  for  this  purpose,  tying  up  loose  ends  and  putting  the 
whole  Part  I  into  perspective.   Also,  examples  of  various  categories  of  environmental  and 
instrumental  parameters  and  their  significance  could  be  usefully  presented  and  discussed.   But 
it's  too  soon  to  evaluate  this  need  at  present. 

Incidentally,  in  preparing  material  for  the  Manual,  we  have  had  the  pleasure  of  redis- 
covering the  fact  that  the  best  way  to  learn  anything  is  to  try  to  teach  it  to  someone  else. 
The  exercise  of  preparing  tutorial  material  for  such  wide  general  application  has  required  us 
to  analyze  our  own  measurement  activities  in  a  different  way,  broadening  our  understanding  and 
resulting  in  improved  methods  and  more  accurate  results.   Note  that  all  references  to  measure- 
ment accuracy  or  uncertainty  in  this  preface  are  concerned  not  only  with  precision  (relating 
to  the  repeatability  of  measurement  results)  but  also  with  accuracy  (relating  to  agreement 
with  the  "truth"  which,  while  unknowable  in  the  last  analysis,  is  approximated  by  analyses  and 
estimates  based  on  the  widest  possible  experience,  including  agreement  with  measurements  of 
the  same  quantities  by  others,  particularly  when  they  have  used  different  instrumentation  and 
methods  of  measurement) . 

ACKNOWLEDGMENTS .   Several  experts  who  read  an  early  draft  version  of  this  chapter  on  coherence 
offered  many  valuable  suggestions.   We  are  particularly  grateful  to:   William  H.  Carter, 
J.  Christopher  Dainty,  Jack  D.  Gaskill,  Arwind  S.  Marathay,  and  Emil  Wolf. 

As  indicated  above,  we  have  been  entirely  supported  by  funds  from  other  federal  agencies 
since  September  30,  1979.   This  is  one  of  the  first  publications  since  that  date,  so  it  is  one 
of  our  first  opportunities  to  publicly  acknowledge  that  most  welcome  support  from  the  following: 

DoD,  CCG  (tri-service  Calibration  Coordination  Group) ,  Working  Group  for 
IR/Laser  (Mr.  Major  L.  Fecteau  -  Army;  Mr.  Robert  Hinebaugh  -  Air  Force; 
Dr.  Felix  Schweizer  -  Navy) 

Naval  Weapons  Center,  China  Lake,  California  (Dr.  Jon  A.  Wunderlich) 

USAF,  Air  Staff,  Pentagon,  Washington,  D.C.  (Lt.  Col.  John  E.  Dunkle)  through 
USAF,  Space  Division,  Los  Angeles,  California  (Maj .  Robert  Chadbourne) 

Naval  Air  Systems  Command,  Washington,  D.C.  (Mr.  E.  T.  Hooper, 

Mr.  Robert  C.  Thyberg,  Mr.  Frank  Daspit,  and  Mr.  Webb  K.  Whiting) 


USAF,  Avionics  Laboratory,  Wright-Patterson  AFB,  Ohio  (Dr.  Paul  J.  Huffman 
and  Mr.  William  J.  Cannon) 

Defense  Advanced  Research  Projects  Agency  (DARPA) ,  Arlington,  Virginia 
(Dr.  Robert  W.  Fossum  through  Dr.  John  Meson) 

U.S.  Army,  BMD  Advanced  Technology  Center,  Huntsville,  Alabama 
(Dr.  William  0.  Davies) 

NASA  Hq. ,  Environmental  Observation  Program,  Washington,  D.C. 
(Dr.  L.  C.  Greenwood  and  Dr.  S.  G.  Tilford) 

Solar  Energy  Research  Institute  (SERI) ,  Golden,  Colorado 
(Dr.  Keith  Masterson) 

U.S.  Army,  Office  of  Missile  Electronic  Warfare,  Electronic  Warfare  Lab- 
oratory, White  Sands  Missile  Range,  New  Mexico  (Mr.  Thomas  A.  Atherton) 

Office  of  Naval  Research,  Arlington,  Virginia  (Dr.  William  J.  Condell) 

U.S.  Army,  Night-Vision  and  Electro-Optics  Laboratory,  Ft.  Belvoir,  Virginia 
(Mr.  Wayne  Grant) . 

There  have  been  many  others  who,  in  these  days  of  tight  budgets,  were  not  able  to  contribute 
financial  support  to  this  project  when  they  wanted  to  but  who  helped  very  substantially  in 
other  ways.   They  recommended  potential  sources  of  financial  support  and  urged  those  who  were 
able  to  do  so  to  participate.   Just  their  words  of  encouragement  and  of  constructive  criticism 
were  tremendously  helpful.   To  all  of  these  we  are  most  sincerely  grateful. 

We  are  also  deeply  indebted  to  a  great  many  individuals  for  invaluable  "feedback",  con- 
cerning the  technical  content  of  the  Manual,  that  has  encouraged  us  and  has  helped  us  to  put 
these  chapters  together  more  effectively.   We  renew  our  strong  invitation  to  all  readers  to 
submit  comments,  criticisms,  and  suggestions.   In  particular,  we  would  welcome  illustrative 
examples  and  problems  from  as  widely  different  areas  of  application  as  possible. 

We  are  especially  grateful  to  Mrs.  Betty  Castle  for  the  skillful  and  conscientious  effort 
that  produced  the  excellent  typing  and  layout  of  still  another  difficult  text. 

Fred  E.  Nicodemus,  Editor 

John  B.  Shumaker,  Project  Head 

Henry  J.  Kostkowski,  Consultant 

August  1982 


vi 


Contents 

Page 

Part  I.    Concepts 1 

Chapter  10.    Introduction  to  Coherence  in  Radiometry   1 

In  this  CHAPTER 1 

INTRODUCTION  1 

THE  CROSS-SPECTRAL  DENSITY  FUNCTION   4 

CROSS-SPECTRAL  DENSITY  and  RADIOMETRIC  QUANTITIES   8 

EXAMPLES  of  CROSS-SPECTRAL  DENSITY  9 

The  PROPAGATION  of  CROSS-SPECTRAL  DENSITY   13 

MODELS  for  the  CROSS-SPECTRAL  DENSITY 

FUNCTION  —  QUASIHOMOGENEITY  2  5 

The  EFFECTS  of  IMAGING  or  FOCUSING  OPTICS   29 

COHERENCE  EFFECTS  on  a  MONOCHROMATOR 

SLIT- SCATTERING  FUNCTION  37 

MORE  on  IMAGING  —  PUPIL  FUNCTIONS  and  ABERRATIONS  44 

THE  CLASSICAL  RADIOMETRIC  QUANTITIES  4  7 

POLARIZATION  5  3 

SUMMARY  of  CHAPTER  10 55 

References 59 


List  of  Figures 

10.1  Alternating  regions  of  constructive  (C-clear)  and  destructive 

(D-dotted)  interference  in  the  radiation  field  about  two  point 
(small)  sources.  These  regions  will  wriggle  and  writhe  as  the 
frequencies  (or  amplitudes  or  phases)  of  the  oscillators  change   2 

10.2  A  train  of  plane  waves  incident  on  the  X-Y  plane  at  an  angle   6 
from  the  Z-axis  in  the  X-Z  plane  and  with  a  Gaussian  amplitude 

distribution  12 

10.3  The  propagation  of  cross-spectral  density  W, (P, ,P2) ,   assumed 
known  for  all  pairs  of  points   P,   and  P_   in  one  cross  section 

of  a  beam  (source  plane) .   Optical  path  lengths   s,   and   s„   extend 

to  corresponding  points  P'   and  PI   in  a  second  plane  where  the 

cross-spectral  density  is  wanted  (detector  plane) 13 

10.4  A  typical  field  point  P'(x',y',z')   is  shown  at  a  distance  R   from 
the  origin  and  a  distance   s   from  a  point   P(x,y,0)   in  the  source 

plane   (X-Y  plane) 14 

10.5  The  sine  function:   sine  x  e  sin  (ttx)  /  (ttx)   18 


Vll 


Page 

10.6  (a)   Cross-spectral  density  propagation  from  source  plane  to 

intermediate  optical  element  (mirror) 

(b)  Cross-spectral  density  propagation  from  intermediate 
optical  element  (mirror)  to  detector  plane 

(c)  Cross-spectral  density  propagation  from  source  plane 

through  intermediate  optical  element  (lens)  to  detector  plane  30 

10.7  Relative  distribution  in  the  x-direction  of  fluence  rate  in  the 
image  of  a   50  [ym]-wide,  uniform,  incoherent  source  for  various 
lens  aperture  widths.   Plots  of  the  part  of  eq.  (10.116)  within 
curly  brackets  for   X  =  514  [nm] ,   a  =  25  [urn] ,   R  =  R"  =  0.8  [m] , 
and   a'  =  0.01  [m]  ; a'  =  0.02  [m]  ; a'  =  0.04  [m]  ; 

a'  =  0.08  [m] 36 

10.8  Relative  distribution  in  the  x-direction  of  fluence  rate  in  the 
image  of  a   50  [ym]-wide,  uniform,  coherent  source  for  various 
angles  of  incidence  of  the  beam.   Plots  of  the  double  integral 
of  eq.  (10.119)  for  the  same  parameter  values  as  in  figure  10.7 
except   a'  =  0.0375  [m]   and   A   =  1.    6  =  0.00  [rad] ; 

6   =  0.03  [rad]; 6   =  0.04  [rad]; 6   =  0.05  [rad] 38 

10.9  Relative  distribution  in  the  x-direction  of  fluence  rate 

in  the  image  of  a   50  [urn] -wide,  uniform,  partially  coherent 

source  at  normal  incidence  for  various  degrees  of  partial 

coherence.   Plots  of  the  double  integral  of  eq.  (10.124) 

for  the  same  parameter  values  as  in  figure  10.8  and  area 

normalized,  in  the  interval  shown,  to  the  area  under  the 

9   =0   curve  in  figure  10.8.   a  /R   =0.00  (coherent); 

o  3  o   o  - 

a  /R   =0.04; a  /R   =0.045;   a  /R   =0.05; 

o   o  o   o  o   o 

a  /R   =  °°  (incoherent 41 

10.10  Two  samples  of  the  relative  distribution  in  the  x-direction  of 
fluence  rate  in  the  image  of  a   50  [ym]-wide  coherent  source  con- 
sisting of   1  [yra]   cells  with  amplitude  and  phase  constant  within 

each  cell  but  random  from  cell  to  cell 43 
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Part  I.    Concepts 

This  is  the  sixth  in  a  series  of  Technical  Notes  (910-  )  entitled 
"Self-Study  Manual  on  Optical  Radiation  Measurements".   It  contains 
Chapter  10  of  Part  I  of  this  Manual.   Additional  chapters  will  continue 
to  be  published,  similarly,  as  they  are  completed.   The  Manual  is  a  com- 
prehensive tutorial  treatment  of  the  measurement  of  optical  radiation 
that  is  complete  enough  for  self  instruction.   Detailed  chapter  summaries 
make  it  also  a  convenient  authoritative  reference  source. 

In  this  chapter  we  introduce  the  reader  to  the  basic  field  quantity 
of  the  theory  of  partial  coherence,  the  cross-spectral  density  function, 
and  show  how  it  is  used  to  describe  radiation  fields.   We  discuss  the 
propagation  of  cross-spectral  density  along  a  beam  and  illustrate  this 
propagation  with  calculations  of  diffraction  effects  in  a  number  of  simple 
aperturing  and  imaging  examples.   For  instance,  we  treat  in  considerable 
detail  one  of  the  most  common  radiometric  situations  in  which  coherence 
effects  can  manifest  themselves,  the  measurement  of  the  slit-scattering 
function  of  a  monochromator .   Among  other  things,  this  treatment  shows 
that  laser  illumination  of  a  monochromator  entrance  slit  must  be  nearly 
on-axis  to  avoid  serious  slit-scattering  function  distortions.   Finally, 
we  present  the  relationships  between  the  cross-spectral  density  function 
and  the  classical  radiometric  quantities,  such  as  spectral  radiance. 

Key    Words:      coherence;  cross-spectral  density;  incoherence;  interference; 
partial  coherence;  radiometry;  spectroradiometry;  wave  optics. 

Chapter  10.    Introduction  to  Coherence  in  Radiometry 

by  John  B.  Shumaker 

In  this  CHAPTER.    We  attempt  to  provide  an  introduction  to  the  theory  of  partial  coherence 
and  to  its  significance  for  radiometry  through  some  illustrations  of  its  application  to 
radiometric  situations.   We  begin  with  the  cross-spectral  density  function  which  is  used  to 
describe  the  radiometric  field,  including  phase  relationships,  much  as  the  spectral  radiance 
function  describes  the  field  in  classical  radiometry.   A  few  examples  of  this  function  are 
presented  and  discussed.   Then  we  give  the  equation  for  the  propagation  of  cross-spectral 
density  from  one  surface  (such  as  a  source)  to  another  (such  as  a  detector) .   This  equation 
is  then  applied  to  such  radiometric  situations  as  the  irradiation  of  a  surface,  image  forma- 
tion, and  the  calculation  of  the  slit-scattering  function  of  a  monochromator.   Although,  in 
principle,  the  cross-spectral  density  function  is  measurable  in  detail,  it  is  usually  more 
practical  to  assume  some  simple  functional  form  or  model  for  this  function;  some  common 
models  are  given.   Finally,  we  discuss  the  definition  of  the  classical  radiometric  quanti- 
ties, such  as  radiance,  in  terms  of  the  cross-spectral  density  function  and  we  mention 
briefly  the  extension  of  the  elementary  theory  presented  in  this  chapter  to  a  vector  theory 
that  includes  polarization. 

INTRODUCTION.  In  treating  the  fundamentals  of  radiometry  (optical  radiation  measurements) 
in  the  first  nine  chapters  of  Part  I  of  this  Manual,  we  have  so  far  limited  ourselves  to  the 
domain  of  geometrical  or  ray  optics,  where  it  is  assumed  that  the  dimensions  of  all  optical 

components  are  large  compared  to  a  wavelength   (^10  '  [m] )   and  all  times  are  long  compared 
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to  a  period   (^10    )  [s] ) ,   etc.    Diffraction  and  interference  phenomena,  which  are  best 

understood  in  terms  of  wave  optics,  rather  than  ray  optics,  have  so  far  been  ignored.   We 


xUnit  symbols  are  all  enclosed  in  square  brackets  to  emphasize  dimensionality  and  encourage 
routine  checks  of  unit-dimension  consistency,  as  discussed  in  reference  [10.1]. 


Figure  10.1.    Alternating  regions  of  constructive  (C-clear)  and  destructive 
(D-dotted)  interference  in  the  radiation  field  about  two  point 
(small)  sources.   These  regions  will  wriggle  and  writhe  as  the 
frequencies  (or  amplitudes  or  phases)  of  the  oscillators  change. 


are,  however,  very  conscious  of  the  fact  that  there  is  a  rapidly  increasing  interest  in 
measurements  of  laser  radiation  which,  being  highly  coherent,  require  a  wave-optics  treat- 
ment.  Also  a  rapidly  increasing  amount  of  material  has  appeared  in  the  literature  on  the 
theory  of  optics  and  its  significance  for  radiometry.   Most  of  this  material,  however,  has 
been  at  a  theoretical  level  where  it  is  not  too  useful  for  those  actually  engaged  in  making 
radiometric  measurements.   It  is  the  objective  of  this  chapter  to  begin  to  bridge  this  gap 
and  to  present  the  results  of  some  of  these  theoretical  treatments  in  a  form  that  will  make 
clear  their  application  to  radiometric  measurements. 

The  theory  of  partial  coherence  is  filled  with  complex  details  and  some  readers  may 
find  that  they  "loose  the  thread"  at  times.   If  so,  the  summary  at  the  end  of  the  chapter 
may  prove  helpful  to  bring  more  clearly  into  focus  the  pattern  of  the  overall  treatment.   In 
addition,  some  recent  publications  [10.2] 1  to  [10.8]  provide  further  reading  and  references 
to  the  growing  literature  on  partial  coherence.   Our  treatment  will  be  found  to  be  different 
in  some  respects  from  these  largely  theoretical  references  due  to  our  preoccupation  with 
radiometric  measurements  and  to  our  recognition  throughout  this  Manual  that  practical 
definitions  of  measurable  radiometric  quantities  (intensity,  radiance,  etc.,  and  cross- 
spectral  density,  introduced  below)  must  include  a  dependence  upon  the  many-dimensional 
resolution  interval  of  the  measuring  radiometer  (e.g.  spectral  pass-band).   This  latter 
difference  is  minor,  being  limited  to  the  purely  formal  connection  between  the  electro- 
magnetic field  and  radiometric  quantities.   Nevertheless  it  has  proven  difficult  to  cite 
specific  references  for  the  ideas  and  equations  presented  here  and  we  have  done  so  only  in 
rare  instances.   However,  we  wish  to  emphasize  that  almost  nothing  in  this  chapter  is 
original;  it  can  all  be  found  somewhere  in  the  literature  in  some  form. 

We  will  begin  by  trying  to  describe  the  nature  of  a  classical  radiometric  field  at  the 
microscopic  level.   When  the  fields  from  separate  oscillators  are  superimposed,  the  electric 
and  magnetic)  field  vectors  add,  producing  a  resultant  field  vector,  and  it  is  the  effect 
that  this  resultant  field  has  on  a  detector  which  we  ultimately  interpret  as  a  measure  of 
incident  radiant  power.   If  two  point  (very  small)  sources  have  roughly  similar  frequencies 
and  are  oriented  so  that  their  fields  have  the  same  polarization,  then  there  will  be  regions 
of  space,  whose  minimum  dimension  is  at  least  half  a  wavelength,  where  the  electric  field 

vectors  will  lie  in  the  same  direction  at  the  same  time  and  will  tend  to  reinforce  each 
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other  over  a  time  period  of  at  least  a  cycle  or  two  of  the  optical  frequency  (^10     [s]). 

Such  regions  will  obviously  alternate  with  regions  where  the  field  vectors  lie  in  opposite 

directions  and  tend  to  cancel  during  this  same  time  period  (see  figure  10.1).   In  general, 

if  the  two  oscillators  are  changing  frequency  or  amplitude  or  phase  these  various  regions 

will  move  about  and  change  shape  and  size  with  time.   When  these  regions  of  reinforcement 

and  cancellation  are  large  enough  and  move  slowly  enough  to  be  detected  by  the  eye  or  by  an 

optical  instrument,  we  call  the  phenomenon  interference  and  we  say  that  the  light  from  the 

two  sources  possesses  coherence.   If  no  such  regions  of  alternating  radiant  power  are 

discerned  by  the  best  instruments  then  we  say  that  the  field  is  incoherent  and  that  there  is 

no  interference.   Presumably  the  interference  (coherence)  is  still  there  but  the  temporal  or 

spatial  response  of  the  instrumentation  is  inadequate  to  resolve  it.   Notice  that  the 

regions  of  positive  and  negative  reinforcement  in  figure  10.1  generally  increase  in  width 

with  increasing  distance  from  the  source.   This  suggests  that,  although  the  field  at  one 

point  may  be  judged  incoherent,  it  may  show  distinct  coherence  effects  at  points  more 

distant  from  the  source.   As  we  shall  see,  this  increase  in  the  degree  of  coherence  with 

distance  from  the  source  is  quite  real;  the  light  from  the  most  incoherent-seeming  source, 


:Figures  in  brackets  indicate  literature  references  listed  at  the  end  of  this  Technical  Note. 
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such  as  a  star,  is  highly  coherent  at  astronomical  distances,  as  is  evidenced  by  stellar 
interferometer  measurements  [10.2]. 

Diffraction  is  the  name  associated  with  departures  from  regular  (geometrical-  or  ray- 
optics)  propagation  of  light.   The  spreading  of  light  after  passage  through  a  pinhole  or 
slit  is  an  example  of  diffraction.   Beyond  the  slit  or  pinhole  the  light  beams  from  the 
different  parts  of  the  opening  interfere  producing  the  diffraction  effect.   Thus,  inter- 
ference and  diffraction  are  essentially  the  same  and  originate  in  the  wave  nature  of  light. 
All  of  these  phenomena  are  included  in  the  theory  of  partial  coherence  which  is  the  subject 
of  this  chapter. 

Fortunately  in  most  optical  radiation  measurement  situations  diffraction  or  inter- 
ference effects  are  quite  negligible  so  that  it  is  adequate  to  use  the  'simple1  geometrical 
optics  and  classical  radiometry  which  has  been  the  subject  of  previous  chapters.   When  these 
effects  cannot  be  neglected,  their  quantitative  treatment  by  the  theory  of  partial  coherence 
can  be  quite  concisely  stated  [eq.  (10.36),  p.  13,  for  example].   However,  the  physical 
detail  and  mathematical  complications  which  are  required  to  handle  any  but  the  most  trivial 
optical  situations  force  the  adoption  of  a  formidable  chain  of  physical  idealizations  and 
mathematical  and  numerical  approximations.   This  is  not  to  say  that  an  accurate  calculation 
could  not  be  carried  out  for  a  real  problem  of  sufficient  importance  but  that  an  approximate 
calculation  coupled  with  suitable  experimental  measurements  will  usually  be  more  efficient 
and  convincing.   Thus,  we  will  approach  the  treatment  of  coherence  as  providing  a  guide  to 
measurement  problem  areas  and  to  estimates  of  errors  and  corrections  to  geometrical  optics 
when  coherence  effects  are  small.   In  this  spirit  we  will  frequently  make  assumptions  and 
approximations  in  this  chapter  which  would  otherwise  be  intolerable. 

THE  CROSS-SPECTRAL  DENSITY  FUNCTION.    In  this  chapter  we  will  treat  what  is  known  as  the 
scalar  theory  of  coherence;   we  deal  with  only  a  single  polarization  component  of  the 
electromagnetic  field  —  say  the  x-component.   In  a  later  section  we  will  discuss  briefly 
the  complications  introduced  by  the  vector  nature  of  the  electromagnetic  field.   Now  we  want 
to  attempt  to  describe  the  connection  between  the  electromagnetic  field  and  the  quantity 
known  as  the  cross-spectral  density  (or  the  mutual  spectral  density) .   The  cross-spectral 
density  can  be  regarded  as  a  radiometric  quantity  somewhat  more  basic  than  spectral  radiance 
in  that  it  provides  an  even  more  detailed  specification  of  the  radiometric  field.   Our 
discussion  of  the  connection  between  the  cross-spectral  density  function  and  the  electro- 
magnetic field  is  not  intended  to  be  rigorous  but  merely  to  provide,  perhaps,  some  quick 
insight  into  the  nature  of  this  function.   (For  a  rigorous  derivation  of  this  relationship 
we  recommend  especially  reference  [10.3].)   It  would  be  entirely  sufficient  to  take  the 
phenomenological  view  that  it  is  convenient  to  define  such  a  radiometric  field  quantity  in 
terms  of  which,  as  we  shall  see,  not  only  can  all  of  the  classical  radiometric  quantities  be 
expressed  but  also  most  interference  and  diffraction  effects  as  well.   Everything  of  signif- 
icance to  radiometry  can  be  viewed  in  this  way  without  any  concern  for  a  more  basic  picture 
of  the  meaning  of  the  cross-spectral  density  function. 

Let  us  denote  the  x-component  of  the  electromagnetic  field  as  a  function  of  time  at  a 
point  P  by  E(P,t).   We  imagine  that  E(P,t)   will  be  roughly  sinusoidal  (as  a  function  of 
time)  but  with  a  continuously  changing  amplitude  and  frequency.   The  optical  frequencies  we 
are  concerned  with  are  in  the  neighborhood  of   10    to   10    [Hz] .   Although  someday 
sensors  may  be  developed  which  are  fast  enough  to  measure  E(P,t),   at  present  all  radio- 
metric detectors  are  square-law  detectors  whose  output  is  at  least  approximately  proportional 

-9 

to  the  integral  of  the  square  of  E(P,t)   over  a  response  time  of   10    [s]   or  more.   The 

square  of  a  sinusoidal  field  is  proportional  to  the  magnitude  of  the  Poynting  vector  of 
electromagnetic  theory,  and  the  Poynting  vector  is  associated  with  the  energy  per  second 


which  crosses  a  unit  area  perpendicular  to  the  vector  direction.   Thus,  a  proportionality 
relationship  is  generally  assumed  between  the  incident  radiometric  power  and  the  response  of 
a  detector  or  the  square  of  the  electric  field:1 

W   -   AtTAa"   f   |  E2(P,t).dP.dt  (10.1) 

At  AS 

where   F  (P  )   is  the  measurable  average  fluence  rate  (see  Chapter  4,  p.  17  [10.9])  and   k 

is  some  proportionality  constant.   The  time  required  for  a  measurement  is   At   and   Aa   is 

2  .  2 

the  cross-sectional  area   (=irr  )   of  a  small  opaque  spherical  detector  surface   AS  (  =  4Trr  ) 

centered  at  the  point  of  interest   P  .    The  integrations  are  over  the  detector  surface  and 

measurement  time.   If  the  electromagnetic  field  passes  through  a  narrow  pass-band  spectral 

filter  centered  at  wavelength   X    before  being  detected  then  the  radiometric  quantity 

incident  on  the  detector  is  the  average  spectral  fluence  rate  over  the  spectral  pass  band   AX 

of  the  filter: 

i       ,    E2 (P,t;X  ) 

Ft,A<VX0>   =   At^Aa-J   j     AX    °   -dP'dt'  <10-4) 

At  AS 

where   E(P,t;X  )   is  the  filter-modified  field  at  the  point   P.   As  we  have  mentioned  in 
previous  chapters  a  radiometric  quantity,  such  as  spectral  fluence  rate,  at  a  particular 
point,  wavelength,  and  instant  should  be  identified  with  an  extrapolation  from  a  sequence  of 
measurements  of  successively  narrower  spatial,  spectral  and  temporal  resolution.   But, 
because  none  of  these  instrument  resolution  intervals  can  ever  be  taken  as  zero  due,  ulti- 
mately, to  the  noise-imposed  limit  of  the  minimum  energy  required  for  meaningful  measurement, 
there  is  in  practice  a  set  of  smallest  values  for   At,  Aa,  and  AX   consistent  with  the  best 


*We  will  frequently  specify  the  function  arguments   x,y,z   simply  as   P.    That  is, 
E(P)  =  E(x,y,z).    An  integral  such  as  in  eq.  (10.1)  containing   dP   means  that  the  point   P 
explores  the  region  defined  by  the  integration  limit.   In  eq.  (10.1),  since   AS   is  a  spheri- 
cal surface,  the  integration  over   P   is  understood  to  be  2-dimensional  and  could  be  written 
out  in  full  as 

2  C2tt    f71      2  2 

E    (P)-dP      =  E    (r,e,cf>)  -r    •  sinO -d9 -d<J> .  (10.2) 


AS  0         0 

Usually  in  this  chapter   P   will  be  confined  to  an  area  in  the   X-Y   plane,  and  an  integral 
would  again  be  understood  to  be  2-dimensional.   Thus 


.  .  .  dP  =  ...  dx-dy  (10.3) 

AA  AA 

or  an  equivalent  in  polar  coordinates.   This  shorthand  use  of   P   and  dP   will  generally  be 
confined  to  generalizations  and  concepts;  equations  intended  for  numerical  evaluation  will 
be  written  out  fully.   Also  note  that  we  use  the  parameter  subscript  to  denote  a  derivative: 
Ffc  =  dF/dt;   F   x  =    d2F/(dt-dX) ;  etc.,  as  in  [10.10]. 


foreseeable  technology  which  can  be  used  to  provide  an  adequate,  unique  definition  of  radio- 
metric quantities.   For  the  purpose  of  establishing  a  correspondence  between  electromagnetic 

theory  and  radiometry  then  we  perform  a  suitable  average  as  in  eq.  (10.4)  of  electromagnetic 
fields  over  resolution  intervals   At,   Aa,   and   AX   with  these  values  and  must  assume  that 
radiometric  quantities  have  been  measured  at,  or  extrapolated  to,  this  same  resolution. 
Normally,  in  radiometric  practice  (intercomparing  beams  by  measurements  with  the  same  instru- 
mentation) these  resolution  intervals  cancel  out  so  that  their  values  are  largely  irrelevant, 

It  is  customary  in  the  treatment  of  oscillating  quantities  in  physics  and  engineering 
to  introduce  complex  variables. 1      Since  a  complex  variable  requires  two  numbers  for  its 
specification  this  notation  provides  a  compact  way  to  keep  track  of  both  amplitude  and  phase 
information  with  a  single  variable.   Whatever  the  form  of   E(t)   it  can  be  expressed  by  a 
Fourier  integral2 


E(t)   =    A(v)  •cos[y(v)  -  2irvt]«dv,  (10.10) 

0 


xWe  will  use  the  following  properties  of  complex  numbers:   A  complex  number,   V,   can  always 
be  expressed  as  the  sum  of  a  real  part,   a,   and  the  product  of   i  by  an  imaginary  part,   b, 
where   a  and  b  are  ordinary  real  numbers: 


(10.5) 


V  e  a  +  i*b 
Re(V)  s   a,     Im(V)  =   b 

2 
and   i,   of  course,  satisfies   i   =  -1.    Equivalently  V  can  be  expressed  as 

V  =  R»exp(i«6)  =  R* (cos6  +  i«sin6)  (10.6) 
where  R  and   6   are  also  ordinary  real  numbers  related  to   a  and  b  by 

R  =  \a2  +  b2,     tan6  -  b/a.  (10.7) 

R   is  called  the  modulus  of  V,   and   9   is  called  the  argument  or  phase  of   V.    The  com- 
plex conjugate  of  V  denoted  by  V*   is  obtained  by  replacing   i   by   -i   throughout  any 
expression  for  V  which  involves  only  real  quantities  and   i 

V*  =  a  -  i«b 

(10.8) 
=  R-exp(-i«  0) . 

We  also  define 

|V|  =  R  =  ^V«V*.  (10.9) 


2The  signs  of  the  phase  terms  y(v)  and  2irvt  in  eqs.  (10.10)  and  (10.11)  are  entirely 
arbitrary;  the  evolution  of  E(t)  or  V(t)  with  time  is  the  same  for  either  sign.  The 
present  choice  leads  to  agreement  with  most  of  the  coherence  literature. 


where  A(v)   and  y (v)      are  uniquely  defined  by  the  function   E(t),  and   v   is  the  frequency. 
The  complex  field  function   V(t)   which  represents   E(t)   is  then  obtained  by  replacing  the 
cosine  function  by  a  complex  exponential  function 


V(t)   =    A(v)  «exp[iY(v)  -  2Trivt]«dv.  (10.11) 

0 

Thus  the  real  part  of   V(t)   is  identically   E(t).   The  imaginary  part  can  be  regarded  as  a 

mathematical  device  for  coding  the  phase  information.   The  substitution  of  V(t)   for   E(t) 

in  equations  which  are  linear  in   E (t)   presents  no  problem;  one  can  always  take  the  real 

part  of  such  equations  to  recover  valid  physical  relationships.   However,  equations  quadratic 

2 
in   E(t)   require  the  substitution  of   JsV*(t)»V(t)   for   E  (t)   and  are  valid  only  in  the 

limit  of  long  time  averages.  Fortunately  in  radiometry  this  limit  is  always  reached;  the 

measurement  time   At   is  long  compared  to  the  natural  time  scales  of  the  field   1/v   and 

1/Av  ^  —  *T-p    In  accord  with  this  customary  notation  then  we  will  write  eq.  (10.4)  in 

complex  notation  as 

K      (         (    V*(P,t;A  ).V(P,t;X  ) 

F.  ,  (P  ,A  )   =     *  .        -^ —  -dP-dt,  (10.12) 

t,A   o   o       2At«Aa  J   J         AX 

At  AS 

with  the  real  part  of  V(P,t;A  )   being  simply  the  electric  field   E(P,t;A  )   at  the  point   P 
and  time   t,   as  transmitted  by  a  very  good  spectral  filter  of  passband   AA   centered  near 
A  .    The  complete  function  V,   including  its  imaginary  part,  is  constructed  from  E   by 
decomposing  E   into  its  Fourier  components  and  reassembling  as  indicated  by  eqs .  (10.10)  and 
(10.11).   V*   is  the  complex  conjugate  of  V.    Actually  we  never  perform  a  detailed  mathe- 
matical averaging  over  time  and  area  because  we  never  know  V(P,t;A)   for  optical  frequencies. 
The  significance  of  eq.  (10.12)  is  primarily  to  show  that  a  definite  connection  exists  between 
the  electromagnetic  field  and  a  measurable  radiometric  quantity.   For  this  reason  it  is 
sufficient  to  write  eq.  (10.12)  as 

Ft,A(W   "   AT  *  <V*(P,t;Ao).V(P,t;Ao)>  ,  (10.13) 

where  the  angular  brackets  represent  the  averaging  over  measurement  time  and  area.   We 
have  also  taken  the  liberty  of  incorporating  the  proportionality  constant   k/2   into  the 
definition  of   V   since  we  will  never  have  occasion  to  assign  numerical  values  or  unit 
dimensions  to  V. 

In  eq.  (10.13)  we  have  shown  explicitly  the  dependence  of  the  fluence  rate  and  the 
complex  field  V(P)   upon  the  field  point  of  interest  in  order  that  we  can  now  introduce  a 
generalization  of  eq.  (10.13)  by  allowing  the  quantities   V*   and   V   to  refer  to  two  dif- 
ferent space  points: 

WX(P1,P2,A)   =  j^    •  <V*(P1  +  AP,t;A.)  -y(P2  +  AP,t;A)>  .  (10.14) 

The  angular  bracket  average  in  eq.  (10.14)  consists  of  an  average  over  a  small  time  interval 
and  over  identical  small  areas  about   P,   and   P?   as  the  small  displacement   AP   covers  the 


i,(P.,P_),   which  may  be  complex,  is  called  the  cross-spectral  density  function.   Usually  we 


measurement  resolution  area.1   We  assume  that  the  radiometric  field  is  steady  in  the  sense 
that,  although  V   depends  upon  time  strongly,  the  average  over  measurement  time  implied  by 
the  angular  brackets  in  eq.  (10.14)  changes  slowly,  if  at  all,  from  one  measurement  time  to 
the  next  so  that  we  can  ignore  any  time  dependence  for   W,  (P,,P_,A).    We  will  also  hereafter 
usually  not  display  the  dependence  upon  wavelength;  we  will  be  dealing  exclusively  with 
spectral  radiometric  quantities  and  will  leave  it  understood  that  they  depend  upon   A. 

apply  this  function  to  pairs  of  points  lying  in  a  particular  surface  such  as  a  source  plane, 
an  image  plane,  or  an  entrance  aperture.   It  is  immaterial  whether  the  points   P,   and   P_ 
are  points  in  a  real  physical  object  or  are  simply  two  points  in  a  field. 

CROSS-SPECTRAL  DENSITY  and  RADIOMETRIC  QUANTITIES.    When  the  two  points  coincide  the  cross- 
spectral  density  function  reduces  to  the  spectral  fluence  rate 

WA(P1,P1)   =   Ft  X(P±) .  (10.15) 

This  quantity,  the  spectral  fluence  rate,  or  its  integral,  fluence  rate,  is  usually  called 
the  optical  intensity,  or  simply  the  intensity,  in  the  coherence  literature.   The  relation- 
ship between  spectral  fluence  rate  from  a  distant  small  source  and  the  quantity  known  in 
radiometry  as  the  spectral  radiant  intensity  of  that  source  is 

Pt  X(P)   =   IA(9,(f>)/R2  (10.16) 

where   R   and   9,(J>   are  the  distance  and  direction  from  the  source  whose  spectral  radiant 
intensity  is   I,   in  the  direction  of  the  point  where  the  description  of  the  field  is 
desired.   The  spectral  radiant  intensity   I,   is  a  property  of  a  (possibly  hypothetical) 
source  while  the  spectral  fluence  rate  F.  , (P)   is  a  measure  of  the  radiant  power  con- 
centration  in  the  field  at  the  point   P. 

Finally,  for  completeness,  we  present  a  relationship  between  the  cross-spectral 
density  in  a  plane  and  the  spectral  radiance  in  that  plane:2 


!By  analogy  to  the  definition  of  fluence  rate  the  area  average  is,  strictly  speaking,  a 
single  surface  integral  in  which  AP  moves  over  a  small  sphere  so  that  P,  +  AP   and 
P_  +  AP   simultaneously  describe  surfaces  centered  on  P-,   and  P„.    The  measurement  resolu- 
tion area  is  the  area  of  a  great  circle  of  this  small  sphere. 


2 Where  it  may  improve  clarity,  we  will  sometimes  separate  the  coordinates  of  P,   and  P_ 
by  a  semicolon.   For  example 

W(P1,P2)  =   W(x1,y1,z1;x2,y2,z2) .  .   (10.17) 

Also  we  will  frequently  omit  the   z   coordinate  in  such  expressions  when  the  points  are 
confined  to  a  plane   z  =  const.    In  this  case  we  distinguish  different  planes  by  coordinates 
labeled  by  primes  ( ' )  and  use  notations  like 

W(x1,y1;x2,y2)  =  W(x1,y1, z;x2 ,y2 , z)  (10.18) 

and       W*(xJfy»jx»fy£)  =  W(x',yJ,z' ;x^,y^,z' ) .  (10.19) 


LA(x,y,e,<)>)       =      S2ji    .11    wA(x+C/2,    y+n/2;    x-^/2,    y-n/2) 

•    expl-   — j— -   •  (g»sin6«cos<(i   +   rv  sine  •  sin<J>)     'd^-dn 


X' 

(10.20) 


Here,   9   and   <J>   are  the  usual  polar  coordinates  of  the  viewing  direction.   The  integrations 
are  over  the  entire  cross-section  of  the  beam  or  over  the  source  if  the  surface  is  a  source 
plane.   We  will  return  to  a  discussion  of  this  equation  later  in  this  chapter.1   For  the 
present  we  will  confine  our  attention  to  the  spectral  fluence  rate  because  of  its  simple, 
direct  connection  to  the  cross-spectral  density  function.   As  is  evident  from  the  discussion 
of  fluence  rate  in  Chapter  4  [10.9],  this  radiometric  quantity  is  closely  related  to  spectral 
irradiance.   At  most  points  in  most  fields  of  interest  in  radiometry  the  fluence  rate  will 
have  a  larger  value  than  the  irradiance  by  a  factor,  of  the  order  of  unity,  which  changes 
very  slowly,  if  at  all,  with  position  in  the  field  and  with  wavelength.   So,  for  the  kind  of 
accuracy  with  which  we  are  concerned  in  this  chapter,  the  two  quantities  are  virtually  inter- 
changeable. 

EXAMPLES  of  CROSS-SPECTRAL  DENSITY.    In  coherence  theory  the  cross-spectral  density  function 
can  be  regarded  as  a  basic  field  entity  much  as  spectral  radiance  is  the  basic  field  entity 
of  conventional  geometric-optics  radiometry.   Unfortunately  the  physical  meaning  of  the 
cross-spectral  density  function  is  more  abstract  than  that  of  radiance.   The  cross-spectral 
density  embodies  simultaneously  information  about  the  radiometric  field  strength  and  infor- 
mation about  the  phase  correlations  at  all  pairs  of  points.   Needless  to  say,  it  is  not  a 
function  which  one  can  tabulate  or  graph  easily,  depending,  as  it  does  upon  at  least 
x, ,  y.  ,  x_ ,  y_ ,   and   A   and  having  both  real  and  imaginary  components. 

To  gain  some  insight  into  the  nature  of  the  cross-spectral  density  function  let  us 
examine  the  behavior  of   V(P,t)   for  a  typical  optical  field  and  calculate   W,  (P-,P_)   from 
the  definition,  eq.  (10.14).   First  let  us  assume  that  the  field  is  generated  by  a  single 
distant  oscillator.   Since  we  are  particularly  interested  in  the  phase  or  argument  of  the 
complex  field  V   let  us  assume  that  the  oscillator  is  a  dipole  rotating  in  a  plane  con- 
taining the  point   P   of  interest.   Then  the  phase  of  the  complex  function   V(P,t)   is  just 
the  angle  of  rotation  of  the  dipole.   If  we  think  of   V   as  a  vector  from  the  origin  to  the 
point  V   in  the  complex  plane  (where  the  X-axis  is  the  real  axis  and  the  Y-axis  is  the 
imaginary  axis)  then  the  vector  V   rotates  in  synchronism  with  the  physical  rotating  dipole. 
Classically,  as  the  dipole  rotates  and  radiates  it  slowly  loses  energy  so  the  amplitude  of 
the  field  from  a  radiating  dipole  will  slowly  decay  with  time.   Therefore  we  can  imagine  the 
vector   V   growing  from  zero  rather  quickly  and  then  gradually  spiralling  in  with  constant 
angular  velocity  until  after  a  long  time  it  vanishes  again  at  the  origin.   A  real  field  will 
be  a  superposition  of  fields  from  a  great  many  such  oscillators  with  different  frequencies, 
starting,  and  sometimes  stopping  due  to  collisions,  at  random  times.   The  resultant  complex 
field  vector  V   will  then  basically  perform  a  rotation  at  some  mean  frequency  but  will 
undergo  frequent  random  angular  accelerations  and  decelerations  and  radial  expansions  and 
contractions.   If  we  write   V(t)   as 

V(t)   =   A(t)  -exptiY  (t)  ]  «exp[-2Triv  t]  (10.21) 

where   v    is  this  mean  frequency,  then  the  complex  quantity  A (t) -exp [iy (t) ]   viewed  as  a 
vector  in  the  complex  plane  will  exhibit  the  motion  of   V(t)   exclusive  of  the  rotation  at 


'See  eq.  (10.149)  on  page  49. 


frequency   v  :   it  will  trace  a  kind  of  complicated  random  tangle  which  contains  only  the 
random  angular  and  radial  motions  of   V(t). 

In  the  product   V*  (P,  ,  t)  *V(P2 ,  t)   the  steady  rotational  motion,   exp(2iTiv  t)  ,   of  the 
complex  vectors  cancels  out  so  that  the  cross-spectral  density  function   W,  (P,  ,P_)   is  just 
the  average  of  the  product  of  two  randomly  dancing  vectors  like   A (P, t) -expfiy (P, t) ] .   If   P 
and   P2   are  very  close  together  then  both  points  will  view  the  source  from  virtually  the 
same  angle  and  will  experience  identical  fields  --  with  the  possible  exception  of  a  small 
nearly  constant  phase  difference,  especially  if  one  point  is  nearer  the  source  than  the 
other.   In  this  case  the  complex  vectors   A (P, , t) -expfiy (P, , t) ]   and   A  (P_ ,t) 'expfiy (P_,t) ] 
will  be  of  equal  length  and  will  move  in  unison  maintaining  their  nearly  constant  angular 
phase  difference  separation   Ay  =  y(P2,t)  ~  Y(P-./t).    The  cross-spectral  density  can  then  be 
written 

WA(P1,P2)   =   iy  •  <A2(P1)  >  -exp(iAY) 

(10.22) 
=   F   A (P1) 'exp(iAy) 

or,  since   A(P..)  =  A(P_),   we  can  write  it  symmetrically  as 


WX(P1,P2)   =  yFt(X(P1).Ft/X(P2).exP(iAY).  (10.23) 

On  the  other  hand,  if  the  points   P,   and   P_   are  very  far  apart  it  may  happen  that,  due  to 
obstructions  in  the  source,  different  oscillators  contribute  to  the  field  at  the  two  points. 
Or,  if  one  point  is  much  closer  to  the  source  than  the  other,  the  random  dance  the  vector 
A  (t) -expfiy (t) ]   is  performing  at   P,   may  be  totally  dissimilar  to  that  simultaneously 
occurring  at   P_   due  simply  to  the  finite  velocity  of  light.   In  such  cases  the  phase 
Ay(P,,P2,t)   of  the  product 


A(P1,t) •exp[-iY(P1,t)] -A(P2,t) -expfiy (P2,t) ] 


=  A(P1,t) -A(P2,t) •exp[iAY(P1,P2,t)] 


(10.24) 


is  a  completely  random  angle  so  that  the  product  complex  vector  is  continuously  pointing  in 
different  random  directions.   The  average  value  of  such  a  randomly  pointing  vector  is  zero 
and  we  have 

WX(P1'P2)  =  °-  (10.25) 

These  represent  the  two  extreme  values  for  the  cross-spectral  density  function.   Usually 
it  lies  between  these  values  and  we  may  write 


where 


VW       =     VFt,A(Pl)*Ft,A(P2)   '^W  (10-26) 


<V*(P.)-V(P_)> 
pIPwPJ       =      =  (10.27: 


1"2' 


AA-V^t7x(Pl)*Ft,A(P2) 
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is  called  the  complex  degree  of  coherence  at  wavelength   X   and  at  points   P,   and   P_. 
When   |vi(P.,P2)|  =  0   the  field  is  said  to  be  incoherent  at  points   P.   and   P_ .    This 
implies  that  knowledge  of  the  phase  of  the  field  at  either  point  provides  no  useful  informa- 
tion about  the  phase  at  the  other  point.   When | y (P, ,P_) |  =  1   the  field  is  coherent  at  these 
two  points,  and  knowledge  of  the  phase  of  the  field  at  either  point  permits  the  phase  at  the 
other  point  to  be  predicted  exactly. 

When  the  field  is  so  incoherent  that  the  only  point  pairs  for  which   W, (P, ,P~)   does  not 
vanish  are  those  within  a  very  short  distance  of  one  another  (a  wavelength  or  two)  —  much 
closer  than  any  other  distance  of  physical  interest  such  as  slit  widths  —  then  it  is  common 
to  refer  to  the  field  as  completely  incoherent  and  to  write  as  an  approximation1 

WX(P1,P2)   =   X2-Aq2(P1) •6(x1-x2)«6(y1-y2)  (10.28) 

where   P,   and   P_   are  both  taken  to  be  in  the   X-Y   plane  with  coordinates   x. ,  y,   and 
x„,y_.2   Equation  (10.28)  is  assumed  to  apply  at  the  surface  of  an  incandescent  solid,  for 
example. 

As  another  example  of  the  cross-spectral  density  function  let  us  consider  a  laser  beam 
of  Gaussian  cross-section  incident  on  the   X-Y   plane  from  a  direction  in  the   X-Z   plane 
making  an  angle  of   8    with  the   Z-axis.   The  essential  geometry  is  shown  in  figure  10.2. 
In  the  plane   z  =  0   the  field  will  be  given  as  a  function  of   x   and   y   by: 

2     2     2 

y   +  x  'cos  8  \       /  27Tix«sin( 


V(x,y)   =   A  •  exp  I-  5 I  •  exp 


4d 
exp(iy)  *  exp  (-27iivt)  , 


2       /     M     X 


(10.31) 


:Any  dependence  of  the  cross-spectral  density  function  or  any  other  radiometric  quantity  upon 
position  —  x,  say  --  can  be  construed  as  being  caused  by  the  interference  of  plane  waves 
propagating  with  a  component  of  their  motion  in  the   x   direction.   If  the  cross-spectral 
density  function  changes  significantly  in  a  distance  I      then  there  must  be  plane  wave  com- 
ponents present  in  the  beam  with  wavelengths   X  <  SL      in  order  to  produce  such  spatial  struc- 
ture.  The  use  of  a  delta  function  in  the  description  of  incoherent  light  is  not  to  be  taken 
as  suggesting  that  wavelengths  approaching  zero  are  present  but  rather  as  a  convenient  mathe- 
matical approximation  valid  if  the  dimensions  I      of  features  of  interest  are  very  large 
compared  to  the  wavelength  of  light  being  observed. 

^(x-.-x-)   is  the  Dirac  delta  function  [10.11].   It  has  the  properties: 

6(x,-x?)  =  0   for   x,  ^  x_ 
and  (10.29) 


/ 


x.  +a 


x.  -a 


(x,-x_)'dx_  =  1   for  any   a  >  0. 


Since   6(0)  =  »  we  see  that  eq.  (10.28)  requires  that   F,  ,  (P)  =  W,  (P,P)  =  °°,   an  obviously 

t  ,  A  A 

unrealistic  feature  of  this  approximation.   We  should  try  to  think  of   A«6(x,-x_)   and 
X'6(y,-y_)   as  convenient  mathematical  approximations  for  sharply  peaked  functions,  of  width 
about   X   and  height  about   1.    Then  we  can  make  the  approximate  indentif ications : 

Ffc  X(P)  ~  A2(P)   and   y(P1#P2)  ~  X2 • 6 (Xj-Xj) • 6 (y1~y2 ) •  (10.30) 
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Figure  10.2.  A  train  of  plane  waves  incident  on  the  X-Y  plane 
at  an  angle  6  from  the  Z-axis  in  the  X-Z  plan 
and  with  a  Gaussian  amplitude  distribution. 


where   d   is  a  measure  of  the  diameter  of  the  beam  and  y      is  the  phase  of  the  wave  at  the 

2itx*  sin6 


origin  at  time  zero.   The  quantity 


X 


is  the  phase  difference  between  a  point  x,y 


and  the  origin  expressed  in  radians.   Assuming  that  the  laser  is  monochromatic  we  obtain  from 
eq.  (10.14) 


[2      2       2      2      2  ~\ 
yl   +  y2   +  <xl   *  x2  ,'coa  °o 
7? J 


2      2       2      2      2 

i  j±     +  y2  +  (*i     +  x2  ) *cos  ' 
Wx(x1,y1;x2,y2)   =  AQ   •  exp | 


exp 


r27ri«  (x.  -x-  )  •  sinO  "1 


(10.32) 


=  exp 


r2fri«  (x,  -x_)  •  sin6  "I     

L 1 ^J  *  VFt,x(xi,yi)*Ft,x(x2'y2)' 


where   A   =  A/ VAX. 
o 


By  comparing  eqs.  (10.32)  and  (10.26)  it  can  be  seen  that  the  complex  degree  of  coherence  is 


exp 


2iri-  (x.-x2) 


•sinO  "1 

°] 


(10.33) 


and  again   |u|  =  1.    In  this  example  the  spectral  fluence  rate  across  the  Gaussian  beam 
cross-section  is  given  by 
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Ft,A(x'y)   =   Ao 


exp 


2     2     2 
y   +  x  'cos  ( 


2d' 


[W'm  2] 


(10.34) 


and  the  phase  difference  is 


AY(x1,y1;x2,y2) 


2tt«  (x,  -x_)  •  sine 
12      o 


[rad] 


(10.35) 


The  PROPAGATION  of  CROSS-SPECTRAL  DENSITY.    If  the  cross-spectral  density  function  is  known 
at  all  points  in  some  cross-section  of  a  beam  it  can  be  calculated  for  any  other  cross-section, 
The  law  of  propagation  of  cross-spectral  density  is 


WA(Pi'P2)   = 


VP1'P2> 


exp[2iri«  (s,-s2)/A]  «cos6,  -cosS- 


sl*s2'r 


dP  -dP- 


(10.36) 


In  this  equation   W,  (P, ,P~)   is  assumed  known  for  all  pairs  of  points   P,   and   P-   in  one 
cross-section  of  the  beam  and   s.   and   s2   are  the  optical  path  lengths   P  P'   and   P?P' 
to  the  points   P,'  and   Pi   in  a  second  cross-section  where  the  cross-spectral  density  is 
wanted  (see  figure  10.3).   We  will  treat  these  two  cross-sections  of  the  beam  for  convenience 
as  a  source  plane  and  a  detector  plane  although,  of  course,  there  need  not  be  a  physical 
source  or  detector  in  either  place.   The  angles   6,   and   9_   are  the  angles  between  the 
normal  to  the  source  plane  and  the  directions  to  points   P'   and   P'   from   P,   and   P2 , 


SOURCE 
PLANE 


Figure  10.3.    The  propagation  of  cross-spectral  density  W, (P,,P2), 
assumed  known  for  all  pairs  of  points   P,   and 
P2   in  one  cross  section  of  a  beam  (source  plane) . 


Optical  path  lengths   s,   and 


ponding  points   P' 


1 
and   PJ 


,   extend  to  corres- 
in  a  second  plane  where 


the  cross-spectral  density  is  wanted  (detector  plane) 
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respectively.   In  the  integrations,   P, 
ently;  that  is 


and 


P„   each  cover  the  entire  source  plane  independ- 


dP. «dP_ 
1    2 


source 


jf    ||    dx1.dy1.dx2.dy2 


(10.37; 


where   x, ,  y,   and   x_,y2   are  the  coordinates  of   P,   and   P_   in  the  source  plane. 

The  factors  appearing  in  eq.  (10.36)  are  easily  understood.   The  exponential  is  just  the 


phase  change  resulting  from  the  difference  in  the  path  lengths   ptp{   and  P2P2 "    This  dif- 
ference is  divided  by   A   to  obtain  the  number  of  periods  and  then  multiplied  by   2it   to 
express  the  phase  difference  in  radians.   The  cosine  factors  express  the  reduction  in  pro- 
jected size  of  the  elements  of  area  dP  =  dx#dy  when  viewed  from  the  direction  of  P'.    The 

distances   s,   and   s~   in  the  denominator  are  in  effect  a  statement  of  the  inverse  square 

2 

law.   The  factor   A    m  the  denominator  xs  required  for  dimensional  homogeneity  and  to  assure 

the  correct  propagation  of  fluence  rate  computed  from  eq.  (10.36)  in  free  space.   Since,  by 
Huygens  principle,  all  points   P   can  contribute  to  each  point   P',   the  integration  covers  the 
entire  source  plane  for  both   P,   and   P_.    Equation  (10.36)  is  valid  for  source  dimensions 
and  distances   s,   and   s?   which  are  large  compared  to  the  wavelength  of  light. 

It  is  in  the  evaluation  of  the  integrals  of  eq.  (10.36)  that  the  difficulties  of  the 
treatment  of  coherence  occur.   To  illustrate  the  kinds  of  approximations  which  are  employed  to 
render  these  integrals  tractable  and  to  better  illustrate  the  meanings  of  the  factors  appearing 
in  the  integrals  let  us  evaluate  eq.  (10.36)  for  a  simple  situation.   In  figure  10.4  we  show  a 
point   P   in  a  source  plane  where  the  cross-spectral  density  is  assumed  known  and  a  typical 
field  point  P'   where  the  cross-spectral  density  is  desired.   The  'source'  plane  need  not  be 
a  real  source,  of  course.   It  can  be  any  plane  through  a  beam  anywhere.   We  take  the  source 
plane  to  be  the   X-Y  plane  and  the  origin  to  lie  in  this  plane  near  the  center  of  the  beam. 


P(x',y',z') 


Figure  10.4.    A  typical  field  point  V  (x' ,y' ,z')      is  shown  at  a 

distance   R   from  the  origin  and  a  distance   s   from 
a  point   P(x,y,0)   in  the  source  plane   (X-Y  plane). 
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The  first  task  is  to  express  the  distance   s   from   P   to   P '   in  terms  of  the  variables  of 
integration  --  the  coordinates   x   and   y   of   P.   This  is  simply 


s   =   V  (x'-x)2  +  (y'-y)2  +  z'2 
=  y   R2  -  2x'x  -  2y'y  +  x2  +  y2 


(10.38) 


V 


R  •  \   1  -  2(x'x  +  y'y)/R2  +  (x2  +  y2)/R2 


2      2      2      2 
where   R   =  x'   +  y'   +  z1    and  the  primed  coordinates  are  the  coordinates  of   P1.    We 

now  assume  that  the  source  dimensions  are  small  compared  to   R.    That  is   x/R  <<  1   and 

y/R  <<  1.    Then  we  can  expand  the  square  root  in  eq.  (10.38)  and  discard  high  order  terms 

to  get: 


s  *  R  -  [(x'x  +  y'y)/R]  +  h(x2   +   y2)/R. 


(10.39) 


Next  we  note  that  only  a  small  error  is  introduced  by  replacing   s. 


and 


by   R, 


and   R- 


in  the  denominator  of  the  integral.   A  similar  substitution  in  the  argument  of  the  exponential 
is  not  permissible  because  this 
This  is  easily  seen  by  writing 


is  not  permissible  because  this  function  changes  quite  rapidly  with  small  changes  in   s,-s_. 


exp[2Tri-  (s,  -S-)  /A]   =   cos  [  2  tt  •  (s,  -s_  )  /A]  +  i- sin[2Tr«  (s. -s_)/A] 


(10.40) 


A  change  of   s,-s?   by  as  little  as  one  wavelength  will  drive  the  trigonometric  functions 
through  an  entire  cycle  while  the  same  change  affects  the  product   s.-s-   in  the  denominator 

r  L  Z 

by  only  a  part  in   10    (assuming  s,       and   s^   are  of  the  order  of  a  meter).   So,  making 
these  approximations  and  also  approximating   cos6   by   z'/R/   which  is  consistent  with 
x,  y  <<  R,   we  obtain 


WA  (x[>y'1>z["'x2'y2'z2' 


2    £ 
R1R2X 


•    exp 


exp 


i_l_   -expp^i    (Wr    fijf      WA(xl'^l;x2'^ 
2  source 

[  2^i  /xixl +  ^yj 
L"  "^  "       Ri 


X2X2    +    Y2Y2 


[/      2    J          2              2    ^         2  \-| 
ttx   (xi  +  yi     x2  +  y2  )    .     .     .     . 
—  '\ R^ R^ )\  •dx1-dy1.dx2.dy 


(10.41) 


In  this  expression  the  integration  limits  describe  the  cross-section  of  the  beam  or  the  area 
of  the  source  and  we  recall  that  we  have  assumed  that  this  is  very  small  compared  to  the 
distance  to  the  points   P'   and   P~ . 

2      2 
The  last  exponential  factor  in  eq.  (10.41)  involving  terms  in   x,  ,   y,  ,   etc.,  identi- 
fies this  integral  as  being  akin  to  diffraction  integrals  in  the  Fresnel  approximation.   In 
coherence  calculations  this  factor  is  frequently  set  to  unity  —  consistent  with  the  stronger 
assumption  that  the  source  dimensions  are  small  compared  to   VAR.    With  this  further  simpli- 
fication eq.  (10.41)  becomes 
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Wx(xi,y',2i;x',y',z')   =   Jg*   -expf™  •  (R^)]  •  [|ff  Wx  (x^y^x^) 


1  2  source 


(10.42) 


exp 


[  2tt1  /xlxl  +  Ylyl    X2X2  +  y2y2\l 

["  HT  (      Rl      "       R2    2jJ'dx1-dy1'dx2.dy 


which  is  in  the  Fraunhofer  approximation.   This  is  the  expression  we  will  use  in  all  our 
examples  although  it  should  be  noted  that  if  W, (x, ,y, ;x_ ,y_)   is  so  complicated  that  the 
integrals  must  be  evaluated  numerically  anyway  then  consideration  should  probably  be  given  to 
including  additional  terms  or  to  calculating   (s,-s_)   exactly  using  eq.  (10.38)  wherever  it 
will  result  in  higher  accuracy.   For  the  most  part  we  can  justify  the  use  of  the  Fraunhofer 
approximation  even  when  the  source  dimensions  are  not  small  compared  to   VXr  on  the  grounds 
that  this  will  give  us  a  worst-case  estimate  of  coherence  and  diffraction  effects.   The  true 
flux  distribution  will  always  lie  somewhere  between  the  predictions  of  geometrical  optics 
theory  and  predictions  in  the  Fraunhofer  approximation  in  the  sense  that  as  the  distance   z ' 
to  the  observing  surface  increases,  sharp  edges  and  shadows  predicted  by  geometrical  optics 
gradually  become  fringed  and  fuzzy,  evolving,  for  sufficiently  large   z'   or  R,   into  the 
Fraunhofer  irradiation  pattern.   This  is  a  smoothly  periodic  pattern  dependent  upon  the 
point  of  observation  only  through  the  angular  parameters   x'/R  and  y'/R.    If  a  calculation 
in  this  approximation  reveals  that  diffraction  will  cause  no  significant  problems  in  a  meas- 
urement situation  then  a  more  exact  computation  generally  is  unnecessary.   This  will  frequently 
be  the  case  if  the  measurement  situation  requires  only  ratios  of  measurements  on  similar 
systems.   However,  if  an  accurate  estimate  of  the  loss  of  flux  from  a  geometrical  beam  is 
required  as,  for  example,  if  the  cross-sectional  area  of  the  beam  doesn't  cancel  out  in  the 
simple  geometrical  ray  analysis  of  the  measurements,  then  it  is  likely  that  a  more  exact  cal- 
culation [eq.  (10.38)]  will  be  necessary.1 

Equation  (10.42)  is  about  as  far  as  eq.  (10.36)  can  be  simplified  without  some  knowledge 
of  W,  (P. ,P_).    Included  among  the  next  twenty-five  pages  or  so  are  some  examples  in  which 
expressions  such  as  eq.  (10.42)  are  further  evaluated  in  considerable,  tedious  detail.   These 
lengthy  passages  are  set  in  a  smaller  type  to  permit  the  reader  to  identify  them  and  skip  over 
them  if  he  wishes.   We  present  such  examples  here  to  show  that  the  propagation  law  of  cross- 
spectral  density  is  not  just  another  pretty  theory  but  that  numbers  can  be  put  in  and  useful 
results  obtained.   The  examples  will  also  indicate  how  the  coherence  state  can,  in  certain 
circumstances,  significantly  affect  the  results  of  a  radiometric  measurement  when  high  spectral, 
spatial,  or  temporal  resolution  measurements  are  attempted. 

As  the  first  example  of  the  use  of  eq.  (10.42)  we  will  apply  it  to  one  of  the  simplest  possible 
sources  —  a  source  in  the  plane  z  =  0  that  is  uniform,  rectangular,  and  incoherent  with  [from 

eq.  (10.28)] 

Wx(x1,y1;x2,y2)   =  A2Aq2  •  6^-x^    •  6(y1~y2)  (10.43) 


for  and 


-b  <  y.  <  b 


-a  <  x„  <  a  -b  <  y2  <  b 


and 


W  (x  , y  ;x2,y  )   =   0   otherwise. 


xSee  reference  [10.12]  for  an  example  of  such  calculations. 
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Then  eq.  (10.42)  will  allow  us  to  calculate  the  cross-spectral  density  for  any  pair  of  points  P' 

2       2     2     2 
(at  x*   y'   zj)   and  P^   (at  x^,  y£,  zp   provided  that  both  EL,   (  =  x|  +  y|  +  z|  )   and 

2       2     2     2  2    2 

R_   (  =  x'   +  y'   +  z'  )   are  much  greater  than  a  +  b  .    If  we  substitute  expression  (10.43) 

into  eq.  (10.42)  we  can  immediately  carry  out  the  integrations  over  x„   and  y„   to  give 


wA^xi'yi'zi;x2'y2'zP 


Rl 


1Z2    .  2      pni   ,_   _  ,"1 


(10.44) 


a   b 


f     T  2*i  (Xrl  +  Vyi   X1X2  +  yly2\1  .    . 
_J  exp|-—  [ g- j- jj-dx1.dy1. 


The  double  integral  in  eq.  (10.44)  factors  into  two  similar  elementary  integrals: 


a        b 


exp 


-a     -b 


[-  ¥•  (^  -  ^H}dv  j exp [-  ¥ •  (^  -  |H. 


.dYl 


=     2a- 


faa     M       X2\l    I           I      .    r2,b     /yl        y2\" 
: — ; —i )  *2b*  (  -—, — t. 


2ira 

\ 


2TTb 


R,        R, 


(10.45) 


It  turns  out  that  these  integrals  are  among  the  few  integrals  likely  to  arise  in  these  problems 
which  can  be  expressed  in  a  closed  form  and  can  be  readily  visualized.   As  a  result,  coherence 
problems  which  lead  to  such  integrals  are  among  the  few  which  are  ever  carried  through  to  com- 
pletion and  so  the  functional  form  sin(u)/u  has  achieved  an  importance  which  is  probably  richly 
undeserved.   It  has  a  special  name 


sinc(u)   =  sin(TTu)/ (iru) 


(See  figure  10.5) 


If  we  employ  this  notation  eq.  (10.44)  becomes 


Vxi'yi'zrx2'y2>zP    ■       70  '  Ao2  •  exp[¥  -(Ri  "  V] 


R1R2 


(10.46) 


4ab-sinc 


|~2a     /!i       fiV 
LX       \R1        V- 


2b 

\        VR 


v2/-l 


Let   us   first   consider   the  case  where  both  points     P'      and     Pi      are  equally   far   from 


the   source   so   that 


R^    =   R     =   R.         Then   eq.     (10.46)    simplifies   a   li 


ittle    to   give 


z'z' 
Wx(x^y|,z|;x^y^Zp     =     4ab-   Aq2    •  -^ 

R 


(10.47) 


2a(x|  -  x'2) 

.  AR 


2b(y[  -  y') 
.      '       AR 
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Figure  10.5.    The  sine  function:   sine  x  =  sin  (ttx)  /  (ttx) 
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There  are  several  interesting  points  to  observe  about  this  cross-spectral  density  function. 
First,  the  radiation  at  a  distance  from  the  source  is  no  longer  incoherent  —  there  are  no  delta 
functions  in  eq.  (10.47)    The  radiation  has  become  partially  coherent  merely  by  the  fact  of 
propagation.   For  a  very  large  source  dimension,   a,   the  sine  function  will  oscillate  and 
decrease  rapidly  as  xi'~xo   increases  (see  figure  10.5)  so  that  eq.  (10.47)  effectively  approaches 
delta  function  behavior1  for  very  large  incoherent  sources.   A  second  feature  of  eq.  (10.47)  is 
that  the  y-dependence  is  totally  independent  of  the  x-dependence  in  that,  if  the  y-dimension  of  the 
source  is  changed,  only  the  y-scale  of  the  cross-spectral  density  function  will  be  changed.   This 
is  a  characteristic  of  solutions  in  the  Fraunhofer  approximation  of  problems  involving  rectangular 
sources  or  apertures:   provided  W. (P. ,P„)   has  this  property  the  solution  factors  into  independent 
double  integrals  —  one  double  integral  involving  only  x.   and  x„,   the  other  double  integral 
involving  only  y1   and  y„ .   Another  interesting  feature  of  eq.  (10.47)  is  that  the  x-dependent 
factor  depends  only  upon  x'   and  x'   through  their  difference  x'-xl,   and  similarly  the  y- 
dependent  factor  depends  only  upon  the  difference  y-J-yA.   When  a  function  depends  only  upon 
differences  in  this  way  it  is  said  to  be  stationary  because  it  looks  the  same  in  the  neighborhood 
of  one  x',y'   point  as  it  does  in  the  neighborhood  of  any  other  point.   In  many  coherence  problems 
where  the  cross-spectral  density  function  is  not  known  the  property  of  stationarity  is  assumed  in 
the  interests  of  simplicity  and  we  see  that  this  property  can  arise  naturally  and  is  therefore  a 
reasonable  assumption.   Finally,  we  can  calculate  the  fluence  rate  from  eq.  (10.47)  by  setting 

pi   =  pi; 

1    2 

,      2 
F_  ,(P')   =  W.(x',yl,j,;x,,.y,,il)   =   4ab  •  A   •  ^p-  (10.50) 

t ,  A  A  O       _Z 

K 

where  6   is  the  angle  between  the  normal  to  the  surface  and  the  viewing  direction.   We  see  that, 
in  contrast  to  expressions  we  shall  derive  below,  there  are  no  diffraction  effects  here:   the 

fluence  rate  at  P'   depends  upon  the  location  of  P'   only  through  the  inverse  square  of  the 

2 
distance  from  the  source  and  through  the  simple  angular  factor  cos  6.   The  reason,  incidentally, 

2 
that  this  latter  factor  in  this  equation  and  in  eq.  (10.51)  below  is  cos  0   and  not   cos6   is  that 

this  source  is  not  lambertian.   Such  a  completely  incoherent  source  as  this  may  not  even  be  physi- 
cally realizable;  however,  it  is  a  useful  limiting  case  which  is  comparatively  easily  treated  and 
is  sufficiently  realistic  for  many  purposes.   We  will  briefly  treat  a  lambertian  source  later  [see, 
e.g.,  eqs.  (10.166a)  to  (10.171a)].   From  eq.  (10.50)  we  can  calculate  the  radiant  intensity  of  the 
source  [eq.  (10.16)]  as 


'One  representation  of  the  Dirac  delta  function  is 


S(x)   =  Lim(S1^xaX).  (10.48) 

In  the  limit  as  a  and  b  become  infinite  the  cross-spectral  density  according  to  eq.  (10.47) 
approaches 

2   .2  /    1  i   ,  ,  «     -  .;  ,  , ,,  • 


X  •  Aq'  •  \-^)     •«(*[  -  x£)-6(yj_  -  yp  (10.49) 

which,  without  the  factor  I  — - — I  ,   would  be  exactly  correct  [compare  eq.  (10.28)].   The  approx- 
imations leading  to  eq.  (10.41)  and  (10.42)  are  of  course,  totally  invalid  in  this  limit  and 
evidently  the  erroneous  factor  is  a  consequence  of  this  inconsistency. 
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2      2 
I,  (9)   =   4ab  •  A   •  cos  9. 
A  o 


(10.51) 


From  eqs.  (10.47),  (10.50),  and  (10.26)  we  obtain 


p(P'  P')   =   sine 


2a 


<*i 


XR 


-]' 


'2b- (y' 


aR 


n 


(10.52) 


for   the  degree  of   coherence.      The  dependence  of   the  degree   of   coherence  upon   the   separation  of 
the  points   in   the  observing   surface   is  essentially   that   given   in  figure   10.5.      The  degree  of 
coherence   is   1   for  points  very   close   together   and  approaches   zero  for  points  very  far  apart. 
Negative  values   for   the  degree  of   coherence   for   certain  point   separations  mean   that   the   fields 
at   such  a  pair  of  points   are  predicted   to  be  out   of   phase.      A  small  value   for      |u(P',P')|      means 
that    the  phase  difference  between  points     P'      and     Pi      is  not   highly  predictable:      an  element   of 
randomness   is  present. 

In   general   the  degree  of   coherence    (like   the   cross-spectral   density   function)    is   defined 
in   terms   of   the   coordinates   of   pairs   of   points    so   that   even    for   a    fixed    frequency   or  wave- 
length  it  would   require   a   7-dimensional   space   to  plot   the   absolute   value  of   this    function.      In 
order   to   extract   some  measure   of   coherence   which   is   easier   to   grapple  with,    the   terms   coher- 
ence  length,    coherence   time,    and   coherence   area   are   often   used.      These   are   all   rough  measures 
of   the   separation  between   the   two   points      P,      and 


P~      at  which   the   degree   of   coherence   is 


reduced  to  some  arbitrary  level1  such  as  %  or  1/e  --  the  exact  level  is  not  important  as 
the  terms  are  only  semi-quantitative  at  best.  Coherence  length  refers  to  a  distance  in  the 
direction   of   propagation   --   our    z-axis   --   and   could,    for   example,    be   defined    in   the   neighbor- 


hood of   any   axial   point      (0,0, z)      as   that   value,    Az 


for  which 


u  (0,0, z    +    ^Az    ;0,0,z 


%Azc) 


=     k 


(10.53) 


The    coherence    time    is   the   time   required    for    light    to   travel    the   distance      Az 


At      =    Az    /c 
c  c' 


(10.54) 


where   c   is  the  velocity  of  light.   Of  more  interest  in  most  radiometric  situations  is  the 
measure  of  the  coherence  distance,  known  as  the  spatial  coherence  interval,  in  directions 
perpendicular  to  the  direction  of  propagation.   These  coherence  intervals,   Ax   and  Ay  , 
could  be  defined  by 


and 


u  (x  +  *sAxc,y,z;x  -  ^Axc,y,z)|  =  h 


y(x,y  +  %Ayc,z;x,y  -  5sAyc,z)|  =  h. 


(10.55) 


The  area   AA  =  Ax  *Ay   known  as  the  coherence  area  is  then  the  region  about  the  point 
(x,y,z)  in  the  plane  normal  to  the  direction  of  propagation  within  which  a  high  degree  of 
coherence  exists;  the  field  outside  this  area  will  be  only  weakly  correlated  with  the  field 
at   (x,y,z)  . 

We  see  from  figure  10.5  and  eq.  (10.52)  that  the  field  from  a  rectangular  incoherent 
source  exhibits  a  coherence  interval  in  the  x-direction  determined  by 


1  In  the  case  of  functions,  such  as   sine,   which  take  on  negative  values,  the  first  zero  of 
the  function  is  often  used. 
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2a«  (JjAx 


=   0.6 


Ax 


0.6AR 


;i0.56) 


where   0.6   is  approximately  the  value  of   u   which  satisfies   sinc(u)  =  h-      The  coherence 
interval  in  the  y-direction  is  the  same  expression  with   b   replacing   a.    The  coherence  area 
is  thus 


0.36  A2R2      ,2.. 

AA   =  c *  A  /Aoj 

c       ab 


;i0-57) 


where      Aw (    =    4ab/R    )       is    the    solid   angle    subtended   by   the    source    from  the   observation   point 
(x,y,z).         In   this   case   the   values    for   the   coherence    interval    and   coherence    area   are    independ- 
ent  of      x      and      y      because      u      is    stationary.      Equation    (10.57)    says   that    for   any    source,    the 

2 
irradiation   in   a   distant   plane  will   be   highly   coherent   over   areas   at    least    as    large    as      X    /Ago. 

If   the    source   is   partially   coherent   the   areas   of   coherence  will   be    larger.       If   an   aperture   of 
this    size    is   placed   in   the    field    from  the    source    the    radiation   passing   through    it   will   be 
highly   polarized   and  will   exhibit    strong   diffraction   effects. 

Returning  now  to  eq.  (10.46),  which  is  the  general  expression  for  the  cross-spectral  density 
due  to  a  distant,  uniform,  rectangular,  incoherent  source,  where  R.  is  not  necessarily  the  same 
as     R„ ,      we   find  the   fluence  rate   still   given  by  eq.    (10.50)   but   the  degree  of   coherence   is 


u(P^,Pp      = 


[la     /Xi        X2\l  .       [2b       M        y2\] 


(10.58) 


This  is  still  a  solution  in  the  Fraunhofer  approximation  [eq.  (10.42)].   In  this  approximation 
the  spatial  structure  of  the  field,  aside  from  the  inverse  square  dependence,  is  a  function  only 
of  the  angular  position  of  the  observation  point.   This  is  evident  in  eqs.  (10.50)  and  (10.58)  in 
the  dependence  upon  P'   and  Pi   only  through  the  ratios  x'/R, ,  y.'  /R,  ,  z'/R-,,  etc.   For  this 
reason  the  field  at  two  points  P'   and  Pi  which  lie  in  the  same  direction  from  the  source  (so 
that  x'/R,  =  x'/R„   and  y-J/R-,  =  Yo/lO   will  be  calculated  to  be  coherent   (|y|  =  1)   in  this 
approximation.   It  is  unlikely  in  radiometry  that  a  higher  approximation  would  be  needed  but,  if 
it  were,  then  the  Fresnel  approximation  of  eq.  (10.41)  is  the  next  step  to  try  if  a  closed-form 
analytic  solution  is  the  goal.   For  the  present  example  of  a  uniform,  rectangular,  incoherent 
source  eq.  (10.41)  can  be  evaluated  exactly  in  terms  of  Fresnel  integrals.   Aside  from  showing 
that  the  degree  of  coherence  does,  after  all,  depend  upon  the  difference  in  radial  distances  when 
these  distances  are  too  small  to  justify  the  Fraunhofer  approximation,  the  solution  is  of  little 
interest.   It  is  too  complicated  to  allow  any  easy  intuitive  interpretation  and,  moreover,  since 
real  sources  are  never  perfectly  uniform,  rectangular,  and  incoherent  the  accuracy  enhancement 
obtained  by  employing  a  better  but  still  incomplete  analytical  approximation  would  seem  to  be  of 
marginal  practical  value. 

Another  example  of  the  propagation  of  cross-spectral  density  that  is  easy  to  evaluate  is 
obtained  by  replacing  the  incoherent  source  we  just  treated  by  a  coherent  one.   This  could,  for 
example,  be  a  small  rectangular  aperture  of  dimensions   2a   by   2b   irradiated  by  a  laser.   We 
will  assume  that  the  angle  of  incidence  of  the  laser  beam  is   6   and  take  the  plane  of  incidence 
to  be  the  X-Z  plane.   We  will  also  assume  that  the  amplitude  is  uniform  over  the  aperture  so  the 
cross-spectral  density  at  the  aperture  has  the  phase  dependence  of  eq.  (10.32)  but  with  a  step- 
function  replacing  the  Gaussian: 
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WA(x1,y1;x2,y2)   =  Aq2  -exp^  •  (x^)  -sii^  J 


-a  <  x  <  a 

-a  <  x„  <  a 

for         Z  ^      (10.59) 

-b  <  y.  <  b 

-b  <  y2  <  b 


=  0     otherwise. 


When  we  insert  eq.  (10.59)  into  eq.  (10.42)  we  obtain 

Wx(x'>y',Z';x2>y2,zp      =       -j-^.exp^-d^   -   R2)J 


RXR2X 


a     b     a     b 


f|||     AQ2    -exppf  .(xrx2).sineo]  (10.60) 


-a  -b  -a  -b 


exp 


[  2*i   /XlXi  +  ylyl   X2X2  +  y2y2\]  .-."., 
["  T  *  V ? ^ jj-dx1.dy1.dx2.dy2 


This  multiple  integral  can  be  factored  into  a  factor  depending  only  upon  the  coordinates  of  P. 
and  a  factor  depending  only  upon  the  coordinates  of  P„,   and  since  the  variables  of  integration 
are  only  dummy  variables,  the  form  of  one  of  these  factors  is  just  the  complex  conjugate  of  the 
form  of  the  other: 

WA(x[,y|,z|;x2,y2,zp   =  V*(x.[,y^,zj)  •  VX2'y2'Z2}  (10.61) 

where 

a  b 


Vx(x.,y.,z.)   =  -|l.exp(^).|  |  Vexp^f  -x-sinsj 


-a  -b 


(10.62) 


•  expL_  ^T  -(XX  ±  yy  )J  -dx-dy. 


This  factorization,  of  course,  is  a  consequence  of  the  fact  that  for  complete  coherence  the 

angular  bracket  average,  <V*(P..  )  •  V  (P_)>,   is  simply  equal  to  the  product  V*(P  )-V  (P2), 

so  that  the  propagation  of  V  (P. )   can  be  treated  independently  of  P„   and  vice  versa.   In  this 

example 


W  =  Vexp(¥  -xrsineo) 


for  -a  <  x:  <  a  and  -b  <  y.  <  b  (10.63) 

and  V  (P  )  =  0   otherwise, 

and  eq.  (10.62)  simply  describes  the  propagation  of  this  field.  In  addition  to  this  factori- 
zation the  integrals  factor  into  separate  x-  and  y-dependent  integrals  because  we  are  con- 
sidering rectangular  apertures.   Thus 
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V^x'.y'.z') 


XR 


/2niR\  _   . 


(10.64) 


where 


IT—  ^2lTi     ■    c\  /  2iri   xx'\   , 


(10.65) 


and 


VA^-exp(-^f^).dy. 


(10.66) 


These  integrals  are  similar  to  those  we  evaluated  for  the  incoherent  source,  so  we  can  write  them 
immediately  as 


I    =   2a> \/A   -sine 
x         v  o 


and 


I   =  2b-  \fk~  -sine  (-^H  • 
y         v  o       \  XR  / 


(10.67) 


Putting  these  expressions  into  eqs.  (10.64)  and  (10.61),  we  finally  obtain 


ziz2Ao2 

WA(x|,y',z|;x^,y',zp   =   ^  2°  2 


.,  2,_2     [2Tri   ,_  _  ."I 
16a  b  -expl—  .(R^R^J 


2a 

X 


h 


sin6 


2a 

A 


sine 


(10.68) 


/2by' 
sinclT-:: — I  -sine 


XR, 


2by' 


XR„ 


for  the  cross-spectral  density  at  the  pair  of  points  P' (x' ,y ' ,z')   and  P'(x',y',z')  produced 

by  a  uniform  coherent  beam  passing  through  a  2a  by  2b  rectangular  aperture  in  the  X-Y  plane 

with  angle  of  incidence  6  .    The  spectral  fluence  rate  at  the  point  P' (x' ,y ' ,z ' ) ,  obtained  by 
setting  x'  =  x'  =  x' ,  etc. ,  is 


Ft)A(X\y',2': 


z'2    u  2,2   .  2    .2 
.  „  •  16a  b  •  A    *smc 
4,2  o 


R  X 


IMf—  •.)]■—  2(¥)-      <»•«> 


If  we  use  this  expression  to  evaluate   \/F<-  a(Pi')   ana  \^+    \  (P?)      and  then  calculate  the 
complex  degree  of  coherence  from  eq.  (10.26)  using  eq.  (10.68)  we  find 


w(P|,P2)   =   expp^:  .(Ri  _  r2)J 


(10.70) 


and 


jy(P^,P2)|  =  1. 


(10.71) 
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The  radiation  is  still  completely  coherent.   The  only  way  to  reduce  the  degree  of  coherence  is  by 
somehow  physically  changing  the  phase  relationships  so  that  some  averaging  will  take  place  in  the 
angular  bracket  expression  of  eq.  (10.14).   Since  this  average  is  over  both  area  and  time,  the 
degree  of  coherence  can  be  reduced  by  randomizing  the  phases  in  either  of  these  radiometric 
dimensions.   If  the  instrumental  spatial  resolution  is  not  too  fine,  coherent  light  scattered 
from  a  matte  surface  will  be  interpreted  as  incoherent  since  many  speckles1  will  be  included  in 
the  instrumental  resolution  area.   On  the  other  hand  if  the  instrumental  spatial  resolution  is 
finer  than  the  speckle  size  but  its  temporal  resolution  is  not  too  fine  then  motion  of  the 
scattering  matte  surface  will  also  introduce  incoherence  into  the  measurement  by  averaging  over 
many  speckles  due  to  speckle  motion  during  the  measurement  time  interval.   In  either  case,  with 
better  instrumental  resolution  the  radiation  would  be  found  to  be  completely  coherent  —  although 
if  the  resolution,  required  to  detect  the  coherence,  results  in  an  unmeasurably  weak  signal  it 
serves  no  useful  purpose  to  think  of  it  as  coherent. 

The  spectral  fluence  rate  of  eq.  (10.69)  describes  the  diffraction  pattern  associated  with 

a  rectangular  aperture.   The  pattern  is  centered  at  x'  =  R-sin6  ,  y'  =  0.    It  is  useful  to 

write  the  argument  of  the  sine  function  of  eqs.  (10.68)  and  (10.69)  in  terms  of  small  angular 

displacements  from  the  beam  direction  9  .    If  we  write 
r  o 

sine  =  x'/R 

and       a  =  9  -  9 
o 

we  have     l^'Yf-  -  sine  \  -  |*-  -(sine  -  sine  ).  (10.72) 

A     \  R  O  /        A  O 

We  can  expand   sine   in  a  Taylor  series  about   9   to  permit  approximating  sin6  by 


sine  =  sin9  +  (9-6  )>cos9  +  .  . 
o       o      o 


sin6  +  a-cos6  . 
o        o 


(10.73) 


Thus  eq.  (10.72)  becomes 


2a 

A 


2a«cos9 


I""    sin6o)     ■  \       °    '   a  (10-74) 

and   the   fluence  rate,    eq.    (10.69)    can  be  written 

,2  „  /2a- cose  \  _/_  \ 

Ft   A(a'6)   =  ~TT  '    16a  b      *A0      'sinc     \ 1 °    -aj-sinc   (jp-  •  3 ) ,  (10.75) 

R   A 

where  we  have  defined  the  angle   B  by  8  =  y'/R.    In  this  form  the  equation  shows  that  the 
angular  distribution  of  radiation  is  the  same  as  if  an  aperture  of  dimensions  2a«cos6   by  2b 
were  irradiated  at  normal  incidence.   Thus,  at  least  for  small  angular  deviations  from  the 
undiffracted  beam,  we  may  replace  tilted  apertures  by  their  projections  normal  to  the  beam. 


Speckle  is  the  random  interference  pattern  produced  by  coherent  radiation  reflected  from  a 
surface  which  is  rough  over  distances  of  wavelengths.   Individual  spots  of  the  speckle  pattern 
differ  from  one  another  in  electric  field  direction  (which  we  ignore  in  this  scalar  theory) , 
amplitude,  relative  phase  of  the  sinusoidal  electric  field  vibration,  and,  of  course,  size  and 

shape. 
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MODELS  for  the  CROSS-SPECTRAL  DENSITY  FUNCTION  —  QUASIHOMOGENEITY.    In  principle  it  is 
possible  to  measure  the  cross-spectral  density  function.   Any  radiometric  measurement  provides 
some  information  about  it.   For  example,  the  expression  [eq.  (10.20)]  for  the  spectral  radi- 
ance as  a  function  of  cross-spectral  density  can  be  inverted  to  give 

2tt  tt/2 


Wx(x+5/2,y+n/2;x-£/2,y-n/2)  =       Lx (x,y, 6 , $)  (10.76) 

0    0 

•  exp[27ri«sin6  •  (£cos<t>+nsin<}>)  /A]  •  sine  *d6  'd<\> . 

This  shows  how  spectral  radiance  measurements  can  be  employed  to  calculate  the  cross-spectral 
density.   All  that  is  needed  is  to  measure  the  spectral  radiance  at  a  point   (x,y)   from  a 
great  distance1  as  a  function  of  direction   (6,<}>)   from  that  point  and  then  evaluate  the 
integrals  for  any  desired  pair  of  points  symmetrically  placed  around   (x,y).   By  varying  E, 
and   n   all  such  pairs  of  points  can  be  covered.   Finally,  by  remeasuring  at  other  points 
(x,y)   the  entire  field  can  be  mapped  out.   From  a  practical  point  of  view  this  requires  an 
almost  prohibitive  number  of  measurements  and,  moreover,  because  of  the  dependence  of  the 
integrand  upon  sine   the  measurements  at  high  angles  become  most  important.   As  far  as  we  are 
aware  this  technique  has  never  been  applied.   More  commonly  the  cross-spectral  density  is 
obtained  from  the  measurement  of  more  conventional  interference  effects.   A  number  of  ref- 
erences to  such  measurements  are  given  in  part  II  of  reference  [10.5]. 

Usually,  if  the  cross-spectral  density  must  be  known,  a  reasonable  functional  form  is 
assumed  and  then  a  small  number  of  measurements  of  diffraction  or  interference  effects  will 
establish  the  values  of  any  unknown  parameters  of  the  functional  form.   We  have  already  seen  a 
couple  of  functional  forms  which  are  commonly  used  for  the  cross-spectral  density  function: 
the  delta  function  for  incoherent  radiation  [eq.  (10.28)]  and  models  for  completely  coherent 
radiation  [eqs.  (10.32)  and  (10.59)].   All  common  models  [10.3],  [10.8],  and  [10.13]  for  the 
cross-spectral  density  function  of  partially  coherent  fields  assume  that  the  complex  degree  of 
coherence  is  a  stationary  function,  depending  only  upon  the  distance  between  the  points   P. 
and  P_.   That  is,  it  is  assumed  that  u(P,  ,P~)   can  be  written2  as  a  function   g(?,n)   where 
5   and   n   are  difference  coordinates  defined  by 

5  =  x1  -  x2   and   n  -  Y1  -  Y2  •  (10.77) 

Thus  we  would  write 


VP1'P2}   =  VFt,X(Pl)-Ft,A(P2)-g(?^)-  (10'78) 

This  is  known  as  the  Schell  model.   A  common  form  for   g   is  given  by 

g(?,n)   =   sin  (up/a) /(ttp/o) 

where  (10.79) 

2    r2  .   2 

p  =  £  +  n  . 


lSee  footnote  on  page  48. 

2 In  general   u   and   g  will  depend  upon   A   as  well  as  upon  the  coordinates  of  the  two  points, 
For  simplicity  we  don't  explicitly  indicate  this  dependence  in  our  notation. 

25 


In  this  expression  the  parameter   a   is  to  be  adjusted  for  best  agreement  with  measured 
diffraction  or  interference  effects.   This  form  for   g,   with  a      taken  as   ~  A/2,   provides 
a  model  for  a  real  'incoherent'  source  which  will  be  lambertian  and  thus  will  describe  the 
field  over  wide  angles  more  realistically  than  the  delta  function  model  can.   Another  com- 
monly used  form  for  the  degree  of  coherence  is  the  Gaussian: 

g(£,n)   =   exp[-(C2  +  n2)/(2a2)]   =   exp  [-p2/ (2a2)]  .  (10.80) 

As  in  the  previous  expression  the  adjustable  parameter   o   is  approximately  equal  to  the 
coherence  interval. 

If  in  the  Schell  model,   u(P..,P_)   decreases  to  zero  with  increasing  separation  of   P, 
and   P-   sufficiently  rapidly  that   F   , (P, )  ~  F   \(po)   over  the  region  where   y   differs 
significantly  from  zero,  or,  putting  it  another  way,  if   F   , (P)   changes  very  slowly  with 
position  compared  with   u»   then,  in  the  expression  (10.78)  for  cross-spectral  density,  we 
can  approximate 


{■ 


Ft.x(pi>-Ft,x<F'2)   -  Ft,»lx-y»  d°-8i> 

where  the  average  coordinates   x   and   y   are  defined  by 

x  =  (x±   +  x2)/2    and    y  =  (y1   +   Y2)/2-  (10.82) 

Then  the  cross-spectral  density  can  be  expressed  as  the  product  of  the  spectral  fluence  rate 
by  a  complex  degree  of  coherence: 

WA(P1,P2)   =   Ft  x(x,y)-g(S,n).  (10.83) 

This  is  known  as  the  quasihomogeneous  model.   This  model  is  important  because  most  of  the 
radiometric  properties  of  a  quasihomogeneous  source  are  in  accord  with  the  ideas  of  classical 
(geometrical  optics)  radiometry.   (For  example,  radiance  is  never  negative  and  is,  to  a  good 
approximation,  invariant  along  a  ray  [10.14]).   This  is  because,  although  a  quasihomogeneous 
source  may  be  highly  coherent  locally,  the  coherence  areas  are  small  compared  to  any  bright- 
ness features  on  the  surface  so  that  the  source,  taken  in  its  entirety,  in  many  respects 
still  behaves  like  a  classical  (incoherent)  source. 

Another  useful  feature  of  the  quasihomogeneous  model  is  that  if  a  source  is  quasihomo- 
geneous then,  the  distant  field  from  that  source  is  also  approximately  quasihomogeneous.   To 
show  this  let  us  consider  eq.  (10.42).   We  must,  however,  remark  before  doing  this  that  the 
validity  of  the  Fraunhofer  approximation  is  even  more  strained  here  than  in  many  coherence 
calculations  because,  in  addition  to  requiring  that  \/A~R  be  large  compared  to  the  source 
dimension,  we  require  that  the  source  dimension  be  large  compared  to  the  coherence  interval 
which,  itself,  can  be  taken  as  larger  than   A   (otherwise  we  could  simply  treat  the  source  as 
incoherent) .   We  will  assume  that  the  points   Pi   and  PI   lie  near  the  z-axis  on  the  surface 
of  a  large  sphere  centered  at  the  origin  so  that  we  may  take   R,  =  R»  =  R   and   z'  =  z'  =  z'. 
Then  we  solve  eqs.  (10.77)  and  (10.82)  for   x,   and   x_   in  terms  of   x   and   £,   etc. 

x±   =  x  +  5/2   yx  =  y  +  n/2 

(10.84) 
x2  =  x  -  C/2    y2  =  y  -  n/2 
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and  substitute  into  eq.  (10.42),  making  analogous  substitutions  for  the  primed  variables1 


VPi'P2>   =   ^2 


wx(x+^c,  y+'sn;  x-hi,   y-'sn) 


source 


,2 


R4A2 


iXpL^XIf  '  (?x,+?,x+ny'+n'y)]  -dx-dydS-dn 
|f   Ft#x(x,y)  -exp[-^i.  (C'x+n'y)]  -dx-dy 


(10.88) 


jf  g(C,n)  -exp[-^i.  Ux'+ny')]  -d?-dn 


source 


where  we  have  used  eq.  (10.83)  for  the  explicit  representation  of  the  quasihomogeneous  source. 
The  limits  of  integration  in  the  space  of  the  average  and  difference  coordinates  are  different 
than  in  the  x, ,  y, ,  x_ ,  y_   space.   For  example,  for  a  rectangular  source  where  the  limits  on 
x 


and   x_ 

to   2a-2|x| 


are   ±a   we  could  have  [from  eqs.  (10.77)  and  (10.82)]  E,      run  from   2|x|-2a 
and   x   run  from   -a   to   +a.    For  a  quasihomogeneous  source  the  range  over 
which   g   differs  significantly  from  zero  is  so  small  compared  to  the  source  dimensions  that 
no  appreciable  error  is  introduced  by  replacing  the  E,      and   n   integration  limits  by  the 


same  limits  that   x. 


and 


y,   are  subject  to  or  even  possibly  by  ± 


other  words,  in  comparison  to  the  size  of  the  source, 


if  convenient.   In 
resembles  roughly  a  delta  function 


*It  is  not  obvious  that  the  elements  of  integration  are  related  by 


dx 


,  •dy-.'dx-'dy-   =   dx'dyd£-dr|. 


(10.85) 


This  equality  follows  from  the  fact  that  the  Jacobian  of  the  coordinate  transformation  given 
by  eqs.  (10.77)  and  (10.82)  is   1.    The  Jacobian  of  a  transformation  such  as  this  is  the 
determinant 


3 (x1,y1,x2,y2) 
3(x,  y,  5,  n  ) 


axj 
3x~ 


3xn 


3xn 


3x 

3x. 
3x~ 

Hi 

3x 


3y    3C 


3yi  3yi 


3y 

3x, 


*Jl 

3? 

3x, 


3y   35 
3y2   3y. 


3x. 

3n 

**1 

3n 

3x2 

3n 

3y2 

3n 

[10.86) 


This  expresses  how  the  multidimensional  element  of  integration  transforms  under  the  coordinate 
transformation,  i.e., 


dx. *dy, -dx-'dy- 


3 (x1,y1,x2,y2) 
3(x,  y,  ?,  n  ) 


■dx'dydC'dn . 


(10.87) 


[For  more  details  see  any  text  on  advanced  calculus.] 

27 


so  that,  beyond  a  certain  point,  the  limits  of  integration  of  any  integral  involving  it  are 
largely  irrelevant.   In  our  rectangular  source  example,  then,  we  could  take  the  limits  of  the 
E,      integration  also  to  be   ±a,   for  instance,  without  significant  error.   Allowing  further 
approximations  of  this  sort  then,  we  see  that  the  left  hand  side  of  eq.  (10.88)  is  a  product 
of  two  independent  functions  defined  by  double  integrals:   the  first,  a  function  only  of  the 
difference  coordinates   £ '   and   n"   in  the  observation  surface,  and  the  second,  a  function 
only  of  the  average  coordinates   x'   and   y1   in  the  observation  surface.   So  we  evidently 
can  write 

WA (x|,y|,z' ;x£,y£,z')   =   F£  A (x' ,y ' ) -g' (5 ' ,n ' ) ,  (10.89) 


where 


g*  11', r)')      =     N 


1_ 

X 


Ft>x(x,y).exp[-|^-  U'x+n'y)]  -dx-dy, 


source 


F^x(x',y')   =  nx-/5_!_\  .  ||  g(C,n).exp[-^i.  (Cx'+ny')] -dC-dn,  (10.90) 


source 


and  Nx   =     ''   j   Ft  x  (x,y)  "dx*dY' 


This  shows  that  the  cross-spectral  density  in  the  observation  surface  is  also  approximately 
quasihomogeneous  at  least  within  the  validity  of  the  Fraunhofer  approximation.   In  deriving 
this  result  we  assumed  that   R.  =  R_  =  R   and   z,'  =  zl  =  z*.    We  assume  that  these  will 
still  be  good  approximations  and  our  conclusions  valid  for  field  points  near  the  z-axis  in 
any  plane  parallel  to  the   x-y  plane  provided  the  distance,   z',   from  the  source  is  large 
compared  to  the  coherence  intervals  there.   The  field  from  a  quasihomogeneous  source  will  not 
be  quasihomogeneous  beyond  an  aperture  which  is  not  large  compared  to  the  coherence  area  at 
that  aperture.   At  such  an  aperture  the  condition  that   F,  , (P)   be  slowly  changing  compared 

t  ,  A 

to   y   is  not  fulfilled  and  this  aperture  then  constitutes  a  non-quasihomogeneous  source  for 
the  field  points  beyond  it. 

Equations  (10.90)  represent  a  simplified  propagation  law  for  quasihomogeneous  fields. 
We  notice  that  the  complex  degree  of  coherence  in  the  source  determines  the  spectral  fluence 
rate  and  hence  the  directionality  of  the  distant  field.   Conversely,  the  distribution  of 
spectral  fluence  rate  or  radiance  in  the  source  determines  the  degree  of  coherence  in  the 
field.   Frequently  only  the  second  of  eqs.  (10.90)  will  be  required  and  this  results  in  a 
further  simplification  of  the  propagation  law.   For  example,  a  rectangular,  uniform,  incohe- 
rent source  [eq.  (10.43)]  is  quasihomogeneous.   If,  within  the  rectangle  of  the  source,  we 
identify 

Pt  A   =  Aq   and   9  =  a2-6(5) -6(n).  (10.91) 

Then,  from  the  second  of  eqs.  (10.90),  we  immediately  obtain  eq.  (10.50)  for  the  spectral 
fluence  rate  at  a  distance  from  such  a  source: 
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F^A(x-,y)  =  v(^2.a2 

(10.92) 

=   4ab-A2-cos26/R2. 
o 

We  will  not  make  much  use  of  the  quasihomogeneous  approximation  in  this  chapter  because 
we  will  normally  wish  to  present  calculations  covering  the  entire  range  from  incoherent  to 
coherent  to  the  same  degree  of  approximation.   However,  it  should  be  apparent  that  this  model 
is  extremely  useful  in  many  problems  involving  fields  originating  in  thermal  sources. 

The  EFFECTS  of  IMAGING  or  FOCUSING  OPTICS.    The  preceding  examples  of  the  propagation  of  the 
cross-spectral  density  function  illustrate  the  meanings  of  the  variables  and  the  kinds  of 
approximations  which  are  frequently  made  in  evaluating  the  integrals.   Usually,  of  course, 
one  is  not  interested  in  simple,  unobstructed  propagation  from  a  source  to  a  detector  but 
rather  propagation  through  a  series  of  lenses  and  apertures  and  what-not  to  a  final  detector. 
To  handle  such  a  situation  one  must  sequentially  apply  the  propagation  equation  [eq.  (10.36)] 
from  the  plane  of  each  aperture  which  might  restrict  the  beam  to  the  plane  of  the  next  one  so 
that  the  cross-spectral  density  calculated  for  one  surface  becomes  the  source  for  the  next 
calculation,  and  so  on  until  finally  the  detector  is  reached.   Clearly,  for  a  complicated 
optical  system,  evaluating  the  multiple  integrals  involved  will  become  a  chore.   Fortunately, 
one  of  the  most  important  optical  systems  —  a  simple  source-lens-image  system  —  is  not 
significantly  more  difficult  than  the  surface  irradiation  examples  we  have  just  treated  so  we 
will  turn  now  to  that  problem. 

We  will  consider  the  propagation  of  the  cross-spectral  density  from  the   X-Y  source 
plane  to  a  spherical  surface  of  radius   R   centered  at  the  origin  (figure  10.6a).   We  can 
think  of  this  surface,  for  example,  as  a  spherical  mirror.   Then  we  will  consider  the  further 
propagation  of  the  cross-spectral  density  from  this  surface  to  the  image  plane.   The  first 
part  of  this  task  is  accomplished  by  setting   R,  =  R2  =  R   in  eq.  (10.42): 

COS9, 'COS6_     t r r r 

"Wi'*2'*V      =  —^2 *  JJJJ  Wx(x1,y1;x2,y2) 


source 


(10.93) 


•  exp[_  fir  ,(xixi  +  yj.yi  -  X2X2  -  Y2Y25]  •dx1.dYl.dx2.dy2. 

In  the  integrations,  points   PWx,,.;^)   and   P2(x2,y2)   each  independently  cover  the  entire 
source  area  in  the   X-Y   plane.   The  points   P'   and   P^,   for  which  the  new  cross-spectral 
density  on  the  left  is  calculated,  lie  in  the  spherical  surface  of  radius   R.    Actually  the 
surface  need  not  be  spherical.   All  we  require  is  that  the  optical    path    length    from  the 
origin  along  all  geometrical  optic  rays  to  the  surface  be  constant.    Therefore  this  surface 
could  be  the  plane  through  the  center  of  a  lens  (figure  10.6c).   On  the  image  side  of  this 
"spherical"  surface  we  shall  assume  the  radius  is   R" ,   thus  allowing  for  the  possibility  of 
imaging  with  a  magnification  other  than  one.   For  the  propagation  from  this  surface  to  the 
image  plane,  then,  we  start  with  eq.  (10.36)  and  approximate  the  distance   s"   (figure  10.6b) 
from  the  surface  to  points   P"   in  the  image  plane  by  [see  eq.  (10.39)] 

s"  «  R"  -  [(x"x'  +  y"y')/R"]  +  h    (x"2  +  y"2)/R"  (10.94) 
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P"(x",y",z") 


P(x,y,0) 


Figure  10.6.    (a)  Cross-spectral  density  propagation  from  source 
plane  to  intermediate  optical  element  (mirror) . 

(b)  Cross-spectral  density  propagation  from  inter- 
mediate optical  element  (mirror)  to  detector  plane, 

(c)  Cross-spectral  density  propagation  from  source 
plane  through  intermediate  optical  element  (lens) 
to  detector  plane. 
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and  finally  obtain 


Wj(xJfy",x3fy5)  =  R||2tx2    -exp^- (x»   +  y^   -  x-   -  yJJ  )J  •  jjjj  W-  (x  ■  ,y ■  ;  x  •  ,y  ■  , 

aperture 

*  SXP[-  T^'{HX'i    +    *!*!  "  X2X2  "  ^J  -dx-'dy'-dx^dy'. 


aperture 

(10.95) 


where  we  have,  as  usual,  approximated   sV*so   iR  tne  denominator  by   R"  .    In  the  integra- 
tions of  eq.  (10.95)  each  point   P'(x',y')   and   P'(x',y')   covers  the  entire  spherical 
surface  independently.   We  will  assume  that  there  is  a  rectangular  aperture  at  this  focusing 
optic  which  limits  this  integration  to  a   2a'   by   2b'   area.   The  points   P"   and   P"   of 
the  cross-spectral  density  function  on  the  left  are  points  in  the  image  plane. 

The  next  task  is  to  substitute  eq.  (10.93)  into  eq.  (10.95).   First  we  will  drop  the 
exponential  factor  outside  the  integral  in  eq.  (10.95).   The  points  of  primary  interest  in 

the  image  plane  lie  near  the  image  of  the  source  and  since  source  points  in  the  Fraunhofer 

2     2 
approximation  satisfy   (x   +  y  )/(XR)  <<  1   the  points  of  interest  in  the  image  plane  will 

satisfy: 

*"2  +,r2  =   (*R"/R)2  +(yR"/R)2  =  x2  hh  y2  (R11/R)  <K   R„/R  (10.96) 

AR  AR  AH 

where   R"/R   is  the  optical  system  magnification.   Thus,  the  quadratic  phase  factor 

exP[|i7r(x»2  +  y»2  -  x»2  -  y»2)]  (10.97) 

will  be  unity  except  for  very  large  system  magnifications  or  image  plane  points  very  distant 
from  the  Z-axis.   Second,  the  cosine  factors,   cose,,  cosO-,  cosG',  and  cos6'   will  prove  to 
be  a  nuisance  in  what  is  to  follow  so  we  will  now  make  one  final  approximation  and  set  all 
these  factors  equal  to  1.   The  angles   6,   and   62   are  the  angles  at  the  object  plane 
between  the  Z-axis  and  the  directions  to  the  points   P,'   and   Pi   in  the  spherical  surface, 
and  similarly   6'   and   61   are  the  corresponding  angles  at  the  spherical  surface  between  the 
normals  to  the  spherical  surface  at   P'   and   P'   and  the  directions  to   P"   and   P" ,   respec- 
tively.  So  replacing  these  factors  by   1   means  that  we  are  assuming  a  small  aperture  at  the 
lens  or  mirror,  i.e.,   a'   and   b'   are  small  compared  to   R   and   R" .    With  these  approxi- 
mations, and  reversing  the  order  to  integrate  first  over  the  aperture  and  then  over  the 
source,  we  obtain: 

W»(x»,y«;x«,y»)  =   2*2.^-  \\\\   Vxl'ri?x2'*2) 

source 

'  J  fflf      eXPf"    IF1'  UlXi    +    Ylyi    "    X2X2'    -    ^2)]  (10-98) 

lens 
aperture 

•    exp[-   x^-(x'x^   +   y'y^    -    x^    -    y^)]  -dx^dy  •  -dx^dy  •  j-dx^dy.,  -d^  -dy2  . 
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Thanks  to  our  having  chosen  a  rectangular  aperture  and  having  eliminated  all  the  other 
factors  which  depend  upon   x' ,  y ' ,  xl ,  and  y'   the  integrations  over  these  singly  primed 
variables  can  be  carried  out  and  (Surprise!)  lead  again  to   sine   functions: 


W^(x»,y»;x^y»)   =    2   *2   4  '  ffff   Vl^l'V^ 

source 

•  l6a^.b^.sinc[^.(^+^)].sinc[^.(^+^)] 

r2a'  (X2     ■    X2\l    .   feb'  (Y2     ,  y2\" 

'hr  •  \s- +  w)}  •  8inobr  *vr-  +  **)_ 


(10.99) 


sine 


•  dx, 'dy. *dx2*dy2. 


Thus,  the  cross- spectral  density  in  the  image  plane  is  given  by  an  integral  over  the  object 
plane  of  the  cross-spectral  density  of  the  object  weighted  by  a  product  of   sine   functions. 

Equations  (10.98)  and  (10.99)  have  been  derived  in  the  Fraunhofer  approximation  which  is 
perhaps  rarely  applicable  in  practical  imaging  situations.   Fortunately,  some  of  the  higher 
order  terms  which  we  have  neglected  actually  cancel  out  in  the  case  of  image  formation  so  the 
approximation  is  not  as  bad  as  it  might  be.   We  will  return  later  to  a  more  general  considera- 
tion of  imaging  and  discuss  a  more  exact  equivalent  of  these  equations  which,  in  principle, 
can  be  applied  to  other  aperture  shapes  and  includes  lens  aberrations  and  focusing  errors. 
The  present  approximation  is  generally  suitable  for  very  small  sources  such  as  pinholes  and 
narrow  slits  but  overstates  diffraction  effects  for  larger  sources.   For  example,  if  the 
radiances  at  the  centers  of  two  circular  sources  of  different  sizes  are  to  be  measured  by 
imaging  them  in  turn  onto  a  small  detector  located  at  the  center  of  the  image  then  the  esti- 
mate of  the  contribution  to  the  detector  signal  of  the  larger  source  caused  by  diffraction 
from  the  extra  source  size  will  be  too  large  in  this  approximation.   Usually,  in  relative 
measurements,  such  errors  will  be  small  and  calculations  in  the  Fraunhofer  approximation  will 
be  adequate  to  establish  the  negligibility  of  the  error  or  to  indicate  the  measurement  modi- 
fications required  to  make  it  negligible. 

If  the  same  derivation  is  carried  out  for  a  quasihomogeneous  source  using  the  propagation 

law  of  eqs.  (10.90)  and  assuming  that  the   2a1  by   2b'   intermediate  aperture  at  the  focusing 

optic  is  large  compared  to  the  coherence  area,  then  we  find  that  the  image  is  quasihomogeneous 
with  cross-spectral  density  given  by 

W£(x£,y|;x2r,y£)   =  wj;  (x"  +  £"/2,y"  +  n"/2;x"  -?"/2,yn  -  n"/2) 

(10.100) 
=   F^x(x",y")-g"(S",n")/ 


where 


g"(5",n")   =  |-  f|  g(?,n)-2a'.sinc[2^-(|  +  j£)]  -2b' -sine  [^--(g  +  j^-)].dS-dn,   (10.101) 
source 

*t.X*m'*m>       -   G'(^)2'(^2-)2-  jjFt,A<^>-2a' 

source 


source  (10.102) 


and 
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^   g(5,n)'2a,-sinc(^^)-2b,-sinc(^l)-d?-dn.  (10.103) 


In  these  equations   x" ,   y" ,  E," ,      and   n"   are  the  average  and  difference  coordinates  in  the 
image  plane  analogous  to  those  defined  in  eqs.  (10.77)  and  (10.82)  for  the  source  plane.   The 
condition  that  the  intermediate  lens  aperture  be  large  compared  to  the  coherence  area  at  the 
aperture  is  almost  invariably  fulfilled  in  optical  radiometry  if  the  imagery  is  good  and  the 
source  is  quasihomogeneous.   We  can  calculate  the  coherence  area  from  the  degree  of  coherence 
in  the  aperture  which,  in  turn,  is  obtained  by  eq.  (10.90)  from  the  fluence  rate  distribution 
in  the  source.   Assuming  a  roughly  uniform  source  of  dimensions   2a   by   2b   we  find  that 
g'(5'/n')   is  proportional  to 

sinc(2a£'/XR)  •  sine  (2bn  '/XR)  (10.104) 

2    2 
and   hence   the   coherence   area   is   roughly      0.4A    R  /ab.       [See   eqs.     (10.56)    and    (10.57)].      The 

2    2  2 

condition  on   the   aperture   then  becomes      4a'b'    >>    0.4A    R  /ab      or      id    >>    1.6A    /(4ab)      where      w 

is   the    solid   angle    subtended   at   the    source   by   the   aperture,    typically      0.01    [sr]       or    larger, 

and      4ab      is   the    source   area. 

We  will  now  illustrate   the  use  of  eq.    (10.99)    by  applying   it   to   the  problem  of   imaging   the 
same    two    sources  we   treated   before,    namely,    a   uniform  rectangular    incoherent    source   and   a   uniform 
rectangular   coherent    source.      We  have   already  seen    [eqs.    (10.47)    and    (10.68)]    what   the  cross- 
spectral  density  produced  by   these   sources   looks   like   in  the  plane  of   the   imaging  optics.      Now  we 
will   see  what   happens  when  this  radiation   is  brought   to   a   focus.      Since  the  equations   are  rapidly 
becoming   inconveniently   long  we  will   treat   only   the  x-component   of   the  problem,    recalling   that    in 
the  Fraunhofer  approximation  with  rectangular   sources  and   apertures,    the     x     and     y     parts   of   the 
integrals  are  entirely   separate   factors  which  can  be  evaluated   independently  of   one   another.      Thus 
we  will  write  just   the  x-component   of  eq.    (10.99)    as 


U..00/V'-  v'M  4a'2         {{     IT(x),  .         .      ["2a'     /Xl       Xl\l 

wx      (xi'x2}    -    ^T'JJ    wx    (xi;x2}   •sincL—VF-  +  lWJ 

source 

.   f"2a'   /X2   X2'V 


(10.105) 


■dx. -dx„ . 


There  is  an  analogous  y-dependent  factor  containing  b   and  y   instead  of   a   and   x   and  the 
complete  cross-spectral  density  is  understood  to  be  a  product  of  the  two  factors: 

Wx(x1,y1;x2,y2)   =  W^ (x)  (x^x^  -W^^  (y1  ;y2>  .  (10.106) 

In  other  coordinate  systems. such  as  polar  coordinates  the  integrals  do  not  factor  this  way  and  it 
is  necessary  to  start  from  the  complete  equation  (10.98). 

For  the  incoherent  uniform  rectangular  source,  we  have  from  eq.  (10.28) 


(x) 
W^   (x  ;x  )   =   XA  -6(x  -x„)       for  -a  <  x1  <  a 


for   | x  |  >  a 


(10.107) 


and  after  substitution  of  this  expression  into  eq.  (10.105)  and  carrying  out  the  integration  over 
x„   we  are  left  with 
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W^(x)(x^;x^) 


4a,2-A 
c 

RR"a 


.   [2a'   /Xl  +  XM" 

:[if.(i4)].dIl 


(10.108) 


This  integral,  if  needed,  can  be  evaluated1  but  frequently  the  fluence  rate  is  of  more  interest: 


F"(x)(x") 

t  ,  A 


W"(x)(x";x") 


4a*2. A 


RR"A 


|_  X    \R    R"/ 


(10.115) 


•dx, 


The  integral  appearing  in  this  equation  can  be  evaluated  in  terms  of  tabulated  functions  with  the 
result 


1When  the  integration  is  carried  out,  the  result  is 

A  A 


W^(x)(x£;x!p 


(L  -  C  +  S) 


where 


2  it 


\2l*L   .(x"  _  x")l 
LXR"    Ul    VJ 


(10.109) 


:os  [|p^-  -(x^  -  x!p]      /x£R  +  aR"   x^R  -  aR1 

__ll  It  ~    *   *»n   I   .11-.         "TTTT   *    _.lln    i     _  ill 


x"R  -  aR"   x"R  +  aR" 


(10.110) 


■|Ci[U2(+)]  -  Ci[U2(-)]  -  C1[U1(+)]  +  Cifu^-)]}      (10.111) 


and 


!in[fR^-(xl-X2>] 


U  (+)  =  Aia^_  (Jc   a 
V+J    A    \R"  +  R 


V  ;     A    \R"    R/  * 


•{Si[u2(+)]  -  Si[u2(-)]  +  Si[ui(+)]  -  SiflJ^-)]}      (10.112) 


(10.113) 
(10.114) 


Most  of  the  integrals  which  arise  in  diffraction  problems  with  simple  rectangular  and  circular 
sources  and  apertures  can  be  evaluated  with  the  help  of  any  of  the  references  [10.15],  [10.16]  or 
[10.17].   Reference  [10.15]  also  includes  tables  of  numerical  values  of  many  of  the  functions  which 
arise  in  these  problems,  such  as  Ci(x)   and  Si(x).    These  tabulations  are  useful  for  isolated 
calculations.   Of  probably  more  value,  however,  are  the  rational  approximations  presented  in  this 
reference  which  permit  efficient  computer  evaluation  of  these  functions  with  good  accuracy.   The 
brute  force  computer  evaluation  of  multiple  integrals  is  often  sufficiently  slow  and  uncertain  as 
to  warrant  the  expenditure  of  considerable  effort  to  find  physical  idealizations  and  mathematical 
approximations  which  permit  reduction  of  the  integrals  to  a  form  which  can  be  evaluated  in  terms 
of  well-known  functions  with  rational  approximations. 
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■-'[¥  ■  (£  -  f )]     "-'ErMS  +  S)] 

2ira'        /x"        a\  2ira'        /x"        a\ 

X  VR"   ~   R^  X         '   VR"        R/ 

x 

The  function  Si(x)  =  J   -dt   is  known  as  the  sine  integral  and  is  tabulated  in  [10.15], 

where  rational  fraction  approximations  suitable  for  computer  calculations  are  also  given.   The 
complete  fluence  rate,  of  course,  is  the  product  of  the  x-component  of  eq.  (10.116)  and  a  similar 
y-component  in  which  x  and  a  have  been  replaced  by  y  and  b.   The  quantities  a  and  b 
are  the  object  half-dimensions  and  a'   and  b'   are  the  aperture  half-dimensions  at  the  imaging 
optic.   The  quantity  R  is  the  object-to-lens  distance  and  R"   is  the  lens-to-image  distance. 
The  coordinates  of  the  point  of  interest  in  the  image  plane  are  x"  and  y".   In  figure  10.7  we 
show  some  graphs  of  eq.  (10.116)  for  several  values  of  the  lens  aperture  halfwidth  a'.    This 
figure  illustrates  how  the  image,  which  ideally  would  be  constant  from  x"  =  -aR"/R  to 
x"  =  aR"/R  and  zero  elsewhere,  broadens  by  diffraction  as  the  aperture  is  decreased.   Of  course 
we  have  assumed  no  imaging  errors  other  than  diffraction  in  eqs.  (10.98)  to  (10.116). 

To  the  extent  that  spectral  fluence  rate  approximates  spectral  irradiance  this  is  the  dis- 
tribution of  spectral  irradiance  we  would  expect  to  find,  for  example,  in  the  image  of  a  tungsten 
strip  filament  lamp. 

At  the  other  extreme  —  a  coherent  uniform  rectangular  source  —  we  have  for  the  source 
cross-spectral  density  [eq.  (10.59)] 


W^x)(xi;x2)   = 


•exp I— —  • (x  -  x„)'Sin6  I   for 


=  0     otherwise 

and  (10.117) 

wjy)(yi;y2)   =   AQ      for  |yi|,  |y2|  <  b 

=  0       otherwise. 

These  can  be  substituted  into  eq.  (10.105)  and  its  y-dependent  analog.   Again  let  us  consider  the 
spectral  fluence  rate  in  the  image  plane1,  which  is  obtained  by  setting  x"  =  x"  =  x"  and 
y"  =  y"  =  y"  and  carrying  out  the  integrations.   The  result  is 


*The  cross  spectral  density  function  is: 


W'^^x^x^1) 


o   R      f2iTi   R   ,  „    „.   .  .  ';; 

7  •  w  -exphr  •  ir  ,(xi  -  x2),sineoJ 

(10.118) 
•  [Su(x")-Su(x")  +  Cu(xp-Cu(x'2')  +  i.Su(x2')-Cu(x'p  -  i-Su(xp -Cu(x2')  ]  . 

W"  y  (y";y")   in  this  example  is  obtained  by  replacing  x",   a  and  a'   by  y",   b  and  b' 

A       1   Z 

by  setting   8   =0.    The  functions   Su(x)   and   Cu(x)   are  defined  in  eq.  (10.120). 
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Figure  10.7.    Relative  distribution  in  the  x-direction  of  fluence  rate  in  the 
image  of  a   50  [ym]-wide,  uniform,  incoherent  source  for  various 
lens  aperture  widths.   Plots  of  the  part  of  eq.  (10.116)  within 
curly  brackets  for  X   =  514  [nm],   a  =  25  [ym]  ,   R  =  R"  =  0.8  [m]  , 


and 


=  0.01  [m] ; a'  =  0.02  [m]  ; 


—  a'  =  0.04  [m] ; 


a'  =  0.08  [m] . 


36 


F'^^(x")  =  -|  •  |jr  -[su2(x")  +  Cu2(x")] 

^•«+i)]--fit-(S-!)]|: 


(10.119) 


Fi.(y)/V..x  _  _°  .  iL 

Ft,X  (y  }      2    R" 

IT 


where 


Su(x")   =  »s{Si[U(++)]  "  Si[U(-+)]  +  Si[U(+-)]  -  S1[U(  — )]} 
Cu(x")   =  h   |C1[U(4-)]  -  C±[U(— )]  -  Ci[U(+f)]  +  Ci[U(-+)]} 


(10.120) 


and 


^      -     f  '  (! 


U<+">   =  f  •(R-  +  t)-(a'  -R'^o) 


U(-+) 


U(— ) 


2tt   /: 

(: 


2tt   /  x 
R' 


!H°' 


+  R-sine 


.) 


(10.121) 


f).(a'  +R.sineo) 
|).(a'  -  R-sineJ 


In  these  expressions  the  source  has  dimensions  2a  by  2b  and  the  aperture  at  the  imaging  optic 
has  dimensions  2a'   by  2b'.    The  distance  from  the  source  to  the  aperture  is  R  and  from  the 

aperture  to  the  image  is  R" .    The  function  Ci(x)   is  the  cosine  integral,  which,  for  our 

b 

/cos  t 
— - —  -dt  and  is  tabulated  in  [10.15]. 

a 
COHERENCE  EFFECTS  on  a  MONOCHROMATOR  SLIT-SCATTERING  FUNCTION.    The  expression  for  the 
fluence  rate  distribution  in  the  image  of  a  monochromatic  rectangular  coherent  source  [eq. 
(10.119)  above]  is  of  some  interest  because  it  is  a  model  for  calculating  the  slit-scattering 
function  to  be  observed  by  laser  irradiation  of  the  entrance  slit  of  a  monochromator .   In  this 
application  the  source  plane  is  the  rectangular  monochromatically  irradiated  entrance  slit  and 
a  single  image  of  the  entrance  slit  is  formed  somewhere  in  the  exit  plane.   The  detector  then 
collects  the  radiation  passing  through  the  exit  slit  as  the  slit  is  scanned  across  this  image 
either  by  motion  of  the  exit  slit  or  by  rotation  of  the  dispersing  element.   To  illustrate 


x) 


(x") 


this  calculation  of  a  slit-scattering  function  we  show  in  figure  10.8  curves  of  F 
for  several  values  of   6  ,   the  angle  of  incidence  of  the  laser  beam  on  the  entrance  slit. 
The  calculation  gives  the  value  of  the  spectral  fluence  rate  at  each  point  of  the  exit  plane. 
This,  of  course,  would  theoretically  be  observable  only  with  an  inf initesimally  narrow  exit 
slit.   In  practice  one  usually  has  entrance  and  exit  slits  of  nearly  the  same  width  so  that  an 
observed  slit-scattering  function  would  correspond  to  a  running  average  (integral)  of  one  of 
the  curves  in  figure  10.8  where  the  width  of  the  average  would  be  comparable  to  the  width  of 
the  entrance  slit  image  —  50  [ym]   in  this  example.   The  resulting  smoothing  would  greatly 
reduce  most  of  the  finer  structure  of  the  curves.   But  it  would  not  eliminate  the  high  level 
of  the  far  wings  of  the  slit-scattering  function  when  the  angle  of  incidence  corresponds  to 
pointing  the  laser  at  the  edge  of  the  imaging  aperture  (the  edge  of  the  grating,   6o  =  0.047 
[rad]   in  figure  10.8)  or  the  splitting  of  the  slit-scattering  function  into  twin  peaks  at 
larger  angles  of  incidence.   The  values  of  the  parameters  chosen  for  this  example  correspond 
to  those  of  a  moderately  high  resolution  single  monochromator  (monochromator  focal  length 
R  =  R"  =  0.8  [m] ,   slit  halfwidth   a  =  25  [ym] ,   grating  projected  halfwidth   a'  =  0.0375  [m]  , 
and  wavelength   A  =  514  [nm] . ) 
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Figure   10.8. 


Relative  distribution  in  the  x-direction  of  fluence  rate  in  the 
image  of  a   50  [ym]-wide,  uniform,  coherent  source  for  various 
angles  of  incidence  of  the  beam.   Plots  of  the  double  integral 
of  eq.  (10.119)  for  the  same  parameter  values  as  in  figure  10.7 


except   a'  =  0.0375  [m]   and  A  =  1. 


=  0.00  [rad] ; 


9   =  0.03  [rad]; 6Q  =  0.04  [rad]; 6Q  =  0.05  [rad] 
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These  curves  show  that  the  measurement  of  the  slit-scattering  function  with  a  laser  beam 
requires  careful  attention  to  the  angle  of  incidence.   With  incoherent  light  this  problem 
doesn't  arise  because  there  is  no  particular  direction  associated  with  incoherent  light  aside 
from  the  obvious  dependence  upon  the  projected  size  of  the  source.   This  point  is  illustrated 
if  we  imagine  trying  to  produce  incoherent  light  in  the  entrance  slit  of  a  monochromator . 
This  can  be  done  by  replacing  the  slit  by  an  incandescent  filament  or  by  illuminating  the 
slit  by  an  infinite  incoherent  source  (an  integrating  sphere  with  one  port  being  the  entrance 
slit  would  probably  do) .   In  neither  case  is  there  any  defined  direction  of  the  sort  exhib- 
ited by  the  laser. 

So  far  we  have  limited  our  examples  to  the  two  extreme  cases  of  complete  incoherence  and 
complete  coherence.   The  reason  for  this  is  that  these  are  the  two  cases  for  which  the 
integrations  are  most  easily  carried  out  —  incoherence,  because  of  the  delta  function  —  and 
coherence,  because  of  the  factorization  into  the  product  of  a  simpler  integral  and  its 
complex  conjugate.   The  theory,  however,  is  perfectly  general  and  applies  to  any  expression 
for  the  cross-spectral  density  with  any  degree  of  coherence. 

As  a  final  numerical  example  of  the  application  of  the  propagation  of  coherence  to 
practical  radiometry  we  will  calculate  for  partially  coherent  radiation  the  slit-scattering 
function  of  the  monochromator  of  the  last  example  (see  also  reference  [10.18]).   We  will 
assume  that  the  entrance  slit  is  irradiated  by  a  monochromatic  incoherent  source  of  width   2aQ 
located  on  the  optic  axis  at  a  distance   R    in  front  of  the  slit  and  will  calculate  the 
distribution  of  spectral  fluence  rate  in  the  exit  plane.   As  we  have  seen  [eq.  (10.52)]  the  x- 
dependent  part  of  the  degree  of  coherence  of  the  light  at  the  plane  of  the  entrance  slit 
depends  upon  the  ratio   a  /R   and  can  be  varied  from  total  coherence  to  total  incoherence  by 
changing  this  ratio  from   0   to   °°.    Based  upon  eq.  (10.47),  we  see  that  we  can  write  the  x- 
dependent  factor  of  the  cross-spectral  density  function  in  the  plane  of  the  entrance  slit  as 


W(x)(x.;x2) 


2a  A 
o  o 


2a  (x, -x_) 
o   1   2 

AR 


(10.122) 


Now,  using  eq.  (10.105),  we  obtain  for  the  cross  spectral  density  in  the  plane  of  the  exit 
slit: 


w^(x) 

(xj 

X") 

= 

8a' 

R 
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2a  A 
o  o 

RR"A2 

. 

sine 

2a' 

A 

2a^ (x,-x_) 
o   1   2 


AR 


2a 
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'   /*2_   *2_\ 
"  \R     R"/ 


(10.1231 


■dx, 'dx_ 


In  this  expression   2a'   is  the  width  of  the  instrument  aperture  —  usually  this  is  the  width 
of  the  grating  --  and   R   and   R"   are  the  distances  from  the  entrance  and  exit  slits  to  the 
collimating  optics.   If  we  confine  our  attention  to  the  spectral  fluence  rate  in  the  exit 
plane  we  can  set  x'l    =  x"  =  x"   to  give  the  x-component  as 


p..  (x) 

Ft,A 

(x") 
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8a'  a  A 
o  o 

R  RR"A2 
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o 
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(10.124) 
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■dx, -dx_  , 


The  results  of  performing  the  indicated  integrations  numerically  for  the  monochromator 
of  the  previous  example   (R  =  R"  =  0.8  [m] ,   a'  =  0.0375  [m] ,   a  =  25  [ym] ,   A  =  514  [nm] ) 
are  shown  in  figure  10.9.   The  two  solid  curves  are  the  extremes  of  complete  coherence  and 
complete  incoherence;  the  incoherent  curve  possesses  the  highest  slit-scattering  function  in 
the  far  wings  and  is  quite  flat  and  structureless  in  the  middle.   The  dashed  curves  show 
cases  of  partial  coherence  and  lie  generally  between  the  two  extremes.   The  curves  in  this 
figure  have  been  area  normalized  in  order  to  emphasize  the  changing  shape  of  the  slit-scat- 
tering function.   Again  the  figure  shows  the  distribution  of  fluence  rate  in  the  exit  plane 
of  the  monochromator.   This  would  be  observable  with  a  very  narrow  exit  slit.   For  a  normal 
exit  slit  the  finer  structure  would  be  largely  smoothed  out  and  only  the  average  far-wing 
behavior  would  distinguish  the  cases  of  differing  coherence. 

Clearly,  if  the  slit-scattering  function  of  a  spectroradiometer  must  be  measured  accu- 
rately for  later  use  (deconvolution)  with  an  unknown  source  it  is  necessary  that  the  state  of 
coherence  of  the  calibration  beam  be  similar  to  that  of  the  unknown  to  be  measured.   For  many 
situations  this  requirement  may  not  be  as  restrictive  as  it  appears.   In  figure  10.9,  for 
example,  the  transition  from  slit-scattering  functions  which  are  barely  distinguishable  from 
completely  coherent  to  those  which  are  barely  distinguishable  from  completely  incoherent 
occurs  as   a  /R   changes  from  about  0.0  3  to  0.06.   Thus  many  radiation  fields  could  probably 
be  treated  as  one  extreme  or  the  other.   In  addition,  by  changing  some  parameter  of  the 
measurement  environment  such  as   R  ,   the  source-radiometer  distance,  it  may  be  possible  to 
shift  the  slit-scattering  functions  for  both  sources  to  the  same  side  of  the  transition 
region,  so  that  for  that  particular  pair  of  measurements,  at  least,  a  common  slit-scattering 
function  will  be  satisfactory. 

More  generally,  we  can  say  that  for  an  imaging  system,  if  the  coherence  interval  in  the 
object  plane  is  very  much  greater  than  the  width  of  the  point  spread  function  of  the  optical 
system  in  the  object  plane,  then  the  image  will  resemble  the  coherent  limit.   If  the  coher- 
ence interval  is  very  much  smaller  than  the  width  of  the  point  spread  function  then  the  image 
will  resemble  the  incoherent  limit.   The  point  spread  function  of  an  imaging  system  is  the 
distribution  of  fluence  rate  in  the  image  of  a  point  source.   Ideally  a  point  source  should 
be  imaged  into  a  point  image  but  due  to  diffraction,  aberrations,  and  focusing  errors  it 
never  will  be.   Here  we  are  only  interested  in  the  scattering  by  diffraction.   For  perfect 
imagery  and  rectangular  apertures  we  can  calculate  the  point  spread  function  from  eq.  (10.105) 
Imagine  the  source  shrinking  to  a  very  small  area  around  x,  =  x_  =  0 .   The  source  cross- 

(x)  (x) 

spectral  density   W,   (x, ;x_)   will  then  approach  a  constant,   F'  ' (0) .   Thus,  the  x-dependent 

A      _L    Z  L. ,  A 

factor  of  the  image  cross-spectral  density  function  will  be  proportional  to: 

sinc(2a'x£/AR") • sine (2a 'x^/AR" )  (10.125) 

2 

and  the  fluence  rate  in  the  image  will  be  proportional  to  sine  (2a'x"/AR").   There  is,  of 

course,  a  similar  y-dependent  factor  and  the  complete  point  spread  function  for  a  rectangular 
aperture  contains  the  product  of  these  two  factors.   In  one  dimension  the  point  spread 
function  is  just  the  slit-scattering  function  for  an  infinitesimal  slit.   To  obtain  the  point 

spread  function  in  the  object  plane  we  interchange  the  roles  of  image  and  object  and  write 

2 

the  point  spread  function  as  proportional  to   sine  (2a'x/AR).   The  width  of  this  function  is 

then,  say,  twice  the  value  of  x  at  which  the  function  is  reduced  to  half  its  peak  value. 

Actually,  with  sine  functions  it  is  more  common  to  use  the  width  between  zeros  of  the  function, 

but  to  be  consistent  with  our  choice  for  coherence  interval  [eq.  (10.56)]  we  will  continue  to 

2 
use  the  full  width  at  half  height.   Since   sine  (u)  =  0.5   for   u  =  0.44,   we  find  that  the 

width  of  the  point  spread  function  is  about   0.4AR/a'.    This  can  also  be  interpreted  as  the 

resolution  limit  of  the  imaging  system;  two  objects  lying  closer  together  than  this  will  not 
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Figure  10.9.    Relative  distribution  in  the  x-direction  of  fluence  rate 

in  the  image  of  a   50  [ym]-wide,  uniform,  partially  coherent 
source  at  normal  incidence  for  various  degrees  of  partial 
coherence.   Plots  of  the  double  integral  of  eq.  (10.124) 
for  the  same  parameter  values  as  in  figure  10.8  and  area 
normalized,  in  the  interval  shown,  to  the  area  under  the 
6   =0   curve  in  figure  10.8.   a^/R^  =  0.00  (coherent); 


a  /R  =  0.04  ; 

o   o 


—  a  /R 
o  o 


—  a  /R 
o   o 

0,045;   • 


- —  a  /R   =0.05 
o   o 


a  /R   =«>  (incoherent) 
o   o 
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be  resolved  in  the  image.   For  our  present  purpose  of  considering  the  effects  of  coherence  in 
radiometry,  the  significance  of  the  point  spread  function  width  is  that  if,  for  two  different 
sources,  the  characteristic  coherence  intervals  at  the  entrance  slit  of  a  monochromator 
radiometer  are  either  both  very  much  greater  than   AR/a'   or  both  very  much  smaller  than 
AR/a'   then  the  monochromator  slit-scattering  function  will  be  similar  for  both  sources.   If 
one  or  both  coherence  intervals  are  comparable  to   AR/a1   or,  worse,  if  one  is  large  and  the 
other  small,  then  the  slit-scattering  functions  may  differ  considerably  for  the  two  sources, 
and,  depending  upon  their  spectral  distributions,  this  may  significantly  affect  the  accuracy 
of  any  direct  comparison  of  the  two  sources  with  this  radiometer.   To  test  if  this  is  a 
problem  or  to  ameliorate  the  problem  if  it  exists  one  can  usually  increase  the  point  spread 
function  width  by  decreasing  the  imaging  aperture  half width   a'   or,  as  mentioned  above,  one 
may  be  able  to  change  the  coherence  interval  by  changing  the  distance  from  the  source  to  the 
slit. 

Any  time  high  resolution  is  sought  in  any  dimension  of  the  measurement  situation  — 
spectral,  temporal,  angular,  or  spatial  —  the  coherence  state  of  the  radiation  can  affect 
the  measurement  results.   In  the  present  example  spectral  resolution  is  achieved  by  means  of 
spatial  and  angular  resolution  of  an  image  within  the  monochromator.   This  image  resolution 
is  limited  by  the  distortion  caused  by  diffraction,  and  the  magnitude  and  nature  of  the 
distortion  due  to  diffraction  are  dependent  upon  the  coherence  state  of  the  radiation  entering 
the  radiometer.   As  a  further  example  let  us  consider  the  image  formed  by  a  lens  of  a  slit 
upon  which  a  laser  speckle  pattern  falls.   If  we  assume  that  the  speckles  consist  of  homo- 
geneous regions  all  of  the  same  size  but  with  random  electromagnetic  field  amplitudes  and 
phases  then  the  fluence  rate  in  the  image  may  be  distributed  somewhat  like  one  of  the  curves 
in  figure  10.10.   The  curves  in  this  figure  were  computed  for  the  same  monochromator  as  in 
the  previous  examples,  but  with  the  slit  coherently  irradiated  with  a  cross-spectral  density 
given  by: 

WA (x) (x1;x2)  =  A-cosfaCx,) ] «cos[a(x2) ] •exp[ig(x1)-i0(x2) ]  (10.126) 

where   a(x)   and   6  (x)   are  angles  between   -it   and   it   which  are  constant  within  each 
interval  of   1  [yra]   but  random  from  interval  to  interval  across  the  width  of  the  slit.   The 
two  curves  are  for  two  different  samples  of  random   a(x)   and   3 (x)   values.   Now  suppose 
that  the  speckle  pattern  moves  with  time  (e.g.,  is  generated  by  a  rotating  ground  glass)  or, 
equivalently ,  suppose  that  the  primary  source  is  polychromatic  in  such  a  way  that  the  phase 
in  each  speckle  changes  in  time  randomly  with  respect  to  its  neighbors.   Then  the  distribu- 
tion of  fluence  rate  in  the  image  will  evolve  continuously  from  curve  to  curve  among  the 
family  of  curves  of  which  those  of  figure  10.10  are  representative.   If  a  small  detector  in 
the  image  plane  possesses  sufficient  temporal  resolution  and  if  the  speckle  pattern  evolves 
sufficiently  slowly  then  the  detector  will  follow  this  evolution  and  an  erratic  signal  output 
will  result.   On  the  other  hand,  if  the  detector  in  the  image  plane  has  poor  temporal  resolu- 
tion and  if  the  speckle  pattern  changes  rapidly  then  a  steady  average  signal  output  will 
result  and,  in  fact,  the  average  distribution  of  fluence  rate  in  the  image  will  be  exactly 
that  of  an  incoherently  irradiated  slit:   the  average  of  a  large  number  of  random  speckle 
pattern  image  distributions  such  as  those  of  figure  10.10  approaches  the  shape  of  the  inco- 
herent curve  of  figure  10.9.   Thus,  a  rapidly  changing,  fine,  speckle  pattern  is  a  model  for 
incoherent  light  and  we  see  that  with  high  temporal  resolution  the  results  of  a  single 
radiometric  measurement  can  depend  upon  the  coherence  state  of  the  radiation.   Although  in 
this  example  we  have  employed  the  time  averaging  of  an  evolving  speckle  pattern  to  simulate 
'natural'  incoherent  light,  the  average  in  the  y-direction  —  along  the  height  of  the  slit  — 
of  a  stationary  speckle  pattern  can  produce  the  same  effect.   We  have  assumed,  in  the 
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Figure    10.10. 


Two  samples  of  the  relative  distribution  in  the  x-direction 
of  fluence  rate  in  the  image  of  a   50  [ym]-wide  coherent 
source  consisting  of   1  [urn]   cells  with  amplitude  and  phase 
constant  within  each  cell  but  random  from  cell  to  cell. 
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calculations  of  figure  10.10,  that  the  height  of  the  slit  spans  one  speckle  width  or  less  and 
have  ignored  the  y-dimension  of  the  problem.   If  the  slit  and  detector  cover  many  speckles  in 
the  vertical  direction  then  again  the  detector  will  integrate  over  many  random  speckle 
patterns  in  this  direction  and  the  observed  fluence  rate  distribution  in  the  image  plane  will 
agree  with  that  of  an  incoherently  irradiated  slit  even  when  the  speckle  pattern  is  stationary 
and  the  irradiation  is  monochromatic. 

MORE  on  IMAGING  —  PUPIL  FUNCTIONS  and  ABERRATIONS.    In  deriving  eq .  (10.98)  for  the  cross- 
spectral  density  function  in  the  image  of  a  source  we  made  two  applications  of  eq.  (10.36)  -- 
the  law  of  propagation  of  cross-spectral  density.   This  gives,  among  other  factors  in  the 
integrand,  a  phase  factor 


[nr  '(si  +  si"  -  s2  -  V)] 


exp  |_ — 5^ —  •(s1  +  s1"  -  s2  -  s2")J  .  (10.127) 

The  distance   s,  +  s, "   is  the  total  optical  path  length  from  point   (x, ,y, )   in  the  source 
plane  to  point   (xV  ,yV)   in  the  image  plane  via  a  point   (x,1  ,y,' )   in  the  intermediate  aperture 
plane.   Similarly   s„  +  s'l      is  the  total  optical  path  length  between  a  second  pair  of  source 
and  image  points  via  a  second  point  in  the  aperture.   In  that  derivation  we  made  some  approx- 
imations for  this  phase  factor  and  then,  in  passing  to  eq.  (10.99)  ,  we  integrated  over  all 
pairs  of  points  in  the  aperture.   It  is  not  necessary,  of  course,  that  there  be  any  focusing 
optic  in  the  aperture  plane  or  that  there  be  a  recognizable  image;  the  "source",  "aperture", 
and  "image"  planes  may  be  any  three  consecutive  cross-sections  of  the  beam  and  the  sequential 
application  of  eq.  (10.36)  is  still  valid.   The  particular  simplifications  which  led  to  eq. 
(10.98)  will  no  longer  apply  in  the  absence  of  perfect  imagery.   Instead,  the  second  four- 
fold integral  in  eq.  (10.98)  will  be  replaced  by 

JJff   exp[^  •(s1  +  s^  -  s2  -  s^j-dx'-dy'-dx'-dy'  (10.128) 

aperture 

where  the  integrations  cover  the  aperture  twice.   Starting  from  expression  (10.12  8)  we  now 
want  to  discuss  a  generalization  of  eq.  (10.98)  which  can  include  arbitrary  apertures  and 
departures  from  perfect  imagery. 

First  of  all,  instead  of  limiting  the  integrations  to  the  aperture  we  can  introduce  a 
pupil  function1   P(x',y')   whose  value  is   1   inside  the  aperture  and   0   outside  the  aper- 
ture and  then  let  the  integrations  cover  the  entire  aperture  plane: 

CO       CO 

ff   if  P(xi'Y;P-P*(x2,y2)-exp  [^  *(si  +  sl  "  s2  "  s^)] -dx^dy^dx^ -dy2 .        (10.129) 

—  CO     —CO 

For  our  rectangular  aperture,  for  example,  we  would  have 

P(x',y')  =  1     for   -a  <  x'  <  a    and    -b  <  y'  <  b  (10.130) 

=   0    otherwise. 


Strictly  speaking  the  pupil  function   P(x',y')   will  also  depend  upon  the  source  point 
(x,y)   and  the  destination  point   (x",y")   of  the  radiation  passing  through  the  aperture. 
This  dependence  is  generally  neglected  and  we  don't  explicitly  show  it  in  our  notation. 
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So  far  this  is  exactly  equivalent  to  expression  (10.128).   However,  we  can  generalize  this 
expression  by  allowing   P(x',y')   to  take  on  values  other  than   0   and   1.    Intermediate 
values  of   P(x',y')   can  describe  optical  systems  involving  partially  absorbing  media  such  as 
the  "aperture"  formed  by  the  slide  in  a  slide  projector.   If  the  aperture  introduces  phase 
variations,  such  as  a  slide  made  of  wavy  glass,  then  this  can  either  be  put  into  the  pupil 
function  by  making   P(x',y')   complex  or  be  accounted  for  in  the  calculation  of   s  +  s". 

The  four- fold  integral  in  expression  (10.129)  can  be  seen  to  be  the  product  of  two 
similar  integrals: 

oo 

If   P(x|,y|).exp[^i  •  (s1  +  s£)]  -dx^-dy^ 

—  OO 

oo 

•  jj   P*(x^,y^) -expT-  ^—   .(s2  +  s^)  •dx^-dy^l  (10.131) 

—  00 

=   K(x1,y1;x^,yp  •  K*  (x2  ,  y2  ;  x^  ,yjj.)  , 
where 


K(x,y;x",y")   =  ||   p  (x '  ,  y '  )  •  exp  [^i  •  (s  +  s"  -  R  -  R")]-dx'-dy'         (10.132) 

and   s  +  s"   is  the  optical  path  length  between   (x,y)   and   (x",y")   via   (x',y').    R   is  the 
perpendicular  optical  path  length  from  the  source  plane  to  the  aperture  plane  and   R"   that  from 
the  aperture  plane  to  the  image  plane.   The  exponential  factors  involving   R   and   R"   which  we 
have  introduced  in  eq.  (10.132)  cancel  out  in  forming  eq.  (10.131).   The  function   K(x,y;x",y") 
is  the  transmission  function  [10.2],  [10.3],  and  [10.5]  of  the  system.   The  transmission  func- 
tion is  similar  to  those  functions  variously  known  as  point-spread  function,  impulse  function, 
Green's  function,  etc.:   if  we  think  of   K(x,y;x"y")   as  a  function  of  the  position  of  the 
point   (x",y")   in  the  image  plane,  then  it  describes  the  field  in  the  image  plane  caused  by 
a  unit  point  source  of  fixed  phase  at  the  source  point   (x,y) .   In  terms  of  the  transmission 
function  we  can  now  write  the  cross-spectral  density  function  in  the  image  plane  as 

W^x^y-jx^y^)   =  r2r„2x4  '  fjff  WX  (xl'  yl;X2  >*2> 

source 

(10.133) 

•  Kx  (x1,y1;x^,y^)  -K*  (x2,y2;x2' ,  y£ )  •dx1-dy1-dx2  -dy2. 

This  is  the  generalization  of  eq.  (10.98)  which  applies  whatever  the  aperture  and  whether  or 
not  imagery  is  good.   Of  course,  it  requires  being  able  to  express  the  optical  path  lengths 
s  +  s"   in  terms  of  the  coordinates  of  the  source,  aperture,  and  image  points. 

When  the  aberrations  or  focusing  errors  are  small  we  can  write 

s  +  s»    =    R  -  Ml  +  yy'  +  R..  _  *'*"  ;,y'y"    +A  uo.i34) 

where  the  first  four  terms  on  the  right  come  from  the  approximations  which  lead  to  eq. 
(10.98)  [see  also  eq.  (10.39)]  and   A   is  a  correction  term.   Notice  that  if   A   is  zero 
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there  exists  a  point  in  the  image  plane  at   x"  =  -x«R"/R/   y"  =  -yR"/R   for  which   s  +  s" 
is  independent  of  the  aperture  point   (x',y') .   This  means  that  for  each  source  point   (x,y) 
there  is  a  conjugate  image  point  for  which  all  rays  through  the  aperture  travel  the  same  optical 
path  length.   Thus,  eq.  (10.134)  with   A  =  0,   and  consequently  also  eqs.  (10.98)  and  (10.99) 
correspond  to  perfect,  aberration-free  focusing.1 

Equation  (10.133)  with  the  transmission  function  defined  by  eq.  (10.132)  is  always  valid 
provided  only  that  the  optical  path  lengths   s,   and   s„   (and   s"   and   s")   which  appear  in 
the  denominator  of  the  propagation  law  [eq.  (10.36)]  are  sufficiently  constant  to  be  approxi- 
mated by   R   (and   R" ) . 2   If  the  departures  from  perfect  imagery  are  small  the  approximation 
of  eq.  (10.134)  may  be  useful.   As  an  example  of  this  approximation  suppose  that  although  the 
true  focal  plane  is  located  at  the  distance   R"   from  the  aperture  the  "image"  plane  of 
interest  is  a  parallel  plane  a  small  distance  t,      beyond  the  focal  plane.   Then  a  first 
approximation  to   A   is  given  [10.2]  by 

A   =   c-fl  -  — V1—)-  (10.137) 

\       2R"Z    / 

Focus  errors  are  treated  in  full  in  reference  [10.2]  ,  Chapter  8,  and  aberrations  are  treated 
in  Chapter  9  of  the  same  reference.   It  will  rarely  be  possible  to  evaluate  in  terms  of 
convenient,  known  functions  the  integrals  which  result  from  the  inclusion  of  focus  errors, 
aberrations  or  exotic  pupil  functions.   Usually,  the  numerical  evaluation  of  the  transmission 
functions  from  eq.  (10.132)  with   s..  ,  s",  s_   and   s"   calculated  directly,  will  be  prefer- 
able to  trying  to  identify  and  work  with  the  terms  of  the  expansion  of  eq.  (10.134).   By 
writing  complex  exponentials  as  the  sum  of  a  cosine  term  and   i   times  a  sine  term,  the 
integrals  of  eqs.  (10.132)  and  (10.133)  can  ultimately  be  expressed  as  a  sum  of  two  integrals  — 
one  real  and  one  imaginary  —  which  can  be  computed  independently.   Eventually,  in  the  chain 
of  calculations  when  an  observable  such  as  fluence  rate  or  radiant  intensity  is  computed, 
these  components  will  combine  to  produce  a  single,  real,  result. 


1  In  this  case  [or  if  we  put  any  phase  error  factor,   exp  — r — ,   into  the  pupil  function 
P(x',y')]   we  can  write  the  transmission  function  as 

00 

K(x,y;x",y")   =JJ   p  (x'  ,  y' )  -exp[-  ffi  -(XX'  {j  YY'  +  X'X"  ^  'y" )]  -dx'  -dy  ■       (10.135) 


Now  by  letting  E,   =  x/R  +  x"/R"   and   n  =  y/R  +  y"/R"   this  integral  can  be  cast  in  the  form 
of  a  two-dimensional  Fourier  transform  of  the  pupil  function  (as  modified  by  including  any 
phase  factors  for  aberrations  or  focus  error) : 

K(x,y;x",y")   =   K(£,n)=ff   P  (x'  ,  y '  )  -exp  [-  2*i  .  (x>  g  +  y  •  n)j  -dx'  -dy  '  .  (10.136) 

—  00 

Many  of  the  equations  of  this  chapter  can  be  expressed  somewhat  more  compactly  by  employing 
Fourier  transform  shorthand.   However,  we  make  no  use  of  this  because  we  suspect  that  many 
readers  never  really  feel  comfortable  with  Fourier  transforms  and  we  fear  that  their  use 
sometimes  tends  to  obscure  the  simple  physics  and  the  essential  approximations  involved. 


2A11  lengths,  including  R   (and   R")»   must  of  course  be   >>  A, 
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The  CLASSICAL  RADIOMETRIC  QUANTITIES.    So  far  in  this  chapter  we  have  avoided  saying  much 
about  spectral  radiance  and  this,  no  doubt,  seems  curious  since  spectral  radiance  lies  at  the 
heart  of  all  the  rest  of  the  chapters  of  this  Manual.   The  reason  is  that  the  classical 
concept  of  spectral  radiance  is  inconsistent  with  a  wave  theory  of  light  because  the  very 
limiting  process  required  to  define  a  ray  and,  hence,  radiance  (letting  the  areas  of  two 
apertures  go  to  zero  --  see  Chapter  2  [10.19])  guarantees  such  severe  diffraction  in  the 
presence  of  coherence  that  the  direction  the  "ray"  came  from  is  totally  indeterminate.   In 
consequence,  no  function  of  the  cross-spectral  density  can  possess  all    the  properties  of  the 
classically  defined  spectral  radiance  for  all  fields  [10.14].   Functions  which  possess  many   of 
the  properties  of  spectral  radiance  can  be  defined  and  that  given  in  eq.  (10.20)  is  commonly 
used. 

Equation  (10.20)  can  be  inferred  from  eq.  (10.36)  for  the  propagation  of  the  cross-spectral 
density  function  using  the  classical  relationship  between  spectral  radiance  and  spectral 
fluence  rate.   Consider  a  source  plane   z  =  0   where  the  spectral  radiance  is  desired.   In 
this  plane  we  assume  that  the  cross-spectral  density   W,  (x.  ,y. ;x_ ,y_)   is  known.   To  measure 
the  spectral  radiance  we  would  place  a  spectroradiometer  at  some  distance   R   away  from  the 
source,  pointing  from  the  direction   0 ,  <j>   to  the  point   x,y  where  the  spectral  radiance  is 
desired.   At  the  entrance  aperture  of  the  radiometer  the  cross-spectral  density  is  given  by 
eq.  (10.36).   Let  us  express  eq.  (10.36)  in  terms  of  the  average  and  difference  coordinates 
introduced  in  eqs.  (10.77)  and  (10.82): 

x  =  h(x1   +  x2)      y  =  h(Y1   +  y2) 

(10.138) 

?  =  xi  "  x2      n  =  y-L  -  y2. 

Then   x,  =  x  +  hKi      y-,  =  y  +  hr\,      x?  =  x  -  %E,i      y2  =  Y   ~    ^n ,   and  the  product  of  the  area 
elements  becomes1   dx*dyd£*dn.    We  also  make  similar  substitutions  for  the  primed  coordi- 
nates in  the   X'  -  Y'   plane  and  define   R   as  the  distance  between  the  source  point   P(x,y,0) 
and  the  distant  point   P' (x',y',z')   where  the  radiometer  is  located: 

R2   =   (x'  -  x)2  +  (y'  -  y)2  +  z'2.  (10.139) 

In  terms  of  these  variables  the  distances   s,   and   s_   of  eq.  (10.36)  become 


!1  =  Vr2  +  (x'-x)-(?'-0  +  (y'-y)  •  (n'-n)  +  %[(£'-U2  +  (n'-n)2] 


S2 


=  ^R2  -  (x'-x)-(C'-U  -  (y'-y)  •  (n'-n)  +  ^[(€'-?)2  +  (n'-n)2] 


(10.140) 


and  we  can  write  eq.  (10.36)  as 


w'(x'+i-2?',  y'+hr]';   x'-kV .   y'-W)  =      wx(x+^,  y+%n;  x-ht,   y-%n) 


(10.141) 


2 

-xPpir  •(S1-S2)]  '    2Z'2,2  'dx-dydC-dn. 
sl  s2  A 


'See  footnote  on  page  27. 
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Now  we  make  the  usual  approximations;  we  expand  the  square  roots  in  eqs.  (10.140)  in  powers  of 
1/R   as  in  going  from  eq.  (10.38)  to  eq.  (10.41) ,  keeping  only  the  lowest  order  terms  to  give 


wMx'+J^',   y'+'sn';   x'-HZ' ,   y'-^n')      =  w, (x+%5,   y+Jgn;   x-ht,   y-'sn) 


-fill 


X 
source 


(10.142) 


.2 


exp  |||i  •  [(x'-x)«  (£'-€)   +    (y'-y)- (n'-n)]  |    •   -f^j  -dx-dydC-dn 


The  spectral  fluence  rate  at  P' (x',y',zl)   is  obtained  by  setting   C'=n'=0: 

F^x(x',y')  =  jjjj  wx(x+^,  y+hx]',   x-kK.   y-^n) 

source  (10.143) 

•  exp<-  j^-   •[C'(x'-x)  +  n*(y'-y)]|'  Z4  2  •dx'dyd^dn. 

The  contribution  to  the  fluence  rate  from  the  element  of  solid  angle   dco  =  cos,   *dx'dy   is  just 

R^ 


dF^A(x,,y',e,(i))  =  II  wA(x+%?/  y+%n;  x-%£,  y-^n) 


source 


(10.144) 


:xp|-  j^-   •[^•(x'-x)  +  n*(y'-y)]>*  ^— j  •d^-dn-do) 


Therefore,  from  the  classical  relation  {eq.  (4.19)  [10.9] }1 


dFMx'^Se,^) 
L*  (x',y',e,<j>)   =   — - 3-r (10.146) 


dw 
we  obtain 


*If  the  expression  for  spectral  radiance  were  derived  from  eq.  (10.141)  without  making  the 
large-distance  approximation  of  eq.  (10.142)  the  formula  would  depend  upon  R.    Since  this  is 
incompatible  with  classical  radiometry  it  is,  strictly  speaking,  necessary  to  define  radiance 
as  the  far  field  limit  of  the  derivative  of  the  fluence  rate: 

Lx(x,y,e,*)   =  IZ^     F'/X(x',y)  (10'145> 

where  R,   as  given  by  eq.  (10.139),  is  the  distance  from  the  source  point   (x,y)   to  the 
detector  point   (x',y'),   and  F'  , (x',y')  =  W, (x* ,y' ;x' ,y' )   as  usual.   In  this  way  L,   will 
depend  only  upon  the  direction  (&,$)      from  the  point   (x,y) .    In  our  derivation  we  have 
effectively  taken  this  limit  when  we  make  the  approximation  leading  to  eq.  (10.142). 
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•  / v  I    t>l  " 


L^(x',y'  ,6,<j)) 


r 


w,  (x+>sC,  y+H>n;  x-4Cr  y-%n) 


rx 

source 

(10.147) 


exp{-  ^  -O(x'-x)  +  n-  (Y'-Y)]}-dC-dn. 


In  this  equation,   6    and   <J>   specify   the  direction  from  the  source  point   (x,y)   to  the 
detector  point   (x',y'): 

cos9  =  z'/R 
sine«cos<J>  =  (x'-x)/R  (10.148) 

sin6«sin<{>  =  (y'-y)/R. 

If  we  use  these  relations  to  eliminate  the  primed  coordinates  and  remember  that  geometrically 
radiance  is  invariant  along  a  ray  we  obtain  eq.  (10.20) 


LA(x,y,e,<t>)  =  £2f^  •  I  wx(x+^c,  y+%n;  x-^e,  y-%n) 


A 

source 


(10.149) 

2-rri    ..._  ..~l._. 

exp 


>  -  — r—  •  sine  •  (C 'COs4>  +  n*sin<|>)  'dC'dn. 


This  is  the  spectral  radiance  at  the  source  point   (x,y)   in  the  direction   (6,$).    The 
integration  covers  the  entire  source  plane,  or  beam  cross-section  if  the  "source"  plane  is  not 
a  primary  source.   We  see  that  the  cross-spectral  density  over  the  entire  source  plane  con- 
tributes to  the  spectral  radiance  at  any  one  point  of  the  source.   This  is  consistent  with  the 
idea  that,  by  diffraction,  all  points  in  the  source  can  contribute  something  to  the  signal 
received  by  a  radiometer  nominally  aimed  at  a  single  point. 

The  spectral  irradiance  at  the  point  (x*  ,y' ,z' ) ,       irradiated  by  a  source  centered  at  the 
origin  can  be  obtained  similarly  from  the  relation  {eq.  (4.20)  [10.9]}  between  irradiance  and 
fluence  rate: 

dE'  (x'  ,y'  ,6,<}>)   =   cos6r-dF^(x'  ,y'  ,&,<$>)  .  (10.150) 

Using  eq.  (10.144)  for  the  element  of  spectral  fluence  rate  we  obtain 


(10.151) 


EA(x\y,,z')  =   J]  [J  wA(x+JsC,  y+^n;  x-hZ,   y-^n) 

2 
•  exp|-  ^li  -O(x'-x)  +  ir(y'-y)]}-  ^-j  -cosG^dC-dn-dx-dy. 

R  A 

which  differs  from  eq.  (10.143)  for  spectral  fluence  rate  only  by  the  extra  factor   cos6  . 
The  angle   9    is  the  angle  between  the  normal  to  the  irradiated  surface  at   (x',y')   and  the 
direction  to  the  source  point   (x,y).    If  the  irradiated  plane  is  parallel  to  the  source 
plane  then   cosG   =  cos6  =  z'/R.    In  both  eqs.  (10.143)  and  (10.151)   R   is  a  function  of   x 
and   y   [eq.  (10.139)]  which  makes  the  integrals  difficult.   If  we  can  assume  that  the  source 

dimensions  are  small  compared  with  Var"  —  the  Fraunhofer  approximation  —  then 

2    2    2  h 

R  «  (x1  +y'  +z '  )    and  these  equations  become: 
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Ffc  A(x' ,y' ,z') 


cos 


2  2 


w,  (x+J^,  y+^n;  x-^,  y-Jgn) 


source 


[10.152) 


•  exp 


r  2zi 

L   XR 


•»] 


(C-x'  +  n-y' )  -dx-dyd^-dn 


and 


Ex(x',y,,z')   =   cosOr-Ft  A(x' ,y',z' 


(10.153) 


where   cos0  =  z '  /R   and  both  the   x,y   and   £ ,  n   integrations  independently  cover  the  source 
plane. 

The  rest  of  the  common  radiometric  quantities  can  now  be  easily  expressed  in  terms  of  the 
cross-spectral  density.   From  eqs.  (10.16),  (10.143),  and  (10.148)  we  find 


ix(e,4>) 


cos  6 


wx(x+Js£,  y+hT]-,   x-JjC,  y-%n: 


source 


(10.154! 


exp  -  — r —  *sin0*  (5*cos(j)  +   n*sin<f>)|  •dx,dy,d5*dn 


for  the  spectral  radiant  intensity  of  the  source  in  the  direction   (6,4>).    The  spectral 
radiant  exitance  from  the  source  plane  is  obtained  by  integrating  the  spectral  radiance  [eq. 
(10.149)]  over  the  hemisphere  with  respect  to  projected  solid  angle   cosO'doj  =  cosG  •  sin0  *d6  *d<J> , 
The  integrations  over   6   and   <))   can  be  carried  out  to  give 


M^^y)      =      ||   wA(x+J-25,   y+^n;   x-^,   y-%n) 


2jn_ 

,2 


/2tt        j'    2,     2\  /2tt       Jr2^    2\ 

sinir  •  V^  +v      cos{—  '  V?  +r)  ) 


(r-v 


r2^     2V 

K  +n  / 


(?  •  v^2: 


(10.155) 


•dC-dn. 


For  highly  incoherent  fields,  these  expressions  for  the  radiometric  quantities  possess 
all  the  properties  of  the  corresponding  classically  defined  radiometric  quantities.   However, 
for  fields  possessing  significant  coherence  they  possess  all  the  right  properties  except  that 
the  spectral  radiance  and  spectral  radiant  exitance  in  some  circumstances  can  be  negative  and 
can  be  non-zero  in  regions  of  a  source  plane  outside  of  the  physically  defined  source  area 
itself.   As  we've  indicated  these  quirks  arise  because  these  are  classical  geometrical  optics 
concepts  which  are  not  strictly  definable  in  the  real  world  of  wave  optics.   Fortunately  it 
appears  that  these  deficiencies  aren't  likely  to  arise  in  most  cases  of  interest  [10.14]  so 
that  eq.  (10.149)  and  eqs.  (10.152)  to  (10.155)  can  be  accepted  as  the  appropriate  extension 
of  classical  radiometry  to  include  coherence. 

We  illustrate  the  calculation  of  the  classical  radiometric  quantities  from  the  cross- 
spectral  density  by  an  example  of  a  source  whose  cross-spectral  density  function  is  given  by 


Wx(x1,y1;x2,y2)   = 


A2     /   x2+y2\   .  (2 

A  «exp  -  ^-  -sine  T 

o    Fl    2a2    I  U 


(M 


r2x  2 


(10.156) 
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where   x,   y,  E,      and   n   are  the  average  and  difference  coordinates  defined  in  eqs.  (10.138). 
The  spectral  fluence  rate  in  the  source  plane  is  a  Gaussian  function  of  the  distance  from  the 
origin: 


Ft/A(x,y) 


A^-exp  (-r2/2o2) 


(10.157) 


2     2     2 
where   r   =  x   +  y  .    The  degree  of  coherence  in  the  source  plane  is  a  sine  function  of  the 

distance  between  the  two  points  whose  coherence  is  being  considered: 


(x1,y1;x2,y2)   =  sine  (  —  • p) 


(10.158) 


2     2     2  2  2 

where   p   =  £   +  n   =  (x1-x0)   +  (y-,-y9) 


determines  the 


1  ~2'       '  ^l"1?'  "    We  see  t'iat  the   parameter 
coherence  interval  in  the  source  and  the  parameter   a   is  a  measure  of  the  size  of  the  source. 
Since  this  is  clearly  a  quasihomogeneous  source,  because  the  cross-spectral  density  is  a  prod- 
uct of  a  function  of  the  average  coordinates,   x   and   y,   and  a  function  of  the  difference 
coordinates,   £   and   n,   (this  implies  that  a    >>  A/g)  ,   let  us  first  simplify  the  expres- 
sions for  the  classical  radiometric  quantities  by  applying  them  to  a  quasihomogeneous  source. 
Upon  substituting   F.  , (x,y) *g(£»nj   for  the  cross-spectral  density  in  eqs.  (10.149),  (10.152), 

t  ,  A 

and  (10.155)  we  obtain: 


Lx  (x,y,6,(J)) 


cos8 


2   't,i 


(x,y) 


g(£,n) 


;i0.159) 


•  exp  -  — r — •  sine  •  (C"Cos<t>+n*sin(())J  *d?-dn  , 


Ft;X(x\y',z') 


=  (^y)   -Nx«   Jf  gU,n)  -exP[-  ^i«  U-x'  +  n-y')]  -dS-dn,      (10.160) 


and 


MX(X'^   =   Ft,A(x'y) 


I  r  ,   2ll 

g(?,n) --y 

A^ 


:io.i6i) 


,i„(^.V?V)   cos(2l.VIW) 


(-.^7/    m-^)\ 


•d5-dn, 


where 


Ft  x  (x,y) -dx'dy, 


=   tan   (y/x) 


-1 


and   6  =  cos   (z'/R)-    The  incident  spectral  irradiance  at  a  distance   R   from  the  source  and 

the  spectral  radiant  intensity  of  that  source  are  obtained  from  the  spectral  fluence  rate  at 

2 
that  distance  by  multiplying  by   cos6r   and   R  ,   respectively  [eqs.  (10.153)  and  (10.16)],  where 

9    is  the  angle  between  the  normal  to  the  irradiated  surface  at   (x',y')   and  the  direction 
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to  the  origin  in  the  source  plane.   If  the  irradiated  surface  is  parallel  to  the  source  plane 

then   0   =6.    The  spectral  flux  is  obtained  by  integrating  the  spectral  exitance   M, (x,y) 

r  A 

over  the  source  area  to  give: 


=   N, 


source 


g(S,n)  •-, 


{2-n       Jr2_,_  2 


COS 


(2-n       ^(TT~2 


&-V 


(2;  ■  vw 


•d5-.dn' 


;i0.162) 


Now  with  F  ^(x,y)  given  by  eq.  (10.157)  and  g(£,n)  =  V  given  by  eq.  (10.158)  the  evalua- 
tion of  the  integrals  occurring  in  eqs.  (10.159)  to  (10.162)  can  be  carried  out  with  the  help 
of  a  table  of  integrals  [10.15],  [10.16],  [10.17]).   We  find: 


N,   =   2TrA2«a2, 
A         o 


(10.163) 


the  integral  appearing  in  eqs.  (10.159)  and  (10.160)  is 


2TTe-ye2-sin26 


for   sin0<g,   and   0    otherwise, 


(10.164) 


and  the  integral  in  eqs.  (10.161)  and  (10.162)  is 


[1 + ^-gsm 


(10.165) 


Thus  we  finally  obtain  the  following  radiometric  quantities  for  the  source  whose  cross-spectral 
density  is  given  by  eq.  (10.156): 


LA  (x,y,e,<j)) 


A  «exp(-r  /2a  ) 


otherwise 


cosO 


>tt6  -\B2- 


sin 


Ft/A(x',y,,z')   = 


o  *2   2 

2  ttA  •  a  ■ 
o 


.2 


>1T6-R4ye2- 


l+z,2/R2 


0    otherwise 


for   sin0<| 


2   2    2 
1-z'VR  )<B 


for   (1-z'VR  ) 


(10.166) 


(10.167) 


(x,y)   =   A2-exp(-r2/2a2)-%[l  +  i^f-* £n  (jT^ff )] 


(10.168) 


K1(x,ly,(z')   = 


.2 


_  _      z'  «cose 

2ttA  «a  •■ 
o 


2TrB-R4'ye2«l+z,2/R'' 


for   (l-z,2/R2)<62 


0     otherwise 


(10.169) 
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ix(e,<M     = 


2      2 
2irA  »o    ' 
o 


!Trg«^62-si 


in26 


otherwise 


for 


sin6<| 


(10.170) 


$ .       =2  ttA„ 
A  O 


°24  -  ^(tSIt)] 


(10.171) 


When   3=1   the  expression  for  spectral  radiance  [eq.  (10.166)]  becomes  independent  of 
direction: 


Lx(x,y,6,<(>)   =  A2-exp(-r2/2a2)/(2TT). 


(10.166a) 


This,  then,  is  a  Lambertian  source  and  we  see  that  the  degree  of  coherence  of  a  Lambertian 
source,  such  as  a  blackbody,  must  resemble  [from  eq.  (10.158)  with   B  =  1] 


M   =   sin(27rp/X)/(2Trp/X) 


(10.172) 


with  coherence  interval  (if  we  use  the  criterion   |u|=^)   of  the  order  of   0.3A.    This  is 
presumably  about  as  incoherent  as  any  real  source  can  be.   Notice  that  the  coherence  interval 
of  a  Lambertian  source  depends  upon  the  observing  wavelength  rather  than  upon  any  parameter 
characteristic  of  the  source.   The  other  radiometric  quantities  for  a  quasihomogeneous 
Lambertian  source  are  readily  obtained  from  eqs.  (10.167)  to  (10.171)  by  setting   8=1: 


Ft  j^x^y^z') 


,2   2   ,  /r>3 

A  • a    • z ' /R 
o 


(10.167a) 


M, (x,y)   = 


A2-exp(-r2/2o2)/2 


(10.168a) 


EA(x,,y*,z') 


A2 • a2 -cosG  -z'/R3 
o        r 


(10.169a) 


ix(e,<fr) 


2   2 

A  «o  -cosf 
o 


:i0.170a) 


2   2 

>,   =  A  -a  'TT. 
A      o 


(10.171a) 


2      2    2    2     2     2   2  -1 

In  these  expressions   R  =  x'  +y'  +z '  ,   r   =  x  +y  ,   and   6  =  cos   (z'/R). 

Although  in  pursuing  this  example  we  have  made  no  mathematical  approximations  there  is 
an  inconsistency  in  that  we  have  treated  a  source  of  infinite  extent  using  equations  derived 
for  a  finite  source  [(eqs.  (10.149)  to  (10.155)].   Fortunately  the  Gaussian  approaches  zero 
so  fast  as   r   increases  that  the  error  is  negligible  as  long  as   a   is  very  much  smaller 
than  the  observation  distance   R. 


POLARIZATION. 


The  theory  of  partial  coherence  given  in  this  chapter  is  known  as  a  scalar 


theory  because  it  treats  a  single  component  of  the  electric  field  vector  as  though  the  field 
in  other  directions  didn't  exist.   A  more  complete  theory  includes  consideration  of  the 
coherence  between,  say,  the  x-component  of  the  field  at  point   P,   and  the  y-component  of  the 
field  at  point  V~.         Thus,  one  replaces  the  cross-spectral  density  function  with  a  2x2 
matrix  of  four  such  functions: 
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w 


(j'k)(P1,P2)   =  -^-    •<V*j  (Px+AP,t)  -Vk(P2  +  AP,t)>  (10.173! 


where  j  and  k  can  each  be  either  x  or  y  and  V.(P,t)  is  the  complex  representation  of 
the  j  component  of  the  electric  field  at  the  point  P  as  seen  through  a  narrow  (AX) 
pass-band  spectral  filter  or  monochromator .  The  angular  brackets  again  denote  an  average  over 
time  and  over  identical  small  areas  around  P,  and  P-  as  the  small  displacement  AP  covers 
the  measurement  resolution  area.  In  terms  of  these  cross-spectral  density  matrix  elements  the 
Stokes-vector  components  [10.20]  of  spectral  fluence  rate  are: 

Ft,x,o(p)  =  w{x,x) (p'p>  +  wxY'Y) (P'P) 


,.  .....     ..f-»<p.p)  -„«*-•'' 


Ft,X,2(P)  =  W^X'y) (P,P)  +  W^Y'X) (P,P) 


Ft,X,3(P)  =  i-  tWAY'X)  (P'P)  "  wj[x'y)  (P,P)] 


(10.174) 


If   V.(P,t)   is  written  as 
3 


V.(P,t)  =  A.  (P,t)  -exptiy.  (P,t)  -  2Trivt]  (10.1751 


(10.177) 


where   v  =  c/X ,   then 

WA(j'k)(P,P)   =   iy  -<A.(P,t) -Ak(P,t) •exp[iyk(P,t)  -iYj(P,t)]>  (10.176! 

and  the  Stokes  components  of  spectral  fluence  rate  become: 

Ft,X,0(P)  "  XX    *<Ax2  +  Ay2> 

Ft,X,l(P)  "  TX    '<Ax2  "  Ay2> 

F*.  ^  o(p)  =  XT  *<A  -A  -cos(y  -Y  )> 
t,  X,  2       AX     x   y      'x  'y  ' 

F*.    \    ^(p)  =  TT  *<A  -A  -sin(Y  -Y  )>• 
t,X,3       AX   N  x   y      'x  'y  ' 

These  equations  are  equivalent  to  eq.  (6.3)  [10.20].   Thus,  polarization  is  included  as  an 
essential  aspect  of  a  more  general  coherence  theory  with  the  Stokes  parameters  being  one 

possible  set  of  four  linear  combinations  of  the  four  cross-spectral  density  function  matrix 

elements.   Although  there  are  circumstances  involving  polarization  in  which  the  scalar  theory 

of  coherence  is  clearly  inadequate  [10.21],  it  appears  that  a  complete  theory  would  probably 

only  be  required  in  exceptional  radiometric  situations  where  precise  calculations  of  signifi- 
cant coherence  effects  are  needed. 
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SUMMARY  of  CHAPTER  10.    The  field  quantity  which,  in  partially  coherent  fields,  plays  a  basic 
role  somewhat  analogous  to  that  of  spectral  radiance  in  incoherent  fields  is  the  cross- spectral 
density,   W,  (x,  ,  y,  ,  z- ;  x~  ,  y~  ,  z~)  .    It  is  a  complex-valued  function  of  the  coordinates  of  two 
field  points.   The  only  direct  interest  radiometrists  have  in  this  function  is  in  the  special 
case  where  both  sets  of  coordinates  describe  the  same  point.   Then  the  cross-spectral  density 
reduces  to  the  spectral  fluence  rate  at  that  point  (Chapter  4,  page  17  [10.9]): 

WA (x,y,z;x,y,z)   =   Ft  A(x,y,z).  (10.15a) 

Otherwise  its  importance  for  radiometry  is  that  if  one  follows  the  propagation  of  cross- 
spectral  density  from  a  source  through  an  optical  system  with  the  calculation  of  fluence  rate 
or  one  of  the  other  radiometric  quantities  of  geometric  optics  only  at  the  final  detector, 
then  all  coherence  effects  such  as  diffraction  will  be  properly  accounted  for.   The  propagation 
law  for  cross-spectral  density  is: 

C  (•  exp  [2iri-  (s, -s-) /A]  • cos9,  -cose,, 

Vpi'p2J    =       VW  • — — -2 ~ •  dprdp2-  (10-36) 

' '  S, • S~ • A 

source  1   2 

Thus  if  the  cross-spectral  density  function  W  (p   p  )   is  known  for  all  pairs  of  points   P, 
and   P_   in  one  surface  or  cross-section  of  a  beam,  which  we  can  regard  as  a  source,  and  if 
the  optical  path  length   s,   for  a  ray  from  P.   to   P'   and  the  optical  path  length   s_   from 
P_   to   Pi   are  known  for  all  points   P,   and   P_   in  the  source,  then  this  multiple  integral 
over  the  source  surface  permits  the  cross-spectral  density  to  be  calculated  at  the  point  pair 
P',P'.    The  angles   6.   are  the  angles  between  the  normal  to  the  source  surface  at   P.   and 
the  direction  to   P|,   where   i   is  either   1   or   2.    For  a  complicated  optical  system  this 
equation  must  be  applied  sequentially  between  each  successive  pair  of  beam  limiting  apertures 
with  the  results  of  one  calculation  serving  as  the  source  for  the  next  stage  of  the  calcula- 
tion.  The  only  approximation  of  any  consequence  made  in  deriving  this  equation  is  that  the 
source  dimensions  as  well  as  all  distances   s,   and   s_   involved  must  be  much  larger  than  the 
wavelength   A. 

In  the  evaluation  of  the  propagation  integrals  [eq.  (10.36)],  it  is  desirable  to  perform 
as  much  of  the  integration  analytically  as  possible.   To  this  end  it  is  usual  to  restrict  the 
surface  containing   P,   and   P_   to  the   X-Y  plane   (z  =  0)   and  to  expand  the  expressions  for 
s,   and   s2   in  power  series  in   1/R,   and   1/R2  where  R,   is  the  distance  from  the  origin 
to  P.!  ,   and  R_   is  the  distance  from  the  origin  to  P'.    If  we  assume  that  the  source 
dimensions  are  small  compared  to   \/A"r.   and  retain  only  the  most  important  terms  of  these 
expansions  we  obtain  the  result: 

W^x-.y-.z'.-x-.y-.z')   =    4^2   .exp[2£i  .  (R^)]  •  \\\\      V*l'yl,JI2'y2> 


R1R2X  source 


exp 


[„^(v^^.^j^)].dXi.ayi.dX2.dy 


(10.42) 


I 


Although  this  condition  on  the  source  size  is  often  not  fulfilled,  the  approximation 
(Fraunhofer  approximation)  may  be  useful  anyway  in  radiometry  if  the  goal  is  merely  to  esti- 
mate the  magnitude  of  small  measurement  errors  which  diffraction  may  introduce. 
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One  situation  in  which  the  sequential  application  of  the  propagation  equation  (10.36)  is 
not  overly  difficult  is  in  simple  lens  (or  focusing  mirror)  image  formation.   If  the  cross- 
spectral  density  function  W,  (P..  ,P„)   is  known  at  all  point  pairs  in  the  object  plane  then  the 
cross-spectral  density  in  the  image  plane  is  given  by: 


Wx(x«,y«,zJ;x«,y-»,z»)   =  J    4  •  JJjj  «x  (x^y^-x^y^ 

object 


(10.133) 


KX  (xi»Yi'"xi»yi)  *K*  (x2,y2;x^,yp  •dx1»dy1'dx2  «dy2 


In  this  expression  R  is  the  distance  from  the  object  to  the  lens  and  R"   from  the  lens  to 
the  image.   The  function  K, (P,P")   is  the  transmission  function  of  the  optical  system.   It  is 
given  by  an  integral  over  the  lens  aperture  plane: 

oo 

K(x,y;x",y")   =  jj   P  (x'  ,y  ■  )  -exp[^i  •  (s  +  s"  -  R  -  R")]  -dx'  -dy '  .  (10.132) 


Here   s   is  the  optical  path  length  for  a  ray  from  point   (x,y)   in  the  object  plane  to  point 
(x',y')   in  the  lens  aperture  plane  and   s"   is,  similarly,  the  optical  path  length  for  a  ray 
from  the  point   (x',y')   to  the  image  plane  point   (x",y").    Since   s   and   s"   depend  upon 
these  points,  to  be  complete  we  should  write  them  as   s (x,y;x ' ,y ' )   and   s" (x* ,y ' ;x" ,y") . 
The  factor   P(x',y')   is  the  pupil  function  for  the  aperture.   In  its  simplest  form  it  merely 
indicates  where  the  aperture  is  transparent  and  where  it  is  opaque  by  taking  on  the  value   1 
within  the  aperture  and  the  value   0   elsewhere.   In  more  complicated  situations  the  pupil 
function  can  indicate  regions  of  partial  transparency  by  assuming  intermediate  values  between 
0   and   1   and  it  can  indicate  regions  of  phase  retardation  by  assuming  complex  values, 
although  this  latter  condition  can  alternatively  be  accounted  for  in  the  calculations  of   s 
and   s".    For  perfect  imagery  and  a  transparent  rectangular  aperture,   K,(P,P")   reduces  to 

.  [2-na*        /x  .  x"\l     •  r2TTb'    /y  ,  y"\"| 

»    .n         a  -k.    Sinb^  •  (r  +  r")J      sinbr-  •   (r  +  h)\ 

KA(x,y;x",y")   =   4a'b'-  ^ •  ^r (  grr  (10.178) 

X      lR    RV  X      VR  +  R") 

where  the  aperture  width  is   2a'   in  the  x-direction  and  2b1   in  the  y-direction.   Eq. 
(10.178)  substituted  into  eq.  (10.133)  yields  eq.  (10.99),  the  equation  for  perfect  imaging 
with  a  transparent  rectangular  lens  aperture. 

All  of  the  classical  radiometric  quantities  can  be  expressed  in  terms  of  the  cross- 
spectral  density.   Equation  (10.15a)  illustrates  the  simplest  of  these  relationships  —  that 
for  spectral  fluence  rate.   Some  of  the  others  are: 

Spectral  radiance: 

LA(x,y,e,*)  =  °^|^  •  Jj  wA(x+J-2c,  y+hT\;  x-hK ,   y-^n) 

source 


(10.149) 
•  exp  -  — j- —  •  sine  •  (£*cos<{>  +  n*sin<f>)  I  •<!$;•  dn. 
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Spectral  irradiance: 

Ex(x\y,,z')  =   J!  ;  J  Wx(x+Jg5,  Y+hn;   x-HZ,   y-^n) 


source 


(10.151) 


,2 


I      2iTi  I   z 

exp<-  j^-   •[^•(x'-x)  +  ir(y'-y)]>'   4  2  •cos6r'd£;-divdx-dy 


R"V 


Spectral  radiant  intensity: 


ix(e,<j>)  =  S2|_ 


•  exp 


L     x 


w,  (x+^,  y+*5n;  x-JsC  y-h-n) 


sine  •  (£-cos<t>  +  n«sin<|>)  I  •  dx •  dy •  d £ •  dn , 


(10.154) 


Spectral  radiant  exitance: 


MA(x,y)   =     yx(x+HZ,    y+hr];    x-hK,    y-bn) 


2JL 
,2 


sin 


(|i.ViW)    oo.gr.  V?*7) 


(fi.V?V)        (|i.V?T?)_ 


(10.155) 


•dC-dn. 


In  these  equations  the  point   (x,y)   is  a  point  in  a  plane  (which  we  call  a  source  but  may  be 
any  cross-section  of  the  beam)  where  the  cross-spectral  density  is  assumed  completely  known. 
£   and  r\      are  displacements  from  this  point.   The  primed  coordinates  refer  to  a  point,  which 
is  irradiated  by  the  source,  at  a  distance   R   from  the  source  point   (x,y) ,   and   6    is  the 
angle  between  the  normal  to  the  irradiated  surface  at   (x',y')   and  the  direction  to  the  point 
(x,y) .    It  is  possible  to  imagine  coherent  fields  for  which  L,   calculated  from  the  expres- 
sion above  will  not  agree  with  our  classical  ideas  of  acceptable  behavior  for  spectral  radi- 
ance --  by  being  negative,  for  example.   The  reason  for  this  is  that  spectral  radiance  is 
inherently  a  geometrical  optics  concept:   its  definition  is  equivalent  to  requiring  that  the 
position  and  momentum  of  a  photon  be  simultaneously  determinable  with  arbitrary  accuracy, 
which,  of  course,  is  forbidden  by  the  Heisenberg  uncertainty  principle.   Fortunately  such  non- 
classical  behavior  of   L>   should  rarely  occur  in  practice.   Although  the  above  expressions 
for  spectral  radiant  intensity  and  spectral  radiant  exitance  are  both  derived  from  that  for 
spectral  radiance,  only  the  latter  suffers  from  this  same  deficiency. 

The  cross-spectral  density  of  a  uniform  rectangular  incoherent  source  is  usually  taken  to 
be  something  like: 


WA(x1,y1;x2,y2: 


=  \ 


2,2 


Aq   .6(Xl-x2)  •  6(y;L-y2) 


for 


x.  |  ,  I  x_  |  <  a 


and 


\y1 I r I y2 I  <  b 


(10.43; 


and 


WX  (xi'Yi?x2'y2^ 


0    otherwise. 
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This  expression  lends  itself  to  easy  evaluation  of  integrals  such  as  eq.  (10.42)  or 
eq.  (10.133)  and  thereforeiis  widely  used  in  spite  of  some  shortcomings  (e.g.,  it  implies  that 
Ft  A  =  °°^  *    The  °PP°s;i-te  extreme  of  a  uniform  rectangular  coherent  source  has  a  cross-spectral 
density  given  by: 

Wx(x1,y1;x2,y2)   =  A^ 

for   |x1|,|x2|  <  a    and    | y1 | , | y2 |  <  b  (10.59a) 

=   0    otherwise, 

where,  for  simplicity,  uniform  phase  across  the  source  is  assumed   [6   =0   in  eq.  (10.59)]. 
An  idea  of  the  consequences  of  coherence  in  a  radiometric  measurement  may  usually  be  obtained 
by  considering  these  two  source  cross-spectral  density  functions  in  the  sequence  of  expressions 
of  the  form  of  eqs.  (10.42)  and  (10.133)  which  describe  the  measurement  geometry  up  to  the 
detector.   At  the  detector  the  irradiance  (or,  simpler,  and  almost  equivalently,  the  fluence 
rate)  can  be  calculated,  and  if  the  results  for  both  kinds  of  sources  when  considered  in  the 
over-all  measurement  process  (e.g.,  including  calibrations)  show  the  same  behavior  over  the 
range  of  source  sizes  and  distances,  spectral  distributions,  etc.,  to  be  encountered  then  it 
is  safe  to  assume  that  coherence  will  not  affect  these  measurements.   On  the  other  hand  if  the 
results  show  that  the  two  types  of  sources  behave  differently  over  some  sub-range  of  source 
characteristics  then  it  may  be  possible  to  change  the  measurement  geometry  to  avoid  the 
problem  or,  as  a  last  resort,  it  may  be  possible  by  more  elaborate  measurements  and  careful 
calculations  to  correct  the  radiometer  results  for  coherence. 

Another  expression  for  the  cross-spectral  density  function,  which  is  frequently  intro- 
duced as  an  approximate  description  of  many  fields  encountered  in  radiometry  —  especially 
those  originating  in  incandescent  sources,  is 

wA(x+^2c,  y+H;n;  x-HZ,   y-^n)   =  Ffc  x  (x,y)  -gU,n) .  (10.83a) 

This  is  known  as  the  quasihomogeneous  model  and  is  appropriate  for  regions  of  a  field  where 
the  fluence  rate  F,  , (x,y)   can  be  assumed  constant  over  areas  the  size  of  the  coherence 

t,  A 

area.   The  function  g(£^n)   is  simply  the  (complex)  degree  of  coherence  of  the  field.   The 
use  of  this  model  considerably  simplifies  the  propagation  integral  eq.  (10.42)  and,  conse- 
quently, also  the  integrals  for  the  classical  radiometric  quantities  [eqs.  (10.149,  10.151, 
10.154,  and  10.155) ]  . 
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final  reports  on  work  performed  by  NBS  for  outside  sponsors 
(both  government  and  non-government).  In  general,  initial  dis- 
tribution is  handled  by  the  sponsor;  public  distribution  is  by  the 
National  Technical  Information  Service  ,  Springfield,  VA  22161, 
in  paper  copy  or  microfiche  form. 
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